/************************************************************** This program computes x such that f(x) = 0 , where f(x) is an N-th order polynominal equation with coefficients A[i] for (i = 0 to N); Guess x[k] for k = 0 and then the program computes recursively x[k+1] = x[k] - h*f(x[k]) until you stop the program by typing-in "s" or you type "g" for the next guess for x[0] ; ************************************************************** Type "s" to stop program. Type "g" to start the next guess. **************************************************************/ N=5 A[0]=-14350.000000 A[1]=14895.000000 A[2]=-5580.000000 A[3]=1060.000000 A[4]=-110.000000 A[5]=5.000000 h=0.200000 x[0]=5.000000 k x[k] fx[k] gx[k] 0 5.000000 0.000000 5.000000 1 5.000000 0.000000 5.000000 2 5.000000 0.000000 5.000000 ******************** Next Guess ****************** h = 0.200000 x[0] = 4.999900 k x[k] fx[k] gx[k] 0 4.999900 0.078001 4.984300 1 4.984300 12.262939 2.531712 2 2.531712 796.245885 -156.717465 3 -156.717465 -543241501456.998962 108648300134.682327 ******************** Next Guess ****************** h = -0.200000 x[0] = 4.999000 k x[k] fx[k] gx[k] 0 4.999000 0.780070 5.155014 1 5.155014 -118.810004 -18.606987 2 -18.606987 -33389489.130194 -6677916.433026 3 -6677916.433026 -66401262260180154620000000000000000.00000 -13280252452036031380000000000000000.00000 ******************** Next Guess ****************** h = 0.200000 x[0] = 5.000000 k x[k] fx[k] gx[k] 0 5.000000 -0.000000 5.000000 1 5.000000 -0.000012 5.000002 2 5.000002 -0.001922 5.000387 3 5.000387 -0.301764 5.060740 4 5.060740 -47.093808 14.479501 5 14.479501 596457.477633 -119277.016026 6 -119277.016026 -120735208232456078400000000.000000 24147041646491217390000000.000000 ******************** Next Guess ****************** h = 0.200000 x[0] = 5.000000 k x[k] fx[k] gx[k] 0 5.000000 0.000000 5.000000 1 5.000000 0.000000 5.000000 2 5.000000 0.000002 5.000000 3 5.000000 0.000247 4.999950 4 4.999950 0.038772 4.992196 5 4.992196 6.091405 3.773915 6 3.773915 878.873286 -172.000742 7 -172.000742 -854536046838.313843 170907209195.662048 8 170907209195.662048 7.290745134752054728000000000000000000000e+56 -1.458149026950411076000000000000000000000e+56 ******************** Next Guess ****************** h = 0.200000 x[0] = 7.000000 k x[k] fx[k] gx[k] 0 7.000000 0.000000 7.000000 1 7.000000 0.000000 7.000000 2 7.000000 0.000000 7.000000 3 7.000000 0.000000 7.000000 ******************** Next Guess ****************** h = 0.200000 x[0] = 2.100000 k x[k] fx[k] gx[k] 0 2.100000 203.274050 -38.554810 1 -38.554810 -738642962.805205 147728554.006231 2 147728554.006231 3.517970344322391444000000000000000000000e+41 -7.035940688644783276000000000000000000000e+40