##############ECHO OF PROBLEM################# ##############temp/sing3postode.ode################# diff ( y , x , 1 ) = m1 * 2.0 / x / x / x ; ! // BEGIN FIRST INPUT BLOCK Digits = 32; max_terms = 30; ! // END FIRST INPUT BLOCK // BEGIN SECOND INPUT BLOCK x_start = -1.0; x_end = 0.7; array_y_init[0 + 1] = exact_soln_y(x_start); glob_h = 0.0001; glob_look_poles = true; glob_max_iter = 100; // END SECOND INPUT BLOCK // BEGIN OVERRIDE BLOCK glob_h = 0.005 ; glob_display_interval = 0.1; glob_look_poles = true; glob_max_iter = 10000; glob_max_minutes = 10; // END OVERRIDE BLOCK ! // BEGIN USER DEF BLOCK double exact_soln_y (double x) { return(1.0/x/x); } // END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion x[1] = -1 y[1] (analytic) = 1 y[1] (numeric) = 1 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.005 TOP MAIN SOLVE Loop Complex estimate of poles used Radius of convergence = 0.5295 Order of pole = 0.6402 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = -0.8999999999999999 y[1] (analytic) = 1.234567901234568 y[1] (numeric) = 1.234567901317548 absolute error = 8.29798452173236e-11 relative error = 6.721367462603211e-09 % Correct digits = 10 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.9014 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = -0.7999999999999998 y[1] (analytic) = 1.562500000000001 y[1] (numeric) = 1.562500000286313 absolute error = 2.863118631069028e-10 relative error = 1.832395923884177e-08 % Correct digits = 9 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.8012 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = -0.6999999999999997 y[1] (analytic) = 2.040816326530614 y[1] (numeric) = 2.040816327375615 absolute error = 8.450018462724529e-10 relative error = 4.140509046735017e-08 % Correct digits = 9 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.7011 Order of pole = 2.082 Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used x[1] = -0.5999999999999996 y[1] (analytic) = 2.777777777777781 y[1] (numeric) = 2.777777780403758 absolute error = 2.625976325276724e-09 relative error = 9.453514770996196e-08 % Correct digits = 9 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.6009 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = -0.4999999999999996 y[1] (analytic) = 4.000000000000007 y[1] (numeric) = 4.000000009566065 absolute error = 9.566058167820302e-09 relative error = 2.391514541955071e-07 % Correct digits = 8 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.5008 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = -0.3999999999999995 y[1] (analytic) = 6.250000000000017 y[1] (numeric) = 6.250000045614327 absolute error = 4.561430966276703e-08 relative error = 7.298289546042704e-07 % Correct digits = 8 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.4006 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = -0.2999999999999994 y[1] (analytic) = 11.11111111111116 y[1] (numeric) = 11.11111144975292 absolute error = 3.386417599671177e-07 relative error = 3.047775839704046e-06 % Correct digits = 7 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.3005 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = -0.1999999999999993 y[1] (analytic) = 25.00000000000018 y[1] (numeric) = 25.00000566617634 absolute error = 5.666176164709213e-06 relative error = 2.266470465883669e-05 % Correct digits = 6 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.2003 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Complex estimate of poles used Real estimate of pole used x[1] = -0.0999999999999992 y[1] (analytic) = 100.0000000000016 y[1] (numeric) = 100.0006799675768 absolute error = 0.0006799675751665291 relative error = 0.0006799675751665182 % Correct digits = 5 h = 0.005 TOP MAIN SOLVE Loop Real estimate of pole used Radius of convergence = 0.1002 Order of pole = 2.082 Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used Real estimate of pole used x[1] = 8.205242041370298e-16 y[1] (analytic) = 1.48531034844043e+30 y[1] (numeric) = 470105018.7434118 absolute error = 1.48531034844043e+30 relative error = 100 % Correct digits = 0 h = 0.005 TOP MAIN SOLVE Loop NO POLE Complex estimate of poles used x[1] = 0.1000000000000008 y[1] (analytic) = 99.99999999999832 y[1] (numeric) = nan absolute error = nan relative error = nan % Correct digits = -2147483648 h = 0.005 TOP MAIN SOLVE Loop NO POLE Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used x[1] = 0.2000000000000009 y[1] (analytic) = 24.99999999999977 y[1] (numeric) = nan absolute error = nan relative error = nan % Correct digits = -2147483648 h = 0.005 TOP MAIN SOLVE Loop NO POLE Complex estimate of poles used Complex estimate of poles used x[1] = 0.300000000000001 y[1] (analytic) = 11.11111111111104 y[1] (numeric) = nan absolute error = nan relative error = nan % Correct digits = -2147483648 h = 0.005 TOP MAIN SOLVE Loop NO POLE Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used x[1] = 0.4000000000000011 y[1] (analytic) = 6.249999999999966 y[1] (numeric) = nan absolute error = nan relative error = nan % Correct digits = -2147483648 h = 0.005 TOP MAIN SOLVE Loop NO POLE Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used x[1] = 0.5000000000000011 y[1] (analytic) = 3.999999999999982 y[1] (numeric) = nan absolute error = nan relative error = nan % Correct digits = -2147483648 h = 0.005 TOP MAIN SOLVE Loop NO POLE Complex estimate of poles used x[1] = 0.6000000000000012 y[1] (analytic) = 2.777777777777767 y[1] (numeric) = nan absolute error = nan relative error = nan % Correct digits = -2147483648 h = 0.005 TOP MAIN SOLVE Loop NO POLE Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used Complex estimate of poles used x[1] = 0.7000000000000013 y[1] (analytic) = 2.040816326530605 y[1] (numeric) = nan absolute error = nan relative error = nan % Correct digits = -2147483648 h = 0.005 Finished! diff ( y , x , 1 ) = m1 * 2.0 / x / x / x ; Iterations = 340 Total Elapsed Time = 19 Seconds Elapsed Time(since restart) = 19 Seconds Time to Timeout = 9 Minutes 41 Seconds Percent Done = 100.3 %