peval.f90 running y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 3 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 polynomial coefficients are: c( 0 )= -0.6500000000000000 c( 1 )= 7.0100000000000007 c( 2 )= -5.2000000000000002 c( 3 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 3 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 3 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 3 ,r( 2 ),c)= 0.0000000000000000E+000 y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 4 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 r( 3 )= 5.0999999999999996 polynomial coefficients are: c( 0 )= 3.3149999999999999 c( 1 )= -36.4009999999999962 c( 2 )= 33.5300000000000011 c( 3 )= -10.3000000000000007 c( 4 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 4 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 4 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 4 ,r( 2 ),c)= 0.0000000000000000E+000 r( 3 )= 5.0999999999999996 , peval( 4 ,r( 3 ),c)= 0.0000000000000000E+000 y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 5 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 r( 3 )= 5.0999999999999996 r( 4 )= 7.6999999999999993 polynomial coefficients are: c( 0 )= -25.5254999999999974 c( 1 )= 2.8360269999999997E+02 c( 2 )= -2.9458199999999999E+02 c( 3 )= 1.1284000000000000E+02 c( 4 )= -18.0000000000000000 c( 5 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 5 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 5 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 5 ,r( 2 ),c)= 0.0000000000000000E+000 r( 3 )= 5.0999999999999996 , peval( 5 ,r( 3 ),c)= 0.0000000000000000E+000 r( 4 )= 7.6999999999999993 , peval( 5 ,r( 4 ),c)= 0.0000000000000000E+000 y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 6 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 r( 3 )= 5.0999999999999996 r( 4 )= 7.6999999999999993 r( 5 )= 12.7999999999999989 polynomial coefficients are: c( 0 )= 3.2672639999999996E+02 c( 1 )= -3.6556400599999997E+03 c( 2 )= 4.0542522999999997E+03 c( 3 )= -1.7389339999999997E+03 c( 4 )= 3.4324000000000001E+02 c( 5 )= -30.7999999999999972 c( 6 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 6 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 6 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 6 ,r( 2 ),c)= 0.0000000000000000E+000 r( 3 )= 5.0999999999999996 , peval( 6 ,r( 3 ),c)= 0.0000000000000000E+000 r( 4 )= 7.6999999999999993 , peval( 6 ,r( 4 ),c)= 0.0000000000000000E+000 r( 5 )= 12.7999999999999989 , peval( 6 ,r( 5 ),c)= 0.0000000000000000E+000 y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 7 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 r( 3 )= 5.0999999999999996 r( 4 )= 7.6999999999999993 r( 5 )= 12.7999999999999989 r( 6 )= 20.5000000000000000 polynomial coefficients are: c( 0 )= -6.6978911999999991E+03 c( 1 )= 7.5267347629999989E+04 c( 2 )= -8.6767812210000004E+04 c( 3 )= 3.9702399299999997E+04 c( 4 )= -8.7753539999999994E+03 c( 5 )= 9.7463999999999999E+02 c( 6 )= -51.2999999999999972 c( 7 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 7 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 7 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 7 ,r( 2 ),c)= 0.0000000000000000E+000 r( 3 )= 5.0999999999999996 , peval( 7 ,r( 3 ),c)= 0.0000000000000000E+000 r( 4 )= 7.6999999999999993 , peval( 7 ,r( 4 ),c)= 0.0000000000000000E+000 r( 5 )= 12.7999999999999989 , peval( 7 ,r( 5 ),c)= 0.0000000000000000E+000 r( 6 )= 20.5000000000000000 , peval( 7 ,r( 6 ),c)= 0.0000000000000000E+000 y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 8 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 r( 3 )= 5.0999999999999996 r( 4 )= 7.6999999999999993 r( 5 )= 12.7999999999999989 r( 6 )= 20.5000000000000000 r( 7 )= 33.2999999999999972 polynomial coefficients are: c( 0 )= 2.2303977695999996E+05 c( 1 )= -2.5131005672789994E+06 c( 2 )= 2.9646354942230000E+06 c( 3 )= -1.4088577089000000E+06 c( 4 )= 3.3192168749999994E+05 c( 5 )= -4.1230865999999995E+04 c( 6 )= 2.6829299999999998E+03 c( 7 )= -84.5999999999999943 c( 8 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 8 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 8 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 8 ,r( 2 ),c)= 0.0000000000000000E+000 r( 3 )= 5.0999999999999996 , peval( 8 ,r( 3 ),c)= 0.0000000000000000E+000 r( 4 )= 7.6999999999999993 , peval( 8 ,r( 4 ),c)= 0.0000000000000000E+000 r( 5 )= 12.7999999999999989 , peval( 8 ,r( 5 ),c)= 0.0000000000000000E+000 r( 6 )= 20.5000000000000000 , peval( 8 ,r( 6 ),c)= 0.0000000000000000E+000 r( 7 )= 33.2999999999999972 , peval( 8 ,r( 7 ),c)= 0.0000000000000000E+000 y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 9 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 r( 3 )= 5.0999999999999996 r( 4 )= 7.6999999999999993 r( 5 )= 12.7999999999999989 r( 6 )= 20.5000000000000000 r( 7 )= 33.2999999999999972 r( 8 )= 53.7999999999999972 polynomial coefficients are: c( 0 )= -1.1999540000447998E+07 c( 1 )= 1.3542785029657015E+08 c( 2 )= -1.6201049015647641E+08 c( 3 )= 7.8761180233043000E+07 c( 4 )= -1.9266244496399999E+07 c( 5 )= 2.5501422782999994E+06 c( 6 )= -1.8557250000000000E+05 c( 7 )= 7.2344099999999999E+03 c( 8 )= -1.3839999999999998E+02 c( 9 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 9 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 9 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 9 ,r( 2 ),c)= 0.0000000000000000E+000 r( 3 )= 5.0999999999999996 , peval( 9 ,r( 3 ),c)= 0.0000000000000000E+000 r( 4 )= 7.6999999999999993 , peval( 9 ,r( 4 ),c)= 0.0000000000000000E+000 r( 5 )= 12.7999999999999989 , peval( 9 ,r( 5 ),c)= 0.0000000000000000E+000 r( 6 )= 20.5000000000000000 , peval( 9 ,r( 6 ),c)= 0.0000000000000000E+000 r( 7 )= 33.2999999999999972 , peval( 9 ,r( 7 ),c)= 0.0000000000000000E+000 r( 8 )= 53.7999999999999972 , peval( 9 ,r( 8 ),c)= 0.0000000000000000E+000 y = c(0) + c(1)*x + c(2)*x^2 +...+ c(n)*x^n, for n 10 roots are: r(0)= 0.1000000000000000 r(1)= 2.5000000000000000 r( 2 )= 2.6000000000000001 r( 3 )= 5.0999999999999996 r( 4 )= 7.6999999999999993 r( 5 )= 12.7999999999999989 r( 6 )= 20.5000000000000000 r( 7 )= 33.2999999999999972 r( 8 )= 53.7999999999999972 r( 9 )= 87.0999999999999943 polynomial coefficients are: c( 0 )= 1.0451599340390205E+09 c( 1 )= -1.1807765300831707E+10 c( 2 )= 1.4246541542925665E+10 c( 3 )= -7.0221092884545212E+09 c( 4 )= 1.7568510758694828E+09 c( 5 )= -2.4138363693632993E+08 c( 6 )= 1.8713507028299998E+07 c( 7 )= -8.1568961099999992E+05 c( 8 )= 1.9289049999999996E+04 c( 9 )= -2.2549999999999997E+02 c( 10 )= 1.0000000000000000 polynomial evaluated at each root: r( 0 )= 0.1000000000000000 , peval( 10 ,r( 0 ),c)= 0.0000000000000000E+000 r( 1 )= 2.5000000000000000 , peval( 10 ,r( 1 ),c)= 0.0000000000000000E+000 r( 2 )= 2.6000000000000001 , peval( 10 ,r( 2 ),c)= 0.0000000000000000E+000 r( 3 )= 5.0999999999999996 , peval( 10 ,r( 3 ),c)= 0.0000000000000000E+000 r( 4 )= 7.6999999999999993 , peval( 10 ,r( 4 ),c)= 0.0000000000000000E+000 r( 5 )= 12.7999999999999989 , peval( 10 ,r( 5 ),c)= 0.0000000000000000E+000 r( 6 )= 20.5000000000000000 , peval( 10 ,r( 6 ),c)= 0.0000000000000000E+000 r( 7 )= 33.2999999999999972 , peval( 10 ,r( 7 ),c)= 0.0000000000000000E+000 r( 8 )= 53.7999999999999972 , peval( 10 ,r( 8 ),c)= 0.0000000000000000E+000 r( 9 )= 87.0999999999999943 , peval( 10 ,r( 9 ),c)= 0.0000000000000000E+000