##############ECHO OF PROBLEM################# ##############temp/mtest6postode.ode################# diff (x1,t,1) = 4.0 * x2 - 2.0 * diff (x2,t ,1) - 2.0 * x1; diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1) + x1; ! // BEGIN FIRST INPUT BLOCK Digits=32; max_terms=30; ! // END FIRST INPUT BLOCK // BEGIN SECOND INPUT BLOCK t_start=1.5; // # did poorly with t_start := 0.5; t_end=5.0; array_x1_init[0 + 1] = exact_soln_x1(t_start); array_x1_init[1 + 1] = exact_soln_x1p(t_start); array_x2_init[0 + 1] = exact_soln_x2(t_start); array_x2_init[1 + 1] = exact_soln_x2p(t_start); glob_h=0.00001; glob_look_poles=true; glob_max_iter=100; // END SECOND INPUT BLOCK // BEGIN OVERRIDE BLOCK glob_desired_digits_correct=10; glob_display_interval=0.001; glob_look_poles=true; glob_max_iter=10000000; glob_max_minutes=3; // END OVERRIDE BLOCK ! // BEGIN USER DEF BLOCK double exact_soln_x1 (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(2.0 * c1 + 6.0 * c3 * exp(-t)); } double exact_soln_x1p (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( - 6.0 * c3 * exp(-t)); } double exact_soln_x2 (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(c1 + c2 * exp(2.0 * t) + c3 * exp(-t)); } double exact_soln_x2p (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( 2.0 * c2 * exp(2.0 * t) - c3 * exp(-t)); } // END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Optimize min_size = 0 min_size = 1 opt_iter = 1 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 3500 step_error = 2.857142857142858e-14 est_needed_step_err = 2.857142857142858e-14 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 1.193323441522736e+50 max_value3 = 1.193323441522736e+50 value3 = 1.193323441522736e+50 best_h = 0 opt_iter = 2 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 7000 step_error = 1.428571428571429e-14 est_needed_step_err = 1.428571428571429e-14 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 4.445774813306187e+41 max_value3 = 4.445774813306187e+41 value3 = 4.445774813306187e+41 best_h = 0 opt_iter = 3 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 14000 step_error = 7.142857142857144e-15 est_needed_step_err = 7.142857142857144e-15 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 1.656402686869081e+33 max_value3 = 1.656402686869081e+33 value3 = 1.656402686869081e+33 best_h = 0 opt_iter = 4 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 28000 step_error = 3.571428571428572e-15 est_needed_step_err = 3.571428571428572e-15 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 6.172239484265773e+24 max_value3 = 6.172239484265773e+24 value3 = 6.172239484265773e+24 best_h = 0 opt_iter = 5 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 56000 step_error = 1.785714285714286e-15 est_needed_step_err = 1.785714285714286e-15 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 2.30057535835239e+16 max_value3 = 2.30057535835239e+16 value3 = 2.30057535835239e+16 best_h = 0 opt_iter = 6 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 112000 step_error = 8.92857142857143e-16 est_needed_step_err = 8.92857142857143e-16 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 85795450.95080923 max_value3 = 85795450.95080923 value3 = 85795450.95080923 best_h = 0 opt_iter = 7 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 224000 step_error = 4.464285714285715e-16 est_needed_step_err = 4.464285714285715e-16 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 0.3203036139898109 max_value3 = 0.3203036139898109 value3 = 0.3203036139898109 best_h = 0 opt_iter = 8 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 448000 step_error = 2.232142857142857e-16 est_needed_step_err = 2.232142857142857e-16 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 1.198409406646213e-09 max_value3 = 1.198409406646213e-09 value3 = 1.198409406646213e-09 best_h = 0 opt_iter = 9 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 896000 step_error = 1.116071428571429e-16 est_needed_step_err = 1.116071428571429e-16 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 4.503649142655383e-18 max_value3 = 4.503649142655383e-18 value3 = 4.503649142655383e-18 best_h = 3.90625e-06 START of Soultion t[1] = 1.5 x1[1] (analytic) = 2.000401634288267 x1[1] (numeric) = 2.000401634288267 absolute error = 0 relative error = 0 % Correct digits = 16 h = 3.90625e-06 x2[1] (analytic) = 1.004084046432682 x2[1] (numeric) = 1.004084046432682 absolute error = 0 relative error = 0 % Correct digits = 16 h = 3.90625e-06 t[1] = 1.5 x1[1] (analytic) = 2.000401634288267 x1[1] (numeric) = 2.000401634288267 absolute error = 0 relative error = 0 % Correct digits = 16 h = 3.90625e-06 x2[1] (analytic) = 1.004084046432682 x2[1] (numeric) = 1.004084046432682 absolute error = 0 relative error = 0 % Correct digits = 16 h = 3.90625e-06 TOP MAIN SOLVE Loop NO POLE NO POLE Finished! Maximum Time Reached before Solution Completed! diff (x1,t,1) = 4.0 * x2 - 2.0 * diff (x2,t ,1) - 2.0 * x1; diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1) + x1; Iterations = 129 Total Elapsed Time = 3 Minutes 1 Seconds Elapsed Time(since restart) = 2 Minutes 48 Seconds Expected Time Remaining = 14 Days 10 Hours 28 Minutes 46 Seconds Optimized Time Remaining = 13 Days 9 Hours 35 Minutes 39 Seconds Expected Total Time = 13 Days 9 Hours 38 Minutes 40 Seconds Time to Timeout Unknown Percent Done = 0.01451 %