##############ECHO OF PROBLEM################# ##############temp/sqrt_sqrtpostode.ode################# diff ( y , x , 1 ) = sqrt(sqrt(0.1 * x + 0.2)); ! /* BEGIN FIRST INPUT BLOCK */ Digits = 32; max_terms = 30; ! /* END FIRST INPUT BLOCK */ /* BEGIN SECOND INPUT BLOCK */ x_start = 0.1; x_end = 0.5 ; array_y_init[0 + 1] = exact_soln_y(x_start); glob_h = 0.05; glob_look_poles = true; glob_max_iter = 1000000; /* 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(0.8 * (x + 2.0) * sqrt(sqrt(0.1 * x + 0.2))); } /* END USER DEF BLOCK */ #######END OF ECHO OF PROBLEM################# START of Soultion x[1] = 0.1 y[1] (analytic) = 1.137271367855683 y[1] (numeric) = 1.137271367855683 absolute error = 0 relative error = 0 % Correct digits = 16 h = 0.005 TOP MAIN SOLVE Loop NO POLE x[1] = 0.2000000000000001 y[1] (analytic) = 1.205364337353195 y[1] (numeric) = 1.205364337353195 absolute error = 8.881784197001252e-16 relative error = 7.368547352665471e-14 % Correct digits = 15 h = 0.005 TOP MAIN SOLVE Loop NO POLE x[1] = 0.3000000000000002 y[1] (analytic) = 1.274235739780028 y[1] (numeric) = 1.274235739780029 absolute error = 6.661338147750939e-16 relative error = 5.227712533711297e-14 % Correct digits = 15 h = 0.005 TOP MAIN SOLVE Loop NO POLE x[1] = 0.4000000000000002 y[1] (analytic) = 1.343860036446944 y[1] (numeric) = 1.343860036446945 absolute error = 1.110223024625157e-15 relative error = 8.261448324339604e-14 % Correct digits = 15 h = 0.005 TOP MAIN SOLVE Loop NO POLE x[1] = 0.5000000000000003 y[1] (analytic) = 1.414213562373095 y[1] (numeric) = 1.414213562373097 absolute error = 1.332267629550188e-15 relative error = 9.420554752102648e-14 % Correct digits = 15 h = 0.005 Finished! diff ( y , x , 1 ) = sqrt(sqrt(0.1 * x + 0.2)); Iterations = 80 Total Elapsed Time = 9 Seconds Elapsed Time(since restart) = 9 Seconds Time to Timeout = 9 Minutes 51 Seconds Percent Done = 101.3 %