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/sing1postode.ode################# diff ( y , x , 1 ) = m1 * 2.0 * x / ( x * x + 0.000001 ) / ( x * x + 0.000001 ) ; ! #BEGIN FIRST INPUT BLOCK # Digits:=32; ELIMINATED in preodein.rb max_terms=30 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK x_start=c(-2.0) x_end=c(-1.5) $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 $glob_type_given_pole=2 # 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.001) # ELIMINATED in preodein.rb $array_given_ord_poles[1][1]=c(1.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(0.000001))) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion TOP MAIN SOLVE Loop x[1] -0.2e1+/-0.50e-30 y[1] (closed_form) 0.2499999375000156249960937509765622e0+/-0.137e-30 y[1] (numeric) 0.2499999375000156249960937509765622e0+/-0.161e-30 absolute error 0.0e0+/-0.298e-30 relative error 0.0e0+/-0.1192000298000000000000000000001306e-27% 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.2000000249999984375001953124694794e1+/-0.504e-28 Order of pole (given) 0.1e1+/-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.19e1+/-0.5e-30 y[1] (closed_form) 0.2770082335157247878878703911716367e0+/-0.159e-30 y[1] (numeric) 0.2770082335157247878845344959007577e0+/-0.302e-30 absolute error 0.3335895270879e-20+/-0.461e-30 relative error 0.12042585263768460879e-17+/-0.166e-27% 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.1900000263157876512613691939990842e1+/-0.503e-28 Order of pole (given) 0.1e1+/-0.50e-31 NO POLE (ratio test) for Equation 1 NO REAL POLE (three term test) for Equation 1 Radius of convergence (six term test) for eq 1 0.1900000659104016953755395913259557e1+/-0.366e-8 Order of pole (six term test) 0.99999999985346711241214372068712e0+/-0.755e-17 TOP MAIN SOLVE Loop x[1] -0.18e1+/-0.5e-30 y[1] (closed_form) 0.3086418800488024540733166440380728e0+/-0.186e-30 y[1] (numeric) 0.3086418800488024540656540509545395e0+/-0.463e-30 absolute error 0.76625930835333e-20+/-0.649e-30 relative error 0.248268092532409755333e-17+/-0.21e-27% 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.1800000277777756344310577866896444e1+/-0.504e-28 Order of pole (given) 0.1e1+/-0.50e-31 NO POLE (ratio test) for Equation 1 NO REAL POLE (three term test) for Equation 1 Radius of convergence (six term test) for eq 1 0.1800000614436112272496228994212047e1+/-0.359e-8 Order of pole (six term test) 0.9999999998534675563933396916136e0+/-0.881e-17 TOP MAIN SOLVE Loop x[1] -0.17e1+/-0.5e-30 y[1] (closed_form) 0.3460206415153489566266586066925236e0+/-0.22e-30 y[1] (numeric) 0.346020641515348956613303472411508e0+/-0.64e-30 absolute error 0.133551342810156e-19+/-0.86e-30 relative error 0.385963514272693650156e-17+/-0.248e-27% 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.1700000294117621616124898492347825e1+/-0.503e-28 Order of pole (given) 0.1e1+/-0.50e-31 NO POLE (ratio test) for Equation 1 NO REAL POLE (three term test) for Equation 1 Radius of convergence (six term test) for eq 1 0.1700000577723087296924815849720853e1+/-0.329e-8 Order of pole (six term test) 0.9999999998534680805548213625531e0+/-0.956e-17 TOP MAIN SOLVE Loop x[1] -0.16e1+/-0.5e-30 y[1] (closed_form) 0.3906248474121689796214923353545565e0+/-0.263e-30 y[1] (numeric) 0.3906248474121689796005272139862725e0+/-0.862e-30 absolute error 0.20965121368284e-19+/-0.112e-29 relative error 0.53670731667928408284e-17+/-0.286e-27% 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.1600000312499969482427835463023206e1+/-0.502e-28 Order of pole (given) 0.1e1+/-0.50e-31 NO POLE (ratio test) for Equation 1 NO REAL POLE (three term test) for Equation 1 Radius of convergence (six term test) for eq 1 0.160000054894105614083684010767411e1+/-0.29e-8 Order of pole (six term test) 0.99999999985346870558264329420594e0+/-0.1e-16 Finished! diff ( y , x , 1 ) = m1 * 2.0 * x / ( x * x + 0.000001 ) / ( x * x + 0.000001 ) ; Iterations 10 Total Elapsed Time 45 Seconds Elapsed Time(since restart) 45 Seconds Time to Timeout 19 Minutes 15 Seconds Percent Done 0.11e3+/-0.424e-27%