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/sqrtpostode.ode#################
diff  (  y  ,  x  ,  1  )   =  sqrt ( x )  ; 
!
#BEGIN FIRST INPUT BLOCK
# Digits:=32; ELIMINATED in preodein.rb
max_terms=40 
!
#END FIRST INPUT BLOCK
#BEGIN SECOND INPUT BLOCK
x_start=c(0.1) 
x_end=c(0.2) 
$array_y_init[0 + 1] = exact_soln_y(x_start)
$glob_look_poles=true 
$glob_min_h=c(0.001) 
$glob_max_h=c(0.01) 
#  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(1.5) 
#  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(2.0)*expt(c(x),c(1.5))/c(3.0)) 
end
#END USER DEF BLOCK
#######END OF ECHO OF PROBLEM#################
START of Soultion
 
 
TOP MAIN SOLVE Loop
x[1]                             0.1e0+/-0.50e-31 
y[1] (closed_form)               0.2108185106778919554665929029632306e-1+/-0.326e-30 
y[1] (numeric)                   0.2108185106778919554665929029632306e-1+/-0.327e-30 
absolute error                   0.0e0+/-0.653e-30 
relative error                   0.0e0+/-0.3097450968134927555692916226852225e-26%
Desired digits                   8 
Estimated correct digits         13 
Correct digits                   32 
h                                0.5e-1+/-0.1e-100 
SETTING H FOR POLE
Radius of convergence (given) for eq  1           0.1000000000000000000000000000001722e0+/-0.184e-30 
Order of pole (given)                             0.15e1+/-0.50e-31 
Radius of convergence (ratio test) for eq  1      0.107692307692307697907692307692665e0+/-0.614e-28 
Order of pole (ratio test)                        0.9999999999999991529729635177658587e0+/-0.5e-31 
Radius of convergence (three term test) for eq 1  0.100000000000000192100000000000796e0+/-0.534e-28 
Order of pole (three term test)                   0.1499999999999940433499999999856542e1+/-0.142e-23 
NO COMPLEX POLE (six term test) for Equation 1
Finished!
diff  (  y  ,  x  ,  1  )   =  sqrt ( x )  ; 
Iterations                      10 
Total Elapsed Time              19 Seconds
Elapsed Time(since restart)     18 Seconds
Time to Timeout                 19 Minutes 41 Seconds
Percent Done                    0.1100000000000002113100000000008756e3+/-0.584e-25%