Initializing... Initialized Initializing(2)... PI = 0.314159265358979323846264338327950271e1+/-0.313e-33 E = 0.271828182845904523536028747135266232e1+/-0.256e-32 LOG_E_10 = 0.230258509299404568401799145468436324e1+/-0.5e-32 Initialized(2) ##############ECHO OF PROBLEM################# ##############temp/mtest7postode.ode################# diff ( y2 , x , 5 ) = y1 ; diff ( y1 , x , 1 ) = neg ( diff ( y2 , x , 4 ) ) ; ! #BEGIN FIRST INPUT BLOCK # Digits:=64; ELIMINATED in preodein.rb max_terms=40 ! #END FIRST INPUT BLOCK #BEGIN SECOND INPUT BLOCK x_start=c(0.1) x_end=c(0.7) # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb # ELIMINATED in preodein.rb $array_y1_init[0 + 1] = exact_soln_y1(x_start) $array_y2_init[0 + 1] = exact_soln_y2(x_start) $array_y2_init[1 + 1] = exact_soln_y2p(x_start) $array_y2_init[2 + 1] = exact_soln_y2pp(x_start) $array_y2_init[3 + 1] = exact_soln_y2ppp(x_start) $array_y2_init[4 + 1] = exact_soln_y2pppp(x_start) $glob_look_poles=true $glob_max_h=c(0.1) # 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=3 #END SECOND INPUT BLOCK #BEGIN OVERRIDE BLOCK $glob_desired_digits_correct=8 $glob_max_minutes=(3.0) $glob_subiter_method=3 $glob_max_iter=100000 $glob_upper_ratio_limit=c(1.000001) $glob_lower_ratio_limit=c(0.999999) $glob_look_poles=true $glob_h=c(0.001) $glob_display_interval=c(0.01) #END OVERRIDE BLOCK ! #BEGIN USER DEF BLOCK def exact_soln_y1 (x) x = c(x) return( cos(c(x)) ) end def exact_soln_y2 (x) x = c(x) return( sin(c(x)) + c(10.0) * c(x) + c(10.0)) end def exact_soln_y2p (x) x = c(x) return( cos(c(x)) + c(10.0)) end def exact_soln_y2pp (x) x = c(x) return( neg(sin(c(x)))) end def exact_soln_y2ppp (x) x = c(x) return( neg(cos(c(x)))) end def exact_soln_y2pppp (x) x = c(x) return( sin(c(x))) end #END USER DEF BLOCK #######END OF ECHO OF PROBLEM################# START of Soultion TOP MAIN SOLVE Loop x[1] 0.1e0+/-0.50e-33 y2[1] (closed_form) 0.11099833416646828152306814198410622e2+/-0.11e-31 y2[1] (numeric) 0.11099833416646828152306814198410622e2+/-0.22e-31 $glob_prec 0.1e-15+/-0.50e-49% absolute error 0.0e0+/-0.33e-31 relative error 0.0e0+/-0.297301759056209682175297451893142239e-30% Desired digits 8 Estimated correct digits 13 Correct digits 34 h 0.1e-2+/-0.5e-33 y1[1] (closed_form) 0.995004165278025766095561987803870286e0+/-0.106e-33 y1[1] (numeric) 0.995004165278025766095561987803870286e0+/-0.11e-32 $glob_prec 0.1e-15+/-0.50e-49% absolute error 0.0e0+/-0.1206e-32 relative error 0.0e0+/-0.12120552275909492467289276061818329380819e-30% Desired digits 8 Estimated correct digits 13 Correct digits 34 h 0.1e-2+/-0.5e-33 NO POLE (given) for Equation 1 NO POLE (ratio test) for Equation 1 Radius of convergence (three term test) for eq 1 0.993346653975306341535444538841937653e-1+/-0.333e-29 Order of pole (three term test) 0.329900332889206208113713132628938511e2+/-0.132e-26 NO COMPLEX POLE (six term test) for Equation 1 NO POLE (given) for Equation 2 NO POLE (ratio test) for Equation 2 Radius of convergence (three term test) for eq 2 0.99334665397530645717423702285909619e-1+/-0.367e-29 Order of pole (three term test) 0.320099667110793791922104186212057034e2+/-0.141e-26 NO COMPLEX POLE (six term test) for Equation 2 Finished! Maximum Time Reached before Solution Completed! diff ( y2 , x , 5 ) = y1 ; diff ( y1 , x , 1 ) = neg ( diff ( y2 , x , 4 ) ) ; Iterations 2 Total Elapsed Time 1 Hours 2 Minutes 16 Seconds Elapsed Time(since restart) 1 Hours 2 Minutes 6 Seconds Expected Time Remaining 8 Days 14 Hours 31 Minutes 4 Seconds Optimized Time Remaining 8 Days 13 Hours 57 Minutes 54 Seconds Expected Total Time 8 Days 15 Hours 0 Minutes 10 Seconds Time to Timeout 0.0 Seconds Percent Done 0.5e0+/-0.417e-30%