/************************************************************** This program computes z such that f(z) = 0 , where f(z) is an N-th order polynominal equation with coefficients C[i] = A[i] + j*B[i] for (i = 0 to N) and g(z) = df(z)/dz ; Guess z[k] = x[k] + j*y[k] for k = 0 and then the program computes recursively z[k+1] = z[k] - f(z[k]) /g(z) until you stop the program by typing-in "s" or you type "g" for the next guess for z[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 x[0]=5.000000 k x[k] fx[k] ffx[k] gx[k] 0 5.000000 0.000000 -780.000000 5.000000 1 5.000000 0.000000 -780.000000 5.000000 2 5.000000 0.000000 -780.000000 5.000000 3 5.000000 0.000000 -780.000000 5.000000 ******************** Next Guess ****************** h = 1.000000 x[0] = 4.000000 k x[k] fx[k] ffx[k] gx[k] 0 4.000000 750.000000 -625.000000 5.200000 1 5.200000 -152.294400 -738.280000 4.993717 2 4.993717 4.903251 -780.866568 4.999997 3 4.999997 0.002706 -780.000486 5.000000 4 5.000000 0.000000 -780.000000 5.000000 5 5.000000 -0.000000 -780.000000 5.000000 ******************** Next Guess ****************** h = 1.000000 x[0] = 1.000000 k x[k] fx[k] ffx[k] gx[k] 0 1.000000 -4080.000000 6500.000000 1.627692 1 1.627692 -1032.943361 3433.023247 1.928577 2 1.928577 -162.933410 2389.459376 1.996765 3 1.996765 -7.051043 2184.457638 1.999993 4 1.999993 -0.015237 2175.020456 2.000000 5 2.000000 -0.000000 2175.000000 2.000000 6 2.000000 -0.000000 2175.000000 2.000000 ******************** Next Guess ****************** h = 1.000000 x[0] = -1.000000 k x[k] fx[k] ffx[k] gx[k] 0 -1.000000 -36000.000000 29700.000000 0.212121 1 0.212121 -11431.632333 12666.663715 1.114619 2 1.114619 -3373.527782 5835.893971 1.692684 3 1.692684 -817.877184 3187.216568 1.949296 4 1.949296 -114.085229 2326.018719 1.998343 5 1.998343 -3.607610 2179.841068 1.999998 6 1.999998 -0.004002 2175.005373 2.000000 7 2.000000 -0.000000 2175.000000 2.000000 8 2.000000 0.000000 2175.000000 2.000000 ******************** Next Guess ****************** h = 1.000000 x[0] = 3.000000 k x[k] fx[k] ffx[k] gx[k] 0 3.000000 1040.000000 180.000000 -2.777778 1 -2.777778 -128876.058866 81351.199512 -1.193584 2 -1.193584 -42115.768422 33544.692498 0.061928 3 0.061928 -13448.731059 14215.974194 1.007958 4 1.007958 -4028.466023 6452.222502 1.632311 5 1.632311 -1017.129185 3415.162414 1.930138 6 1.930138 -159.205996 2384.642583 1.996901 7 1.996901 -6.753638 2184.059071 1.999994 8 1.999994 -0.013983 2175.018773 2.000000 9 2.000000 -0.000000 2175.000000 2.000000 ******************** Next Guess ****************** h = 1.000000 x[0] = 4.000000 k x[k] fx[k] ffx[k] gx[k] 0 4.000000 750.000000 -625.000000 5.200000 1 5.200000 -152.294400 -738.280000 4.993717 2 4.993717 4.903251 -780.866568 4.999997 3 4.999997 0.002706 -780.000486 5.000000 4 5.000000 0.000000 -780.000000 5.000000 5 5.000000 -0.000000 -780.000000 5.000000 ******************** Next Guess ****************** h = 1.000000 x[0] = 6.000000 k x[k] fx[k] ffx[k] gx[k] 0 6.000000 -580.000000 -225.000000 3.422222 1 3.422222 1016.944999 -260.128499 7.331617 2 7.331617 744.091309 2839.863798 7.069600 3 7.069600 125.685208 1913.746701 7.003926 4 7.003926 6.696381 1711.717962 7.000013 5 7.000013 0.022855 1700.040064 7.000000 6 7.000000 0.000000 1700.000000 7.000000 7 7.000000 0.000000 1700.000000 7.000000 8 7.000000 0.000000 1700.000000 7.000000 ******************** Next Guess ****************** h = 1.000000 x[0] = 8.000000 k x[k] fx[k] ffx[k] gx[k] 0 8.000000 3690.000000 6255.000000 7.410072 1 7.410072 979.226856 3157.575327 7.099952 2 7.099952 185.240160 2011.006929 7.007839 3 7.007839 13.417951 1723.439412 7.000053 4 7.000053 0.090727 1700.159036 7.000000 5 7.000000 0.000004 1700.000007 7.000000 6 7.000000 0.000000 1700.000000 7.000000 7 7.000000 0.000000 1700.000000 7.000000 8 7.000000 0.000000 1700.000000 7.000000