##############ECHO OF PROBLEM################# ##############temp/arctan_sqrtpostode.ode################# diff ( y , x , 1 ) = arctan(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(10.0 * (0.1 * x + 0.2) * arctan(sqrt ( 0.1 * x + 0.2)) - 10.0 * sqrt( 0.1 * x + 0.2) + 10.0 * arctan(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) = 0.6167902654763967 y[1] (numeric) = 0.6167902654763967 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) = 0.6602063571030989 y[1] (numeric) = 0.6602063571030993 absolute error = 3.33066907387547e-16 relative error = 5.044890946657981e-14 % Correct digits = 15 h = 0.005 TOP MAIN SOLVE Loop NO POLE x[1] = 0.3000000000000002 y[1] (analytic) = 0.7044963719892491 y[1] (numeric) = 0.7044963719892502 absolute error = 1.110223024625157e-15 relative error = 1.575910208721556e-13 % Correct digits = 14 h = 0.005 TOP MAIN SOLVE Loop NO POLE x[1] = 0.4000000000000002 y[1] (analytic) = 0.7496341409653187 y[1] (numeric) = 0.7496341409653199 absolute error = 1.221245327087672e-15 relative error = 1.629121808026313e-13 % Correct digits = 14 h = 0.005 TOP MAIN SOLVE Loop NO POLE x[1] = 0.5000000000000003 y[1] (analytic) = 0.795595112510077 y[1] (numeric) = 0.7955951125100781 absolute error = 1.110223024625157e-15 relative error = 1.395462349086634e-13 % Correct digits = 14 h = 0.005 Finished! diff ( y , x , 1 ) = arctan(sqrt(0.1 * x + 0.2)); Iterations = 80 Total Elapsed Time = 6 Seconds Elapsed Time(since restart) = 6 Seconds Time to Timeout = 9 Minutes 54 Seconds Percent Done = 101.3 %