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/sin_sqrt_linpostode.ode################# diff ( y , x , 1 ) = sin ( sqrt ( 2.0 * x + 3.0 ) ) ; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=40 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK x_start=c(2.0) x_end=c(3.0) $array_y_init[0 + 1] = exact_soln_y(x_start) $glob_look_poles=true $glob_max_h=c(0.0001) $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(-1.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.5) # ELIMINATED in preodein.rb $array_given_ord_poles[1][2]=c(0.0) # 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(sin(sqrt(c(2.0)*c(x)+c(3.0)))-sqrt(c(2.0)*c(x)+c(3.0))*cos(sqrt(c(2.0)*c(x)+c(3.0)))) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion TOP MAIN SOLVE Loop x[1] 0.2e1+/-0.50e-31 y[1] (closed_form) 0.2802891969591020065097816066030565e1+/-0.11e-27 y[1] (numeric) 0.2802891969591020065097816066030565e1+/-0.11e-27 absolute error 0.0e0+/-0.22e-27 relative error 0.0e0+/-0.7849036009479201667626776731971292e-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.3500000000000000000000000000004443e1+/-0.109e-26 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.21e1+/-0.5e-31 y[1] (closed_form) 0.2848803743596289649235895872757115e1+/-0.119e-27 y[1] (numeric) 0.2849638774800133183177766968286002e1+/-0.11e-27 absolute error 0.835031203843533941871095528887e-3+/-0.229e-27 relative error 0.2931164372837429413254941020362485e-1+/-0.803e-26% Desired digits 8 Estimated correct digits 13 Correct digits 4 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.3600000000000000000000000000004026e1+/-0.11e-26 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.22e1+/-0.5e-31 y[1] (closed_form) 0.289137377664734129030921453408717e1+/-0.251e-27 y[1] (numeric) 0.2893047669197997212802394974705833e1+/-0.11e-27 absolute error 0.1673892550655922493180440618663e-2+/-0.361e-27 relative error 0.5789263789328770421994692395044082e-1+/-0.124e-25% Desired digits 8 Estimated correct digits 13 Correct digits 4 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.3700000000000000000000000000003928e1+/-0.11e-26 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.23e1+/-0.5e-31 y[1] (closed_form) 0.2930589588687642502509775385854991e1+/-0.263e-27 y[1] (numeric) 0.2933104889779928495350854515402621e1+/-0.11e-27 absolute error 0.251530109228599284107912954763e-2+/-0.373e-27 relative error 0.8582918270082227908575740172396917e-1+/-0.127e-25% Desired digits 8 Estimated correct digits 13 Correct digits 4 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.3800000000000000000000000000003889e1+/-0.11e-26 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.24e1+/-0.5e-31 y[1] (closed_form) 0.2966443639457176103385369729339699e1+/-0.278e-27 y[1] (numeric) 0.2969801694171113540129420103320507e1+/-0.11e-27 absolute error 0.3358054713937436744050373980808e-2+/-0.388e-27 relative error 0.1132013657455471093531826583763222e0+/-0.13e-25% Desired digits 8 Estimated correct digits 13 Correct digits 3 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.3900000000000000000000000000003747e1+/-0.109e-26 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.25e1+/-0.5e-31 y[1] (closed_form) 0.2998933019437576951297737468831758e1+/-0.294e-27 y[1] (numeric) 0.3003134045866608071939642512270419e1+/-0.11e-27 absolute error 0.4201026429031120641905043438661e-2+/-0.404e-27 relative error 0.1400840366157622803921311789009955e0+/-0.134e-25% Desired digits 8 Estimated correct digits 13 Correct digits 3 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4000000000000000000000000000003566e1+/-0.11e-26 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.26e1+/-0.5e-31 y[1] (closed_form) 0.3028059162193812079676282163953467e1+/-0.308e-27 y[1] (numeric) 0.3033102321415447792575565977731711e1+/-0.11e-27 absolute error 0.5043159221635712899283813778244e-2+/-0.418e-27 relative error 0.1665475788782796435008796420368955e0+/-0.138e-25% Desired digits 8 Estimated correct digits 13 Correct digits 3 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4100000000000000000000000000003566e1+/-0.109e-26 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.27e1+/-0.5e-31 y[1] (closed_form) 0.3053827576101420111789666290720214e1+/-0.323e-27 y[1] (numeric) 0.3059711037470738108256923857152375e1+/-0.11e-27 absolute error 0.5883461369317996467257566432161e-2+/-0.433e-27 relative error 0.1926585972096350570971185590020772e0+/-0.141e-25% Desired digits 8 Estimated correct digits 13 Correct digits 3 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4200000000000000000000000000003443e1+/-0.109e-26 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.28e1+/-0.5e-31 y[1] (closed_form) 0.3076247593690195119964165585010672e1+/-0.342e-27 y[1] (numeric) 0.3082968595879279498053736467941938e1+/-0.11e-27 absolute error 0.6721002189084378089570882931266e-2+/-0.452e-27 relative error 0.2184805346250436253983673658202787e0+/-0.147e-25% Desired digits 8 Estimated correct digits 13 Correct digits 3 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4300000000000000000000000000003059e1+/-0.109e-26 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.29e1+/-0.5e-31 y[1] (closed_form) 0.3095332137042335591296190900056747e1+/-0.354e-27 y[1] (numeric) 0.3102887045199594413988861117074052e1+/-0.11e-27 absolute error 0.7554908157258822692670217017305e-2+/-0.464e-27 relative error 0.2440742325144376718992980236191267e0+/-0.15e-25% Desired digits 8 Estimated correct digits 13 Correct digits 3 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4400000000000000000000000000003145e1+/-0.109e-26 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 ) = sin ( sqrt ( 2.0 * x + 3.0 ) ) ; Iterations 20 Total Elapsed Time 1 Minutes 14 Seconds Elapsed Time(since restart) 1 Minutes 13 Seconds Time to Timeout 18 Minutes 46 Seconds Percent Done 0.105e3+/-0.257e-28%