##############ECHO OF PROBLEM################# ##############temp/expt_sin_sin_newpostode.ode################# diff ( y , x , 1 ) = expt ( sin ( 2.0 * x + 3.0 ) , sin ( x + 1.0 ) ) * ( ( cos ( x + 1.0 ) * ln ( sin ( 2.0 * x + 3.0 ) ) ) + ( ( 2.0 * sin ( x + 1.0 ) * cos ( 2.0 * x + 3.0 ) ) / sin ( 2.0 * x + 3.0 ) ) ) ; ! /* BEGIN FIRST INPUT BLOCK */ Digits=32; max_terms=40; ! /* END FIRST INPUT BLOCK */ /* BEGIN SECOND INPUT BLOCK */ x_start=1.8; x_end=1.81; array_y_init[0 + 1] = exact_soln_y(x_start); glob_look_poles=true; /* END SECOND INPUT BLOCK */ /* BEGIN OVERRIDE BLOCK */ glob_desired_digits_correct=10; glob_display_interval=0.01; glob_look_poles=true; glob_max_iter=100000000; glob_max_minutes=10.0; glob_subiter_method=3; /* END OVERRIDE BLOCK */ ! /* BEGIN USER DEF BLOCK */ double exact_soln_y (double x) { return(expt(sin(2.0*x+3.0),sin(x+1.0))); } /* END USER DEF BLOCK */ #######END OF ECHO OF PROBLEM################# START of Optimize min_size = 0 min_size = 1 glob_desired_digits_correct = 10 estimated_h = 1e-06 estimated_answer = 1 desired_abs_gbl_error = 1e-10 range = 0.01000000000000001 estimated_steps = 10000.00000000001 step_error = 2.499999999999999e-14 est_needed_step_err = 2.499999999999999e-14 opt_iter = 1 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 NO REAL POLE (three term test) for Equation 1 NO COMPLEX POLE (six term test) for Equation 1 min_size = 0.6766027626496453 min_size = 1 glob_desired_digits_correct = 10 estimated_h = 1e-06 estimated_answer = 1 desired_abs_gbl_error = 1e-10 range = 0.01000000000000001 estimated_steps = 10000.00000000001 step_error = 2.499999999999999e-14 est_needed_step_err = 2.499999999999999e-14 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 8.329430301256176e-190 estimated_step_error = 8.329430301256176e-190 Double H and LOOP opt_iter = 2 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 NO REAL POLE (three term test) for Equation 1 NO COMPLEX POLE (six term test) for Equation 1 min_size = 0.6766027626496453 min_size = 1 glob_desired_digits_correct = 10 estimated_h = 2e-06 estimated_answer = 1 desired_abs_gbl_error = 1e-10 range = 0.01000000000000001 estimated_steps = 5000.000000000005 step_error = 3.535533905932736e-14 est_needed_step_err = 3.535533905932736e-14 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 max_estimated_step_error = 5.723923339350064e-179 estimated_step_error = 5.723923339350064e-179 Double H and LOOP opt_iter = 3 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 = 0.1584 Order of pole (three term test) = 5.783e-05 NO COMPLEX POLE (six term test) for Equation 1 SETTING H FOR POLE ACCURACY START of Soultion TOP MAIN SOLVE Loop x[1] = 1.8 y[1] (analytic) = 0.6766027626496453 y[1] (numeric) = 0.6766027626496453 absolute error = 0 relative error = 0 % Desired digits = 10 Estimated correct digits = 12 Correct digits = 16 h = 2e-06 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 NO REAL POLE (three term test) for Equation 1 NO COMPLEX POLE (six term test) for Equation 1 Finished! diff ( y , x , 1 ) = expt ( sin ( 2.0 * x + 3.0 ) , sin ( x + 1.0 ) ) * ( ( cos ( x + 1.0 ) * ln ( sin ( 2.0 * x + 3.0 ) ) ) + ( ( 2.0 * sin ( x + 1.0 ) * cos ( 2.0 * x + 3.0 ) ) / sin ( 2.0 * x + 3.0 ) ) ) ; Iterations = 5000 Total Elapsed Time = 0.0 Seconds Elapsed Time(since restart) = 0.0 Seconds Time to Timeout = 10 Minutes 0.0 Seconds Percent Done = 100 %