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/sqrt_sqrtpostode.ode################# diff ( y , x , 1 ) = sqrt ( sqrt ( 0.1 * x + 0.2 ) ) ; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=40 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK x_start=c(0.1) x_end=c(0.5) $array_y_init[0 + 1] = exact_soln_y(x_start) $glob_look_poles=true $glob_min_h=c(0.001) # 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(0.8) * (c(x) + c(2.0)) * sqrt(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.1e0+/-0.50e-31 y[1] (closed_form) 0.1137271367855683019111487491957184e1+/-0.835e-29 y[1] (numeric) 0.1137271367855683019111487491957184e1+/-0.84e-29 absolute error 0.0e0+/-0.1675e-28 relative error 0.0e0+/-0.1472823503117114825097315474273738e-26% Desired digits 8 Estimated correct digits 14 Correct digits 32 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.209999999999999999999999999999996e1+/-0.554e-28 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 TOP MAIN SOLVE Loop x[1] 0.2e0+/-0.5e-31 y[1] (closed_form) 0.120536433735319444236116605375305e1+/-0.896e-29 y[1] (numeric) 0.1205364337353194442361283558715041e1+/-0.899e-29 absolute error 0.117504961991e-21+/-0.179e-28 relative error 0.9748501623087993688813205605040399e-20+/-0.148e-26% Desired digits 8 Estimated correct digits 13 Correct digits 32 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.2199999999999999999999999999999949e1+/-0.554e-28 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 TOP MAIN SOLVE Loop x[1] 0.3e0+/-0.5e-31 y[1] (closed_form) 0.1274235739780027645416240166918823e1+/-0.983e-29 y[1] (numeric) 0.1274235739780027645416466088072426e1+/-0.967e-29 absolute error 0.225921153603e-21+/-0.195e-28 relative error 0.1772993383798832647164375566509566e-19+/-0.153e-26% Desired digits 8 Estimated correct digits 13 Correct digits 32 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.2299999999999999999999999999999943e1+/-0.554e-28 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 TOP MAIN SOLVE Loop x[1] 0.4e0+/-0.5e-31 y[1] (closed_form) 0.1343860036446943959883905164110318e1+/-0.106e-28 y[1] (numeric) 0.1343860036446943959884231470635185e1+/-0.102e-28 absolute error 0.326306524867e-21+/-0.208e-28 relative error 0.2428128793305944029138249556411736e-19+/-0.154e-26% Desired digits 8 Estimated correct digits 13 Correct digits 32 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.2399999999999999999999999999999937e1+/-0.554e-28 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 ) = sqrt ( sqrt ( 0.1 * x + 0.2 ) ) ; Iterations 8 Total Elapsed Time 22 Seconds Elapsed Time(since restart) 21 Seconds Time to Timeout 19 Minutes 38 Seconds Percent Done 0.1125e3+/-0.586e-28%