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=(5.0) $glob_subiter_method=3 $glob_max_iter=1000 $glob_upper_ratio_limit=c(1.1) $glob_lower_ratio_limit=c(0.9) # 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 Optimize min_size 0.0e0+/-0.50e-31 min_size 0.1e1+/-0.50e-31 $glob_desired_digits_correct 8 estimated_h 0.1e-2+/-0.5e-31 estimated_answer 0.1e1+/-0.50e-31 desired_abs_gbl_error 0.9999999999999999999999999999751083e-8+/-0.815e-35 range 0.3e0+/-0.100e-29 estimated_steps 0.3e3+/-0.16e-25 step_error 0.2886751345948128822545743902436858e-10+/-0.152e-34 est_needed_step_err 0.2886751345948128822545743902436858e-10+/-0.152e-34 opt_iter 1 SETTING H FOR MIN H SETTING H FOR DISPLAY INTERVAL 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 18 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 Radius of convergence (ratio test) for eq 1 0.1071428571428571246428571428571437e1+/-0.206e-28 Order of pole (ratio test) 0.1000000000000000000085052455426806e1+/-0.5e-31 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.95e0+/-0.5e-30 y[1] (closed_form) -0.116635077999708412305000728969236e1+/-0.128e-29 y[1] (numeric) -0.11663505981928060096255835684548e1+/-0.222e-29 absolute error 0.18180427811342442372123756e-6+/-0.350e-29 relative error 0.155874442947497265287996053005002e-4+/-0.3e-27% Desired digits 8 Estimated correct digits 17 Correct digits 7 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.950000000000000000000000000000013e0+/-0.498e-29 Order of pole (given) 0.4e1+/-0.50e-31 Radius of convergence (ratio test) for eq 1 0.1017857142857142684107142857142869e1+/-0.154e-28 Order of pole (ratio test) 0.1000000000000000000054931395215897e1+/-0.5e-31 Radius of convergence (three term test) for eq 1 0.9500000000000040755000000000172834e0+/-0.425e-27 Order of pole (three term test) 0.530400000000002248792523716909783e-13+/-0.761e-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 17 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 Radius of convergence (ratio test) for eq 1 0.9642857142857141217857142857143003e0+/-0.176e-28 Order of pole (ratio test) 0.999999999999999999701301779979624e0+/-0.5e-31 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.85e0+/-0.5e-30 y[1] (closed_form) -0.162833299409729289639731325055973e1+/-0.294e-29 y[1] (numeric) -0.162833204740408209433271507193887e1+/-0.381e-29 absolute error 0.94669321080206459817862086e-6+/-0.675e-29 relative error 0.5813879680838179213564455356475033e-4+/-0.414e-27% Desired digits 8 Estimated correct digits 16 Correct digits 7 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.8500000000000000000000000000000193e0+/-0.497e-29 Order of pole (given) 0.4e1+/-0.50e-31 Radius of convergence (ratio test) for eq 1 0.9107142857142855594642857142857217e0+/-0.157e-28 Order of pole (ratio test) 0.9999999999999999978508110335037892e0+/-0.5e-31 Radius of convergence (three term test) for eq 1 0.8500000000000036465000000000153669e0+/-0.443e-27 Order of pole (three term test) 0.530400000000002235152299704839314e-13+/-0.879e-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 16 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 Radius of convergence (ratio test) for eq 1 0.8571428571428569971428571428571554e0+/-0.186e-28 Order of pole (ratio test) 0.9999999999999999894510723468686163e0+/-0.5e-31 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 TOP MAIN SOLVE Loop x[1] -0.75e0+/-0.5e-30 y[1] (closed_form) -0.237037037037037037037037037037036e1+/-0.326e-29 y[1] (numeric) -0.237036723293478861798742924339543e1+/-0.654e-29 absolute error 0.313743558175238294112697493e-5+/-0.980e-29 relative error 0.1323605636051786553287942548593755e-3+/-0.413e-27% Desired digits 8 Estimated correct digits 15 Correct digits 6 h 0.5e-1+/-0.1e-100 Radius of convergence (given) for eq 1 0.7500000000000000000000000000000226e0+/-0.498e-29 Order of pole (given) 0.4e1+/-0.50e-31 Radius of convergence (ratio test) for eq 1 0.8035714285714284348214285714285818e0+/-0.136e-28 Order of pole (ratio test) 0.9999999999999999521802576238774613e0+/-0.5e-31 Radius of convergence (three term test) for eq 1 0.7500000000000032175000000000137081e0+/-0.355e-27 Order of pole (three term test) 0.530400000000002260956605939312705e-13+/-0.806e-26 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 21 Seconds Elapsed Time(since restart) 18 Seconds Time to Timeout 4 Minutes 39 Seconds Percent Done 0.1166666666666666666666666666666666e3+/-0.725e-27%