##############ECHO OF PROBLEM################# ##############temp/expt_sin_lin_new_sm_hpostode.ode################# diff ( y , x , 1 ) = ( expt ( sin ( x ) , ( 2.0 * x + 1.0 ) ) ) * ( ( 2.0 * ln ( sin ( x ) ) ) + ( ( ( 2.0 * x + 1.0 ) * cos ( x ) ) / sin ( x ) ) ) ; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=30 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK x_start=0.1 x_end=0.11 $array_y_init[0 + 1] = exact_soln_y(x_start) $glob_look_poles=true $glob_max_iter=1000000 $glob_max_h=0.00001 #END SECOND INPUT BLOCK #BEGIN OVERRIDE BLOCK $glob_desired_digits_correct=10 $glob_display_interval=0.01 $glob_look_poles=true $glob_max_iter=1000000000 $glob_max_minutes=30.0 $glob_subiter_method=3 #END OVERRIDE BLOCK ! #BEGIN USER DEF BLOCK def exact_soln_y (x) return(expt(sin(x),2.0*x+1.0)) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Optimize min_size 0.0 min_size 1.0 $glob_desired_digits_correct 10 estimated_h 1.0e-06 estimated_answer 1.0 desired_abs_gbl_error 1.0e-10 range 0.009999999999999995 estimated_steps 9999.999999999996 step_error 1.0000000000000005e-14 est_needed_step_err 1.0000000000000005e-14 opt_iter 1 hn_div_ho 0.5 hn_div_ho_2 0.25 hn_div_ho_3 0.125 max_estimated_step_error 3.7769694633364675e-133 estimated_step_error 3.7769694633364675e-133 best_h 2.0e-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.5397974895984017e-125 estimated_step_error 1.5397974895984017e-125 best_h 4.0e-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.0333306545976978e-117 estimated_step_error 1.0333306545976978e-117 best_h 8.0e-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 6.934431113027482e-110 estimated_step_error 6.934431113027482e-110 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 4.653438742818478e-102 estimated_step_error 4.653438742818478e-102 best_h 1.0e-05 START of Soultion TOP MAIN SOLVE Loop x[1] 0.1 y[1] (analytic) 0.06296962708006738 y[1] (numeric) 0.06296962708006738 absolute error 0.0 relative error 0.0% Correct digits 16 h 1.0e-05 NO INFO (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 0.10000068455822053 Order of pole (three term test) 0.004681940735909185 NO COMPLEX POLE (six term test) for Equation 1 Finished! diff ( y , x , 1 ) = ( expt ( sin ( x ) , ( 2.0 * x + 1.0 ) ) ) * ( ( 2.0 * ln ( sin ( x ) ) ) + ( ( ( 2.0 * x + 1.0 ) * cos ( x ) ) / sin ( x ) ) ) ; Iterations 1001 Total Elapsed Time 3 Seconds Elapsed Time(since restart) 3 Seconds Time to Timeout 29 Minutes 57.0 Seconds Percent Done 100.19999999996122%