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/sing5postode.ode################# diff ( y , x , 1 ) = m1 * 3.0 / x / x / x / x ; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=20 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK x_start=c(-1.0) x_end=c(-0.7) $array_y_init[0 + 1] = exact_soln_y(x_start) $glob_look_poles=true $glob_max_h=c(0.001) $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 # 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.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(4.0) # 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(1.0)/c(x)/c(x)/c(x)) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion TOP MAIN SOLVE Loop x[1] -0.1e1+/-0.50e-30 y[1] (closed_form) -0.1e1+/-0.155e-29 y[1] (numeric) -0.1e1+/-0.165e-29 absolute error 0.0e0+/-0.320e-29 relative error 0.0e0+/-0.320000000000000000000000000000992e-27% 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.1e1+/-0.1e-31 Order of pole (given) 0.4e1+/-0.50e-31 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 0.1000000000000004290000000000018223e1+/-0.553e-27 Order of pole (three term test) 0.530400000000002253139380210077595e-13+/-0.939e-26 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] -0.9e0+/-0.5e-30 y[1] (closed_form) -0.137174211248285322359396433470507e1+/-0.235e-29 y[1] (numeric) -0.137174164060658904766435016479489e1+/-0.291e-29 absolute error 0.47187626417592961416991018e-6+/-0.526e-29 relative error 0.3439977965842526887298645212200013e-4+/-0.383e-27% Desired digits 8 Estimated correct digits 14 Correct digits 7 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.9000000000000000000000000000000154e0+/-0.499e-29 Order of pole (given) 0.4e1+/-0.50e-31 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 0.9000000000000038610000000000164529e0+/-0.478e-27 Order of pole (three term test) 0.530400000000002262107035729195166e-13+/-0.906e-26 NO COMPLEX POLE (six term test) for Equation 1 TOP MAIN SOLVE Loop x[1] -0.8e0+/-0.5e-30 y[1] (closed_form) -0.1953125e1+/-0.375e-29 y[1] (numeric) -0.195312325355389037499363750504986e1+/-0.498e-29 absolute error 0.174644610962500636249495014e-5+/-0.873e-29 relative error 0.89418040812800325759741447168e-4+/-0.446e-27% Desired digits 8 Estimated correct digits 13 Correct digits 7 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.8000000000000000000000000000000211e0+/-0.5e-29 Order of pole (given) 0.4e1+/-0.50e-31 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 0.8000000000000034320000000000146081e0+/-0.509e-27 Order of pole (three term test) 0.530400000000002258842240472724302e-13+/-0.107e-25 NO COMPLEX POLE (six term test) for Equation 1 Finished! diff ( y , x , 1 ) = m1 * 3.0 / x / x / x / x ; Iterations 6 Total Elapsed Time 7 Seconds Elapsed Time(since restart) 7 Seconds Time to Timeout 19 Minutes 53 Seconds Percent Done 0.1166666666666666666666666666666666e3+/-0.725e-27%