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/nonlinear2postode.ode################# diff ( y , x , 1 ) = y * y ; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=20 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb x_start=c(0.0) x_end=c(0.2) $array_y_init[0 + 1] = exact_soln_y(x_start) $glob_look_poles=true # 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(0.5) # 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.0) # ELIMINATED in preodein.rb $array_given_ord_poles[1][2]=c(0.0) # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb #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(2.0)/(c(1.0) - c(2.0)*c(x))) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion TOP MAIN SOLVE Loop x[1] 0.0e0+/-0.50e-31 y[1] (closed_form) 0.2e1+/-0.15e-30 y[1] (numeric) 0.2e1+/-0.35e-30 absolute error 0.0e0+/-0.50e-30 relative error 0.0e0+/-0.2500000000000000000000000000000375e-28% 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.500000000000000000000000000000312e0+/-0.5e-30 Order of pole (given) 0.0e0+/-0.50e-31 Radius of convergence (ratio test) for eq 1 0.4999999999999999300449272949272976e0+/-0.54e-29 Order of pole (ratio test) 0.9999999999999919999999999999999167e0+/-0.5e-31 Radius of convergence (three term test) for eq 1 0.5000000000000021953306415806511096e0+/-0.16e-27 Order of pole (three term test) 0.100000000000005850000000000025395e1+/-0.578e-26 NO COMPLEX POLE (six term test) for Equation 1 SETTING H FOR POLE SETTING H FOR POLE TOP MAIN SOLVE Loop x[1] 0.1355000000000003556568132125965441e0+/-0.198e-25 y[1] (closed_form) 0.2743484224965709124112089301328961e1+/-0.149e-24 y[1] (numeric) 0.2743484224965709123146909202062172e1+/-0.148e-24 absolute error 0.965180099266789e-18+/-0.297e-24 relative error 0.351808146182744247227121718556234e-16+/-0.108e-22% Desired digits 8 Estimated correct digits 13 Correct digits 32 h 0.4050000000000015807232358586253256e-1+/-0.198e-25 SETTING H FOR POLE Radius of convergence (given) for eq 1 0.364499999999999644343186787403917e0+/-0.393e-24 Order of pole (given) 0.0e0+/-0.50e-31 Radius of convergence (ratio test) for eq 1 0.3644999999999995934741728546891158e0+/-0.566e-23 Order of pole (ratio test) 0.999999999999922395132779839221742e0+/-0.5e-31 Radius of convergence (three term test) for eq 1 0.3645000000000012448674585689801679e0+/-0.168e-21 Order of pole (three term test) 0.100000000000005850000000000025383e1+/-0.828e-20 NO COMPLEX POLE (six term test) for Equation 1 SETTING H FOR POLE Finished! diff ( y , x , 1 ) = y * y ; Iterations 5 Total Elapsed Time 9 Seconds Elapsed Time(since restart) 9 Seconds Time to Timeout 19 Minutes 51 Seconds Percent Done 0.1187800000000003361093399726150681e3+/-0.142e-16%