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/div_c_linpostode.ode################# diff ( y , x , 1 ) = 2.0 / ( 0.2 * x + 0.3 ) ; ! #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.5) x_end=c(3.1) $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 # 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(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(10.0) * ln(c(0.2) * c(x) + c(0.3))) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion TOP MAIN SOLVE Loop x[1] 0.25e1+/-0.50e-31 y[1] (closed_form) -0.223143551314209755766295090309808e1+/-0.232e-28 y[1] (numeric) -0.223143551314209755766295090309808e1+/-0.234e-28 absolute error 0.0e0+/-0.466e-28 relative error 0.0e0+/-0.2088341774859640200692160017648203e-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.4000000000000000000000000000003566e1+/-0.11e-26 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.26e1+/-0.5e-31 y[1] (closed_form) -0.198450938723838254751987414873122e1+/-0.231e-28 y[1] (numeric) -0.198450938722625270359539845663631e1+/-0.24e-28 absolute error 0.1212984392447569209491e-10+/-0.471e-28 relative error 0.6112263314287177617133616220506285e-9+/-0.237e-26% Desired digits 8 Estimated correct digits 13 Correct digits 12 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4100000000000000000000000000003566e1+/-0.109e-26 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.27e1+/-0.5e-31 y[1] (closed_form) -0.17435338714477775270092068608665e1+/-0.232e-28 y[1] (numeric) -0.174353387142517637511768905078467e1+/-0.24e-28 absolute error 0.2260115189151781008183e-10+/-0.472e-28 relative error 0.1296284073492102598477873411680456e-8+/-0.27e-26% Desired digits 8 Estimated correct digits 13 Correct digits 11 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4200000000000000000000000000003443e1+/-0.109e-26 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.28e1+/-0.5e-31 y[1] (closed_form) -0.150822889734583635145908274564685e1+/-0.232e-28 y[1] (numeric) -0.15082288973141637967747114322438e1+/-0.24e-28 absolute error 0.3167255468437131340305e-10+/-0.472e-28 relative error 0.2099983281059546483585437456262074e-8+/-0.312e-26% Desired digits 8 Estimated correct digits 13 Correct digits 11 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4300000000000000000000000000003059e1+/-0.109e-26 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.29e1+/-0.5e-31 y[1] (closed_form) -0.127833371509884895722342967029053e1+/-0.232e-28 y[1] (numeric) -0.127833371505929135445248939999573e1+/-0.24e-28 absolute error 0.395576027709402702948e-10+/-0.472e-28 relative error 0.3094466046206207029163122009820466e-8+/-0.369e-26% Desired digits 8 Estimated correct digits 13 Correct digits 11 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4400000000000000000000000000003145e1+/-0.109e-26 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.3e1+/-0.5e-31 y[1] (closed_form) -0.105360515657826301227500980839296e1+/-0.233e-28 y[1] (numeric) -0.105360515653182959596189362647028e1+/-0.24e-28 absolute error 0.4643341631311618192268e-10+/-0.473e-28 relative error 0.4407098429919942588806036352529719e-8+/-0.448e-26% Desired digits 8 Estimated correct digits 13 Correct digits 11 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.4500000000000000000000000000002758e1+/-0.109e-26 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 Finished! diff ( y , x , 1 ) = 2.0 / ( 0.2 * x + 0.3 ) ; Iterations 12 Total Elapsed Time 13 Seconds Elapsed Time(since restart) 13 Seconds Time to Timeout 19 Minutes 47 Seconds Percent Done 0.1083333333333333333333333333333333e3+/-0.4e-28%