##############ECHO OF PROBLEM################# ##############temp/h3sinpostode.ode################# diff ( y , x , 3 ) = sin(x); ! /* BEGIN FIRST INPUT BLOCK */ Digits=32; max_terms=30; ! /* END FIRST INPUT BLOCK */ /* BEGIN SECOND INPUT BLOCK */ x_start=0.1; x_end=1.5; /* # */ /* # trouble about Pi/2 */ /* # */ array_y_init[0 + 1] = exact_soln_y(x_start); array_y_init[1 + 1] = exact_soln_yp(x_start); array_y_init[2 + 1] = exact_soln_ypp(x_start); glob_look_poles=true; glob_max_iter=20; /* END SECOND INPUT BLOCK */ /* BEGIN OVERRIDE BLOCK */ glob_desired_digits_correct=10; glob_display_interval=0.1; glob_look_poles=true; glob_max_iter=10000000; glob_max_minutes=3; glob_subiter_method=3; /* END OVERRIDE BLOCK */ ! /* BEGIN USER DEF BLOCK */ double exact_soln_y (double x) { return(1.0 + cos(x)); } double exact_soln_yp (double x) { return(-sin(x)); } double exact_soln_ypp (double x) { return(-cos(x)); } /* END USER DEF BLOCK */ #######END OF ECHO OF PROBLEM################# START of Optimize min_size = 0 min_size = 1 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 1.4 estimated_steps = 1400000 step_error = 7.142857142857143e-17 est_needed_step_err = 7.142857142857143e-17 opt_iter = 1 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 2.467208605532302e-183 estimated_step_error = 2.467208605532302e-183 best_h = 2e-06 opt_iter = 2 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 1.655715664606566e-175 estimated_step_error = 1.655715664606566e-175 best_h = 4e-06 opt_iter = 3 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 1.111131969458436e-167 estimated_step_error = 1.111131969458436e-167 best_h = 8e-06 opt_iter = 4 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 7.456680367024187e-160 estimated_step_error = 7.456680367024187e-160 best_h = 1.6e-05 opt_iter = 5 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 5.004093412037817e-152 estimated_step_error = 5.004093412037817e-152 best_h = 3.2e-05 opt_iter = 6 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 3.358190142481625e-144 estimated_step_error = 3.358190142481625e-144 best_h = 6.4e-05 opt_iter = 7 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 2.25364312158092e-136 estimated_step_error = 2.25364312158092e-136 best_h = 0.000128 opt_iter = 8 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 1.512394117654307e-128 estimated_step_error = 1.512394117654307e-128 best_h = 0.000256 opt_iter = 9 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 1.014950270158056e-120 estimated_step_error = 1.014950270158056e-120 best_h = 0.000512 opt_iter = 10 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 6.811212724408108e-113 estimated_step_error = 6.811212724408108e-113 best_h = 0.001024 opt_iter = 11 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 4.570923134360026e-105 estimated_step_error = 4.570923134360026e-105 best_h = 0.002048 opt_iter = 12 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 3.067488750234022e-97 estimated_step_error = 3.067488750234022e-97 best_h = 0.004096 opt_iter = 13 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 2.058549011618091e-89 estimated_step_error = 2.058549011618091e-89 best_h = 0.008192 opt_iter = 14 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 1.381458319740245e-81 estimated_step_error = 1.381458319740245e-81 best_h = 0.016384 opt_iter = 15 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 9.270668118428106e-74 estimated_step_error = 9.270668118428106e-74 best_h = 0.032768 opt_iter = 16 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 6.221249002016116e-66 estimated_step_error = 6.221249002016116e-66 best_h = 0.065536 opt_iter = 17 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 4.174750875079242e-58 estimated_step_error = 4.174750875079242e-58 best_h = 0.131072 opt_iter = 18 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 2.801274754399603e-50 estimated_step_error = 2.801274754399603e-50 best_h = 0.1 START of Soultion TOP MAIN SOLVE Loop x[1] = 0.1 y[1] (analytic) = 1.995004165278026 y[1] (numeric) = 1.995004165278026 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.09933 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.2 y[1] (analytic) = 1.980066577841242 y[1] (numeric) = 1.980066577841242 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.1947 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.3 y[1] (analytic) = 1.955336489125606 y[1] (numeric) = 1.955336489125606 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.2823 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.4 y[1] (analytic) = 1.921060994002885 y[1] (numeric) = 1.921060994002885 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.3587 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.5 y[1] (analytic) = 1.877582561890373 y[1] (numeric) = 1.877582561890373 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.4207 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.6 y[1] (analytic) = 1.825335614909678 y[1] (numeric) = 1.825335614909678 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.466 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.7 y[1] (analytic) = 1.764842187284489 y[1] (numeric) = 1.764842187284488 absolute error = 2.220446049250313e-16 relative error = 1.258155581982573e-14 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.4927 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.7999999999999999 y[1] (analytic) = 1.696706709347165 y[1] (numeric) = 1.696706709347165 absolute error = 2.220446049250313e-16 relative error = 1.308679948642782e-14 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.4998 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.8999999999999999 y[1] (analytic) = 1.621609968270664 y[1] (numeric) = 1.621609968270664 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.4869 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 0.9999999999999999 y[1] (analytic) = 1.54030230586814 y[1] (numeric) = 1.54030230586814 absolute error = 2.220446049250313e-16 relative error = 1.44156510107854e-14 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.4546 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 1.1 y[1] (analytic) = 1.453596121425577 y[1] (numeric) = 1.453596121425577 absolute error = 2.220446049250313e-16 relative error = 1.527553641979085e-14 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.4042 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 1.2 y[1] (analytic) = 1.362357754476674 y[1] (numeric) = 1.362357754476673 absolute error = 2.220446049250313e-16 relative error = 1.629855331284299e-14 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.3377 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 1.3 y[1] (analytic) = 1.267498828624587 y[1] (numeric) = 1.267498828624587 absolute error = 2.220446049250313e-16 relative error = 1.751832821541781e-14 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.2578 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] = 1.4 y[1] (analytic) = 1.169967142900241 y[1] (numeric) = 1.169967142900241 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.1675 Order of pole (three term test) = 28 NO COMPLEX POLE (six term test) for Equation 1 Finished! diff ( y , x , 3 ) = sin(x); Iterations = 14 Total Elapsed Time = 0.0 Seconds Elapsed Time(since restart) = 0.0 Seconds Time to Timeout = 3 Minutes 0.0 Seconds Percent Done = 107.1 %