Initializing... Initialized Initializing(2)... PI = 0.3141592653589793238462643383279486e1+/-0.313e-31 E = 0.2718281828459045235360287471352649e1+/-0.256e-30 LOG_E_10 = 0.2302585092994045684017991454684277e1+/-0.5e-30 Initialized(2) ##############ECHO OF PROBLEM################# ##############temp/cosh_sqrtpostode.ode################# diff ( y , x , 1 ) = cosh ( sqrt ( 0.1 * x + 0.2 ) ) ; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=30 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK x_start=c(-1.9) x_end=c(-1.8) $array_y_init[0 + 1] = exact_soln_y(x_start) $glob_look_poles=true $glob_max_h=0.01 $glob_upper_ratio_limit=c(1.001) $glob_lower_ratio_limit=c(0.999) # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb $glob_type_given_pole=1 # ELIMINATED in preodein.rb $array_given_rad_poles[1][1]=c(-2.0) # ELIMINATED in preodein.rb $array_given_rad_poles[1][2]=c(0.0) # ELIMINATED in preodein.rb $array_given_ord_poles[1][1]=c(0.5) # ELIMINATED in preodein.rb $array_given_ord_poles[1][2]=c(0.0) #END SECOND INPUT BLOCK #BEGIN OVERRIDE BLOCK $glob_desired_digits_correct=8 $glob_max_minutes=(20.0) $glob_subiter_method=3 $glob_max_iter=1000 $glob_upper_ratio_limit=c(1.11) $glob_lower_ratio_limit=c(0.99) # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb #END OVERRIDE BLOCK ! #BEGIN USER DEF BLOCK def exact_soln_y (x) x = c(x) return(c(20.0) * sqrt(c(0.1) * c(x) + c(0.2)) * sinh( sqrt(c(0.1) * c(x) + c(0.2))) - c(20.0) * cosh( sqrt(c(0.1) * c(x) + c(0.2)))) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion TOP MAIN SOLVE Loop x[1] -0.19e1+/-0.50e-30 y[1] (closed_form) -0.198997498610763839281121100923226e2+/-0.493e-28 y[1] (numeric) -0.198997498610763839281121100923226e2+/-0.511e-28 absolute error 0.0e0+/-0.1004e-27 relative error 0.0e0+/-0.5045289548909401789680683057031565e-27% Desired digits 8 Estimated correct digits 14 Correct digits 32 h 0.5e-1+/-0.1e-100 SETTING H FOR POLE Radius of convergence (given) for eq 1 0.1000000000000000000000000000001722e0+/-0.223e-29 Order of pole (given) 0.5e0+/-0.50e-31 Radius of convergence (ratio test) for eq 1 0.1114338328878997723530664462819e0+/-0.504e-27 Order of pole (ratio test) 0.9999999999999728521112889094602485e0+/-0.5e-31 Radius of convergence (three term test) for eq 1 0.996766451942924141692344241566037e-1+/-0.2e-25 Order of pole (three term test) 0.1637424832938150180327818351643668e1+/-0.574e-23 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] -0.18003233548057075858307655758433e1+/-0.20e-25 y[1] (closed_form) -0.197993254794571046461248672106749e2+/-0.364e-24 y[1] (numeric) -0.197993506688569563232438813510075e2+/-0.206e-25 absolute error 0.251893998516771190141403326e-4+/-0.384e-24 relative error 0.1272235252550018164648399610373748e-3+/-0.193e-23% Desired digits 8 Estimated correct digits 13 Correct digits 6 h 0.996766451942924141692344241566037e-2+/-0.2e-26 Radius of convergence (given) for eq 1 0.1996766451942924141692344241567207e0+/-0.398e-25 Order of pole (given) 0.5e0+/-0.50e-31 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 ) = cosh ( sqrt ( 0.1 * x + 0.2 ) ) ; Iterations 11 Total Elapsed Time 40 Seconds Elapsed Time(since restart) 37 Seconds Time to Timeout 19 Minutes 20 Seconds Percent Done 0.11961197423315089700308130898804e3+/-0.24e-22%