##############ECHO OF PROBLEM#################
##############temp/mtest6_rev_sm_hpostode.ode#################
diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1)+x1;
diff (x1,t,1) = 4.0 * x2 -  2.0 * diff (x2,t ,1) - 2.0 * x1;
!
#BEGIN FIRST INPUT BLOCK
# Digits:=64; ELIMINATED in preodein.rb
max_terms=30 
!
#END FIRST INPUT BLOCK
#BEGIN SECOND INPUT BLOCK
## problem from Boyce DePrima -
## _Elementary Differential Equations and Boundary Value Problems_
## page 269
##
t_start=1.5 
## did poorly with t_start := 0.5;
t_end=5.0 
$array_x1_init[0 + 1] = exact_soln_x1(t_start)
## I think following line should be omitted
## diff(x1,1,exact_soln_x1p(t_start));
$array_x2_init[0 + 1] = exact_soln_x2(t_start)
$array_x2_init[1 + 1] = exact_soln_x2p(t_start)
$glob_max_h=0.0000001 
$glob_look_poles=true 
$glob_max_iter=100 
#END SECOND INPUT BLOCK
#BEGIN OVERRIDE  BLOCK
$glob_desired_digits_correct=10 
$glob_display_interval=0.001 
$glob_look_poles=true 
$glob_max_iter=10000000 
$glob_max_minutes=3 
$glob_subiter_method=3 
#END OVERRIDE BLOCK
!
#BEGIN USER DEF BLOCK
def exact_soln_x1 (t)
 
c1 = 1.0;
c2 = 0.0002;
c3 = 0.0003;
return(2.0 * c1 + 6.0 * c3 * exp(-t)) 
end
def exact_soln_x1p (t)
 
c1 = 1.0;
c2 = 0.0002;
c3 = 0.0003;
return(   - 6.0 * c3 * exp(-t)) 
end
def exact_soln_x2 (t)
 
c1 = 1.0;
c2 = 0.0002;
c3 = 0.0003;
return(c1 + c2 * exp(2.0 * t) + c3 * exp(-t)) 
end
def exact_soln_x2p (t)
 
c1 = 1.0;
c2 = 0.0002;
c3 = 0.0003;
return(   2.0 * c2 * exp(2.0 * t) - c3 * exp(-t)) 
end
#END USER DEF BLOCK
#######END OF ECHO OF PROBLEM#################
START of Optimize
min_size                      0.0
min_size                      1.0
opt_iter                        1
$glob_desired_digits_correct  10
desired_abs_gbl_error         1.0e-10
range                         3.5
estimated_steps               35000000.0
step_error                    2.8571428571428573e-18
est_needed_step_err           2.8571428571428573e-18
hn_div_ho                     0.5
hn_div_ho_2                   0.25
hn_div_ho_3                   0.125
value3                        7.130136767447487e-204
max_value3                    7.130136767447487e-204
value3                        7.130136767447487e-204
best_h                        1.0e-07
START of Soultion
 
TOP MAIN SOLVE Loop
t[1]                             1.5 
x2[1] (analytic)                 1.0040840464326821 
x2[1] (numeric)                  1.0040840464326821 
absolute error                   0.0 
relative error                   0.0%
Correct digits                   16 
h                                1.0e-07 
x1[1] (analytic)                 2.0004016342882673 
x1[1] (numeric)                  2.0004016342882673 
absolute error                   0.0 
relative error                   0.0%
Correct digits                   16 
h                                1.0e-07 
Complex estimate of poles used for equation 1
NO POLE for equation 2
Radius of convergence            32770.8950325607 
Order of pole                    59662864.67507763 
Finished!
Maximum Time Reached before Solution Completed!
diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1)+x1;
diff (x1,t,1) = 4.0 * x2 -  2.0 * diff (x2,t ,1) - 2.0 * x1;
Iterations                       4859 
Total Elapsed Time              3 Minutes 0.0 Seconds
Elapsed Time(since restart)     2 Minutes 59.0 Seconds
Expected Time Remaining         15 Days 0 Hours 1 Minutes 56.29553943150677 Seconds
Optimized Time Remaining        14 Days 22 Hours 1 Minutes 55.64945310121402 Seconds
Expected Total Time             14 Days 22 Hours 4 Minutes 55.64945310121402 Seconds
Time to Timeout                 0.0 Seconds
Percent Done                    0.013885714293821698%