##############ECHO OF PROBLEM################# ##############temp/sing5postode.ode################# diff ( y , x , 1 ) = m1 * 3.0 / x / 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.00001; 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/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 Real estimate of pole used Radius of convergence = 1.002 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.8999999999999999 y[1] (analytic) = -1.371742112482854 y[1] (numeric) = -1.371742112788488 absolute error = 3.056339625828741e-10 relative error = 2.228071587229151e-08 % Correct digits = 9 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 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.7999999999999998 y[1] (analytic) = -1.953125000000001 y[1] (numeric) = -1.953125001144094 absolute error = 1.144092598437396e-09 relative error = 5.857754103999465e-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 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 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.915451895043735 y[1] (numeric) = -2.915451898805259 absolute error = 3.761523537093581e-09 relative error = 1.290202573223097e-07 % Correct digits = 8 h = 0.005 Finished! diff ( y , x , 1 ) = m1 * 3.0 / x / x / x / x ; Iterations = 60 Total Elapsed Time = 5 Seconds Elapsed Time(since restart) = 4 Seconds Time to Timeout = 9 Minutes 55 Seconds Percent Done = 101.7 %