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/mtest7postode.ode################# diff ( y2 , x , 5 ) = y1 ; diff ( y1 , x , 1 ) = m1 * 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=(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_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( $glob__m1 * sin(c(x))) end def exact_soln_y2ppp (x) x = c(x) return( $glob__m1 * 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-31 y2[1] (closed_form) 0.1109983341664682815230681419841062e2+/-0.11e-29 y2[1] (numeric) 0.1109983341664682815230681419841062e2+/-0.22e-29 absolute error 0.0e0+/-0.33e-29 relative error 0.0e0+/-0.2973017590562096821752974518932006e-28% Desired digits 8 Estimated correct digits 13 Correct digits 32 h 0.5e-1+/-0.1e-100 y1[1] (closed_form) 0.9950041652780257660955619878038697e0+/-0.106e-31 y1[1] (numeric) 0.9950041652780257660955619878038697e0+/-0.11e-30 absolute error 0.0e0+/-0.1206e-30 relative error 0.0e0+/-0.1212055227590949246728927606181859e-28% Desired digits 8 Estimated correct digits 13 Correct digits 32 h 0.5e-1+/-0.1e-100 NO POLE (given) for Equation 1 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 NO POLE (given) for Equation 2 NO POLE (ratio test) for Equation 2 NO REAL POLE (three term test) for Equation 2 NO COMPLEX POLE (six term test) for Equation 2 TOP MAIN SOLVE Loop x[1] 0.2e0+/-0.5e-31 y2[1] (closed_form) 0.1219866933079506121545941262711838e2+/-0.116e-29 y2[1] (numeric) 0.1219866920234743486023291042186967e2+/-0.708e-29 absolute error 0.12844762635522650220524871e-6+/-0.824e-29 relative error 0.1052964244476776811563582461587852e-5+/-0.675e-28% Desired digits 8 Estimated correct digits 13 Correct digits 8 h 0.5e-1+/-0.1e-100 y1[1] (closed_form) 0.9800665778412416311241965167481681e0+/-0.182e-31 y1[1] (numeric) -0.860096699880529201028603416069507e2+/-0.138e-23 absolute error 0.869897365658941617339845381236988e2+/-0.138e-23 relative error 0.8875900732938308295279188941910969e4+/-0.14e-21% Desired digits 8 Estimated correct digits 13 Correct digits 0 h 0.5e-1+/-0.1e-100 NO POLE (given) for Equation 1 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 NO POLE (given) for Equation 2 NO POLE (ratio test) for Equation 2 NO REAL POLE (three term test) for Equation 2 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 ) = m1 * diff ( y2 , x , 4 ) ; Iterations 3 Total Elapsed Time 33 Minutes 1 Seconds Elapsed Time(since restart) 32 Minutes 57 Seconds Expected Time Remaining 1 Hours 6 Minutes 2 Seconds Optimized Time Remaining 1 Hours 5 Minutes 54 Seconds Expected Total Time 1 Hours 38 Minutes 55 Seconds Time to Timeout 0.0 Seconds Percent Done 0.3333333333333333333333333333333333e2+/-0.238e-28%