(%i1) batch(diffeq.max)
read and interpret file: /home/dennis/mastersource/mine/omnisode/diffeq.max
(%i2) load(stringproc)
(%o2) /usr/share/maxima/5.27.0/share/stringproc/stringproc.mac
(%i3) check_sign(x0, xf) := block([ret],
if xf > x0 then ret : 1.0 else ret : - 1.0, ret)
(%o3) check_sign(x0, xf) := block([ret],
if xf > x0 then ret : 1.0 else ret : - 1.0, ret)
(%i4) est_size_answer() := block([min_size], min_size : glob_large_float,
if omniabs(array_y ) < min_size then (min_size : omniabs(array_y ),
1 1
omniout_float(ALWAYS, "min_size", 32, min_size, 32, "")),
if min_size < 1.0 then (min_size : 1.0,
omniout_float(ALWAYS, "min_size", 32, min_size, 32, "")), min_size)
(%o4) est_size_answer() := block([min_size], min_size : glob_large_float,
if omniabs(array_y ) < min_size then (min_size : omniabs(array_y ),
1 1
omniout_float(ALWAYS, "min_size", 32, min_size, 32, "")),
if min_size < 1.0 then (min_size : 1.0,
omniout_float(ALWAYS, "min_size", 32, min_size, 32, "")), min_size)
(%i5) test_suggested_h() := block([max_value3, hn_div_ho, hn_div_ho_2,
hn_div_ho_3, value3, no_terms], max_value3 : 0.0, no_terms : glob_max_terms,
hn_div_ho : 0.5, hn_div_ho_2 : 0.25, hn_div_ho_3 : 0.125,
omniout_float(ALWAYS, "hn_div_ho", 32, hn_div_ho, 32, ""),
omniout_float(ALWAYS, "hn_div_ho_2", 32, hn_div_ho_2, 32, ""),
omniout_float(ALWAYS, "hn_div_ho_3", 32, hn_div_ho_3, 32, ""),
value3 : omniabs(array_y hn_div_ho_3 + array_y hn_div_ho_2
no_terms no_terms - 1
+ array_y hn_div_ho + array_y ),
no_terms - 2 no_terms - 3
if value3 > max_value3 then (max_value3 : value3,
omniout_float(ALWAYS, "value3", 32, value3, 32, "")),
omniout_float(ALWAYS, "max_value3", 32, max_value3, 32, ""), max_value3)
(%o5) test_suggested_h() := block([max_value3, hn_div_ho, hn_div_ho_2,
hn_div_ho_3, value3, no_terms], max_value3 : 0.0, no_terms : glob_max_terms,
hn_div_ho : 0.5, hn_div_ho_2 : 0.25, hn_div_ho_3 : 0.125,
omniout_float(ALWAYS, "hn_div_ho", 32, hn_div_ho, 32, ""),
omniout_float(ALWAYS, "hn_div_ho_2", 32, hn_div_ho_2, 32, ""),
omniout_float(ALWAYS, "hn_div_ho_3", 32, hn_div_ho_3, 32, ""),
value3 : omniabs(array_y hn_div_ho_3 + array_y hn_div_ho_2
no_terms no_terms - 1
+ array_y hn_div_ho + array_y ),
no_terms - 2 no_terms - 3
if value3 > max_value3 then (max_value3 : value3,
omniout_float(ALWAYS, "value3", 32, value3, 32, "")),
omniout_float(ALWAYS, "max_value3", 32, max_value3, 32, ""), max_value3)
(%i6) reached_interval() := block([ret],
if glob_check_sign array_x >= glob_check_sign glob_next_display
1
then ret : true else ret : false, return(ret))
(%o6) reached_interval() := block([ret],
if glob_check_sign array_x >= glob_check_sign glob_next_display
1
then ret : true else ret : false, return(ret))
(%i7) display_alot(iter) := block([abserr, analytic_val_y, ind_var,
numeric_val, relerr, term_no], if reached_interval()
then (if iter >= 0 then (ind_var : array_x ,
1
omniout_float(ALWAYS, "x[1] ", 33, ind_var, 20,
" "), analytic_val_y : exact_soln_y(ind_var),
omniout_float(ALWAYS, "y[1] (analytic) ", 33, analytic_val_y,
20, " "), term_no : 1, numeric_val : array_y ,
term_no
abserr : omniabs(numeric_val - analytic_val_y),
omniout_float(ALWAYS, "y[1] (numeric) ", 33, numeric_val,
20, " "), if omniabs(analytic_val_y) # 0.0
abserr 100.0
then (relerr : -----------------------,
omniabs(analytic_val_y)
if relerr > 1.0E-34 then glob_good_digits : 2 - floor(log10(relerr))
else glob_good_digits : 16) else (relerr : - 1.0, glob_good_digits : - 1),
if glob_iter = 1 then array_1st_rel_error : relerr
1
else array_last_rel_error : relerr, omniout_float(ALWAYS,
1
"absolute error ", 4, abserr, 20, " "),
omniout_float(ALWAYS, "relative error ", 4, relerr, 20,
"%"), omniout_int(INFO, "Correct digits ", 32,
glob_good_digits, 4, " "), omniout_float(ALWAYS,
"h ", 4, glob_h, 20, " "))))
(%o7) display_alot(iter) := block([abserr, analytic_val_y, ind_var,
numeric_val, relerr, term_no], if reached_interval()
then (if iter >= 0 then (ind_var : array_x ,
1
omniout_float(ALWAYS, "x[1] ", 33, ind_var, 20,
" "), analytic_val_y : exact_soln_y(ind_var),
omniout_float(ALWAYS, "y[1] (analytic) ", 33, analytic_val_y,
20, " "), term_no : 1, numeric_val : array_y ,
term_no
abserr : omniabs(numeric_val - analytic_val_y),
omniout_float(ALWAYS, "y[1] (numeric) ", 33, numeric_val,
20, " "), if omniabs(analytic_val_y) # 0.0
abserr 100.0
then (relerr : -----------------------,
omniabs(analytic_val_y)
if relerr > 1.0E-34 then glob_good_digits : 2 - floor(log10(relerr))
else glob_good_digits : 16) else (relerr : - 1.0, glob_good_digits : - 1),
if glob_iter = 1 then array_1st_rel_error : relerr
1
else array_last_rel_error : relerr, omniout_float(ALWAYS,
1
"absolute error ", 4, abserr, 20, " "),
omniout_float(ALWAYS, "relative error ", 4, relerr, 20,
"%"), omniout_int(INFO, "Correct digits ", 32,
glob_good_digits, 4, " "), omniout_float(ALWAYS,
"h ", 4, glob_h, 20, " "))))
(%i8) adjust_for_pole(h_param) := block([hnew, sz2, tmp],
block(hnew : h_param, glob_normmax : glob_small_float,
if omniabs(array_y_higher ) > glob_small_float
1, 1
then (tmp : omniabs(array_y_higher ),
1, 1
if tmp < glob_normmax then glob_normmax : tmp),
if glob_look_poles and (omniabs(array_pole ) > glob_small_float)
1
array_pole
1
and (array_pole # glob_large_float) then (sz2 : -----------,
1 10.0
if sz2 < hnew then (omniout_float(INFO, "glob_h adjusted to ", 20, h_param,
12, "due to singularity."), omniout_str(INFO, "Reached Optimal"),
return(hnew))), if not glob_reached_optimal_h
then (glob_reached_optimal_h : true, glob_curr_iter_when_opt :
glob_current_iter, glob_optimal_clock_start_sec : elapsed_time_seconds(),
glob_optimal_start : array_x ), hnew : sz2), return(hnew))
1
(%o8) adjust_for_pole(h_param) := block([hnew, sz2, tmp],
block(hnew : h_param, glob_normmax : glob_small_float,
if omniabs(array_y_higher ) > glob_small_float
1, 1
then (tmp : omniabs(array_y_higher ),
1, 1
if tmp < glob_normmax then glob_normmax : tmp),
if glob_look_poles and (omniabs(array_pole ) > glob_small_float)
1
array_pole
1
and (array_pole # glob_large_float) then (sz2 : -----------,
1 10.0
if sz2 < hnew then (omniout_float(INFO, "glob_h adjusted to ", 20, h_param,
12, "due to singularity."), omniout_str(INFO, "Reached Optimal"),
return(hnew))), if not glob_reached_optimal_h
then (glob_reached_optimal_h : true, glob_curr_iter_when_opt :
glob_current_iter, glob_optimal_clock_start_sec : elapsed_time_seconds(),
glob_optimal_start : array_x ), hnew : sz2), return(hnew))
1
(%i9) prog_report(x_start, x_end) := block([clock_sec, opt_clock_sec,
clock_sec1, expect_sec, left_sec, percent_done, total_clock_sec],
clock_sec1 : elapsed_time_seconds(), total_clock_sec :
convfloat(clock_sec1) - convfloat(glob_orig_start_sec),
glob_clock_sec : convfloat(clock_sec1) - convfloat(glob_clock_start_sec),
left_sec : - convfloat(clock_sec1) + convfloat(glob_orig_start_sec)
+ convfloat(glob_max_sec), expect_sec :
comp_expect_sec(convfloat(x_end), convfloat(x_start),
convfloat(glob_h) + convfloat(array_x ),
1
convfloat(clock_sec1) - convfloat(glob_orig_start_sec)),
opt_clock_sec : convfloat(clock_sec1)
- convfloat(glob_optimal_clock_start_sec),
glob_optimal_expect_sec : comp_expect_sec(convfloat(x_end),
convfloat(x_start), convfloat(glob_h) + convfloat(array_x ),
1
convfloat(opt_clock_sec)), glob_total_exp_sec :
total_clock_sec + glob_optimal_expect_sec,
percent_done : comp_percent(convfloat(x_end), convfloat(x_start),
convfloat(glob_h) + convfloat(array_x )), glob_percent_done : percent_done,
1
omniout_str_noeol(INFO, "Total Elapsed Time "),
omniout_timestr(convfloat(total_clock_sec)),
omniout_str_noeol(INFO, "Elapsed Time(since restart) "),
omniout_timestr(convfloat(glob_clock_sec)),
if convfloat(percent_done) < convfloat(100.0)
then (omniout_str_noeol(INFO, "Expected Time Remaining "),
omniout_timestr(convfloat(expect_sec)),
omniout_str_noeol(INFO, "Optimized Time Remaining "),
omniout_timestr(convfloat(glob_optimal_expect_sec)),
omniout_str_noeol(INFO, "Expected Total Time "),
omniout_timestr(convfloat(glob_total_exp_sec))),
omniout_str_noeol(INFO, "Time to Timeout "),
omniout_timestr(convfloat(left_sec)), omniout_float(INFO,
"Percent Done ", 33, percent_done, 4, "%"))
(%o9) prog_report(x_start, x_end) := block([clock_sec, opt_clock_sec,
clock_sec1, expect_sec, left_sec, percent_done, total_clock_sec],
clock_sec1 : elapsed_time_seconds(), total_clock_sec :
convfloat(clock_sec1) - convfloat(glob_orig_start_sec),
glob_clock_sec : convfloat(clock_sec1) - convfloat(glob_clock_start_sec),
left_sec : - convfloat(clock_sec1) + convfloat(glob_orig_start_sec)
+ convfloat(glob_max_sec), expect_sec :
comp_expect_sec(convfloat(x_end), convfloat(x_start),
convfloat(glob_h) + convfloat(array_x ),
1
convfloat(clock_sec1) - convfloat(glob_orig_start_sec)),
opt_clock_sec : convfloat(clock_sec1)
- convfloat(glob_optimal_clock_start_sec),
glob_optimal_expect_sec : comp_expect_sec(convfloat(x_end),
convfloat(x_start), convfloat(glob_h) + convfloat(array_x ),
1
convfloat(opt_clock_sec)), glob_total_exp_sec :
total_clock_sec + glob_optimal_expect_sec,
percent_done : comp_percent(convfloat(x_end), convfloat(x_start),
convfloat(glob_h) + convfloat(array_x )), glob_percent_done : percent_done,
1
omniout_str_noeol(INFO, "Total Elapsed Time "),
omniout_timestr(convfloat(total_clock_sec)),
omniout_str_noeol(INFO, "Elapsed Time(since restart) "),
omniout_timestr(convfloat(glob_clock_sec)),
if convfloat(percent_done) < convfloat(100.0)
then (omniout_str_noeol(INFO, "Expected Time Remaining "),
omniout_timestr(convfloat(expect_sec)),
omniout_str_noeol(INFO, "Optimized Time Remaining "),
omniout_timestr(convfloat(glob_optimal_expect_sec)),
omniout_str_noeol(INFO, "Expected Total Time "),
omniout_timestr(convfloat(glob_total_exp_sec))),
omniout_str_noeol(INFO, "Time to Timeout "),
omniout_timestr(convfloat(left_sec)), omniout_float(INFO,
"Percent Done ", 33, percent_done, 4, "%"))
(%i10) check_for_pole() := block([cnt, dr1, dr2, ds1, ds2, hdrc, m, n, nr1,
nr2, ord_no, rad_c, rcs, rm0, rm1, rm2, rm3, rm4, found, h_new, ratio, term],
n : glob_max_terms, m : - 1 - 1 + n, while (m >= 10)
and ((omniabs(array_y_higher ) < glob_small_float)
1, m
or (omniabs(array_y_higher ) < glob_small_float)
1, m - 1
or (omniabs(array_y_higher ) < glob_small_float)) do m :
1, m - 2
array_y_higher
1, m
m - 1, if m > 10 then (rm0 : ----------------------,
array_y_higher
1, m - 1
array_y_higher
1, m - 1
rm1 : ----------------------, hdrc : convfloat(m - 1) rm0
array_y_higher
1, m - 2
- convfloat(m - 2) rm1, if omniabs(hdrc) > glob_small_float
glob_h convfloat(m - 1) rm0
then (rcs : ------, ord_no : 2.0 - convfloat(m) + --------------------,
hdrc hdrc
array_real_pole : rcs, array_real_pole : ord_no)
1, 1 1, 2
else (array_real_pole : glob_large_float,
1, 1
array_real_pole : glob_large_float))
1, 2
else (array_real_pole : glob_large_float,
1, 1
array_real_pole : glob_large_float), n : - 1 - 1 + glob_max_terms,
1, 2
cnt : 0, while (cnt < 5) and (n >= 10) do (if omniabs(array_y_higher ) >
1, n
glob_small_float then cnt : 1 + cnt else cnt : 0, n : n - 1), m : cnt + n,
if m <= 10 then (rad_c : glob_large_float, ord_no : glob_large_float)
elseif ((omniabs(array_y_higher ) >= glob_large_float)
1, m
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 1
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 2
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 3
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 4
or (omniabs(array_y_higher ) >= glob_large_float))
1, m - 5
or ((omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float))
1, m 1, m - 1 1, m - 2 1, m - 3 1, m - 4 1, m - 5
then (rad_c : glob_large_float, ord_no : glob_large_float)
array_y_higher array_y_higher
1, m 1, m - 1
else (rm0 : ----------------------, rm1 : ----------------------,
array_y_higher array_y_higher
1, m - 1 1, m - 2
array_y_higher array_y_higher
1, m - 2 1, m - 3
rm2 : ----------------------, rm3 : ----------------------,
array_y_higher array_y_higher
1, m - 3 1, m - 4
array_y_higher
1, m - 4
rm4 : ----------------------, nr1 : convfloat(m - 3) rm2
array_y_higher
1, m - 5
- 2.0 convfloat(m - 2) rm1 + convfloat(m - 1) rm0,
nr2 : convfloat(m - 4) rm3 - 2.0 convfloat(m - 3) rm2 + convfloat(m - 2) rm1,
- 1.0 2.0 - 1.0 - 1.0 2.0 - 1.0 5.0 8.0 3.0
dr1 : ----- + --- + -----, dr2 : ----- + --- + -----, ds1 : --- - --- + ---,
rm3 rm2 rm1 rm4 rm3 rm2 rm3 rm2 rm1
5.0 8.0 3.0
ds2 : --- - --- + ---, if (omniabs(nr1 dr2 - nr2 dr1) <= glob_small_float)
rm4 rm3 rm2
or (omniabs(dr1) <= glob_small_float) then (rad_c : glob_large_float,
ord_no : glob_large_float) else (if omniabs(nr1 dr2 - nr2 dr1) >
dr1 dr2 - ds2 dr1 + ds1 dr2
glob_small_float then (rcs : ---------------------------,
nr1 dr2 - nr2 dr1
rcs nr1 - ds1 convfloat(m)
ord_no : ------------- - ------------,
2.0 dr1 2.0
if omniabs(rcs) > glob_small_float then (if rcs > 0.0
then rad_c : sqrt(rcs) omniabs(glob_h) else rad_c : glob_large_float)
else (rad_c : glob_large_float, ord_no : glob_large_float))
else (rad_c : glob_large_float, ord_no : glob_large_float)),
array_complex_pole : rad_c, array_complex_pole : ord_no),
1, 1 1, 2
found : false, if (not found) and ((array_real_pole = glob_large_float)
1, 1
or (array_real_pole = glob_large_float))
1, 2
and ((array_complex_pole # glob_large_float) and (array_complex_pole # glob_large_float))
1, 1 1, 2
and ((array_complex_pole > 0.0) and (array_complex_pole > 0.0))
1, 1 1, 2
then (array_poles : array_complex_pole ,
1, 1 1, 1
array_poles : array_complex_pole , found : true, array_type_pole : 2,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Complex estimate of poles used"))),
if (not found) and ((array_real_pole # glob_large_float)
1, 1
and (array_real_pole # glob_large_float) and (array_real_pole > 0.0)
1, 2 1, 1
and (array_real_pole > 0.0) and ((array_complex_pole = glob_large_float) or (array_complex_pole = glob_large_float) or (array_complex_pole <= 0.0) or (array_complex_pole <=
1, 2 1, 1 1, 2 1, 1 1, 2
0.0)))
then (array_poles : array_real_pole ,
1, 1 1, 1
array_poles : array_real_pole , found : true, array_type_pole : 1,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Real estimate of pole used"))),
if (not found) and (((array_real_pole = glob_large_float)
1, 1
or (array_real_pole = glob_large_float))
1, 2
and ((array_complex_pole = glob_large_float) or (array_complex_pole = glob_large_float)))
1, 1 1, 2
then (array_poles : glob_large_float, array_poles : glob_large_float,
1, 1 1, 2
found : true, array_type_pole : 3, if reached_interval()
1
then omniout_str(ALWAYS, "NO POLE")),
if (not found) and ((array_real_pole < array_complex_pole )
1, 1 1, 1
and (array_real_pole > 0.0) and (array_real_pole >
1, 1 1, 2
0.0))
then (array_poles : array_real_pole ,
1, 1 1, 1
array_poles : array_real_pole , found : true, array_type_pole : 1,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Real estimate of pole used"))),
if (not found) and ((array_complex_pole # glob_large_float)
1, 1
and (array_complex_pole # glob_large_float)
1, 2
and (array_complex_pole > 0.0) and (array_complex_pole >
1, 1 1, 2
0.0))
then (array_poles : array_complex_pole ,
1, 1 1, 1
array_poles : array_complex_pole , array_type_pole : 2, found : true,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Complex estimate of poles used"))),
if not found then (array_poles : glob_large_float,
1, 1
array_poles : glob_large_float, array_type_pole : 3,
1, 2 1
if reached_interval() then omniout_str(ALWAYS, "NO POLE")),
array_pole : glob_large_float, array_pole : glob_large_float,
1 2
if array_pole > array_poles then (array_pole : array_poles ,
1 1, 1 1 1, 1
array_pole : array_poles ), if array_pole glob_ratio_of_radius <
2 1, 2 1
omniabs(glob_h) then (h_new : array_pole glob_ratio_of_radius, term : 1,
1
ratio : 1.0, while term <= glob_max_terms do (array_y :
term
array_y ratio, array_y_higher : array_y_higher ratio,
term 1, term 1, term
ratio h_new
array_x : array_x ratio, ratio : ---------------, term : 1 + term),
term term omniabs(glob_h)
glob_h : h_new), if reached_interval() then display_pole())
(%o10) check_for_pole() := block([cnt, dr1, dr2, ds1, ds2, hdrc, m, n, nr1,
nr2, ord_no, rad_c, rcs, rm0, rm1, rm2, rm3, rm4, found, h_new, ratio, term],
n : glob_max_terms, m : - 1 - 1 + n, while (m >= 10)
and ((omniabs(array_y_higher ) < glob_small_float)
1, m
or (omniabs(array_y_higher ) < glob_small_float)
1, m - 1
or (omniabs(array_y_higher ) < glob_small_float)) do m :
1, m - 2
array_y_higher
1, m
m - 1, if m > 10 then (rm0 : ----------------------,
array_y_higher
1, m - 1
array_y_higher
1, m - 1
rm1 : ----------------------, hdrc : convfloat(m - 1) rm0
array_y_higher
1, m - 2
- convfloat(m - 2) rm1, if omniabs(hdrc) > glob_small_float
glob_h convfloat(m - 1) rm0
then (rcs : ------, ord_no : 2.0 - convfloat(m) + --------------------,
hdrc hdrc
array_real_pole : rcs, array_real_pole : ord_no)
1, 1 1, 2
else (array_real_pole : glob_large_float,
1, 1
array_real_pole : glob_large_float))
1, 2
else (array_real_pole : glob_large_float,
1, 1
array_real_pole : glob_large_float), n : - 1 - 1 + glob_max_terms,
1, 2
cnt : 0, while (cnt < 5) and (n >= 10) do (if omniabs(array_y_higher ) >
1, n
glob_small_float then cnt : 1 + cnt else cnt : 0, n : n - 1), m : cnt + n,
if m <= 10 then (rad_c : glob_large_float, ord_no : glob_large_float)
elseif ((omniabs(array_y_higher ) >= glob_large_float)
1, m
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 1
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 2
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 3
or (omniabs(array_y_higher ) >= glob_large_float)
1, m - 4
or (omniabs(array_y_higher ) >= glob_large_float))
1, m - 5
or ((omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float) or (omniabs(array_y_higher ) <= glob_small_float))
1, m 1, m - 1 1, m - 2 1, m - 3 1, m - 4 1, m - 5
then (rad_c : glob_large_float, ord_no : glob_large_float)
array_y_higher array_y_higher
1, m 1, m - 1
else (rm0 : ----------------------, rm1 : ----------------------,
array_y_higher array_y_higher
1, m - 1 1, m - 2
array_y_higher array_y_higher
1, m - 2 1, m - 3
rm2 : ----------------------, rm3 : ----------------------,
array_y_higher array_y_higher
1, m - 3 1, m - 4
array_y_higher
1, m - 4
rm4 : ----------------------, nr1 : convfloat(m - 3) rm2
array_y_higher
1, m - 5
- 2.0 convfloat(m - 2) rm1 + convfloat(m - 1) rm0,
nr2 : convfloat(m - 4) rm3 - 2.0 convfloat(m - 3) rm2 + convfloat(m - 2) rm1,
- 1.0 2.0 - 1.0 - 1.0 2.0 - 1.0 5.0 8.0 3.0
dr1 : ----- + --- + -----, dr2 : ----- + --- + -----, ds1 : --- - --- + ---,
rm3 rm2 rm1 rm4 rm3 rm2 rm3 rm2 rm1
5.0 8.0 3.0
ds2 : --- - --- + ---, if (omniabs(nr1 dr2 - nr2 dr1) <= glob_small_float)
rm4 rm3 rm2
or (omniabs(dr1) <= glob_small_float) then (rad_c : glob_large_float,
ord_no : glob_large_float) else (if omniabs(nr1 dr2 - nr2 dr1) >
dr1 dr2 - ds2 dr1 + ds1 dr2
glob_small_float then (rcs : ---------------------------,
nr1 dr2 - nr2 dr1
rcs nr1 - ds1 convfloat(m)
ord_no : ------------- - ------------,
2.0 dr1 2.0
if omniabs(rcs) > glob_small_float then (if rcs > 0.0
then rad_c : sqrt(rcs) omniabs(glob_h) else rad_c : glob_large_float)
else (rad_c : glob_large_float, ord_no : glob_large_float))
else (rad_c : glob_large_float, ord_no : glob_large_float)),
array_complex_pole : rad_c, array_complex_pole : ord_no),
1, 1 1, 2
found : false, if (not found) and ((array_real_pole = glob_large_float)
1, 1
or (array_real_pole = glob_large_float))
1, 2
and ((array_complex_pole # glob_large_float) and (array_complex_pole # glob_large_float))
1, 1 1, 2
and ((array_complex_pole > 0.0) and (array_complex_pole > 0.0))
1, 1 1, 2
then (array_poles : array_complex_pole ,
1, 1 1, 1
array_poles : array_complex_pole , found : true, array_type_pole : 2,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Complex estimate of poles used"))),
if (not found) and ((array_real_pole # glob_large_float)
1, 1
and (array_real_pole # glob_large_float) and (array_real_pole > 0.0)
1, 2 1, 1
and (array_real_pole > 0.0) and ((array_complex_pole = glob_large_float) or (array_complex_pole = glob_large_float) or (array_complex_pole <= 0.0) or (array_complex_pole <=
1, 2 1, 1 1, 2 1, 1 1, 2
0.0)))
then (array_poles : array_real_pole ,
1, 1 1, 1
array_poles : array_real_pole , found : true, array_type_pole : 1,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Real estimate of pole used"))),
if (not found) and (((array_real_pole = glob_large_float)
1, 1
or (array_real_pole = glob_large_float))
1, 2
and ((array_complex_pole = glob_large_float) or (array_complex_pole = glob_large_float)))
1, 1 1, 2
then (array_poles : glob_large_float, array_poles : glob_large_float,
1, 1 1, 2
found : true, array_type_pole : 3, if reached_interval()
1
then omniout_str(ALWAYS, "NO POLE")),
if (not found) and ((array_real_pole < array_complex_pole )
1, 1 1, 1
and (array_real_pole > 0.0) and (array_real_pole >
1, 1 1, 2
0.0))
then (array_poles : array_real_pole ,
1, 1 1, 1
array_poles : array_real_pole , found : true, array_type_pole : 1,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Real estimate of pole used"))),
if (not found) and ((array_complex_pole # glob_large_float)
1, 1
and (array_complex_pole # glob_large_float)
1, 2
and (array_complex_pole > 0.0) and (array_complex_pole >
1, 1 1, 2
0.0))
then (array_poles : array_complex_pole ,
1, 1 1, 1
array_poles : array_complex_pole , array_type_pole : 2, found : true,
1, 2 1, 2 1
if glob_display_flag then (if reached_interval()
then omniout_str(ALWAYS, "Complex estimate of poles used"))),
if not found then (array_poles : glob_large_float,
1, 1
array_poles : glob_large_float, array_type_pole : 3,
1, 2 1
if reached_interval() then omniout_str(ALWAYS, "NO POLE")),
array_pole : glob_large_float, array_pole : glob_large_float,
1 2
if array_pole > array_poles then (array_pole : array_poles ,
1 1, 1 1 1, 1
array_pole : array_poles ), if array_pole glob_ratio_of_radius <
2 1, 2 1
omniabs(glob_h) then (h_new : array_pole glob_ratio_of_radius, term : 1,
1
ratio : 1.0, while term <= glob_max_terms do (array_y :
term
array_y ratio, array_y_higher : array_y_higher ratio,
term 1, term 1, term
ratio h_new
array_x : array_x ratio, ratio : ---------------, term : 1 + term),
term term omniabs(glob_h)
glob_h : h_new), if reached_interval() then display_pole())
(%i11) get_norms() := block([iii], if not glob_initial_pass
then (iii : 1, while iii <= glob_max_terms do (array_norms : 0.0,
iii
iii : 1 + iii), iii : 1, while iii <=
glob_max_terms do (if omniabs(array_y ) > array_norms
iii iii
then array_norms : omniabs(array_y ), iii : 1 + iii)))
iii iii
(%o11) get_norms() := block([iii], if not glob_initial_pass
then (iii : 1, while iii <= glob_max_terms do (array_norms : 0.0,
iii
iii : 1 + iii), iii : 1, while iii <=
glob_max_terms do (if omniabs(array_y ) > array_norms
iii iii
then array_norms : omniabs(array_y ), iii : 1 + iii)))
iii iii
(%i12) atomall() := block([kkk, order_d, adj2, adj3, temporary, term, temp,
temp2], array_tmp1 : array_const_3D0 array_x ,
1 1 1
array_tmp2 : array_const_1D0 + array_tmp1 ,
1 1 1
array_tmp3_a1 : sinh(array_tmp2 ), array_tmp3_a2 : cosh(array_tmp2 ),
1 1 1 1
array_tmp3_a1
1
array_tmp3 : --------------, array_tmp4 : array_tmp3 + array_const_0D0 ,
1 array_tmp3_a2 1 1 1
1
if not array_y_set_initial then (if 1 <= glob_max_terms
1, 2
then (temporary : array_tmp4 expt(glob_h, 1) factorial_3(0, 1),
1
array_y : temporary, array_y_higher : temporary,
2 1, 2
temporary 1.0
temporary : -------------, array_y_higher : temporary, 0)), kkk : 2,
glob_h 2, 1
array_tmp1 : array_const_3D0 array_x , array_tmp2 : array_tmp1 ,
2 1 2 2 2
array_tmp3_a2 array_tmp2
1 2
array_tmp3_a1 : --------------------------,
2 1
array_tmp3_a1 array_tmp2
1 2
array_tmp3_a2 : --------------------------,
2 1
array_tmp3_a1 - ats(2, array_tmp3_a2, array_tmp3, 2)
2
array_tmp3 : -----------------------------------------------------,
2 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
2 2 1, 3
then (if 2 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(1, 2), array_y : temporary,
2 3
temporary 2.0
array_y_higher : temporary, temporary : -------------,
1, 3 glob_h
array_y_higher : temporary, 0)), kkk : 3,
2, 2
array_tmp3_a2 array_tmp2
2 2
array_tmp3_a1 : --------------------------,
3 2
array_tmp3_a1 array_tmp2
2 2
array_tmp3_a2 : --------------------------,
3 2
array_tmp3_a1 - ats(3, array_tmp3_a2, array_tmp3, 2)
3
array_tmp3 : -----------------------------------------------------,
3 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
3 3 1, 4
then (if 3 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(2, 3), array_y : temporary,
3 4
temporary 3.0
array_y_higher : temporary, temporary : -------------,
1, 4 glob_h
array_y_higher : temporary, 0)), kkk : 4,
2, 3
array_tmp3_a2 array_tmp2
3 2
array_tmp3_a1 : --------------------------,
4 3
array_tmp3_a1 array_tmp2
3 2
array_tmp3_a2 : --------------------------,
4 3
array_tmp3_a1 - ats(4, array_tmp3_a2, array_tmp3, 2)
4
array_tmp3 : -----------------------------------------------------,
4 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
4 4 1, 5
then (if 4 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(3, 4), array_y : temporary,
4 5
temporary 4.0
array_y_higher : temporary, temporary : -------------,
1, 5 glob_h
array_y_higher : temporary, 0)), kkk : 5,
2, 4
array_tmp3_a2 array_tmp2
4 2
array_tmp3_a1 : --------------------------,
5 4
array_tmp3_a1 array_tmp2
4 2
array_tmp3_a2 : --------------------------,
5 4
array_tmp3_a1 - ats(5, array_tmp3_a2, array_tmp3, 2)
5
array_tmp3 : -----------------------------------------------------,
5 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
5 5 1, 6
then (if 5 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(4, 5), array_y : temporary,
5 6
temporary 5.0
array_y_higher : temporary, temporary : -------------,
1, 6 glob_h
array_y_higher : temporary, 0)), kkk : 6,
2, 5
while kkk <= glob_max_terms do (array_tmp3_a1 :
kkk
array_tmp3_a2 array_tmp2
kkk - 1 2
--------------------------------, array_tmp3_a2 :
kkk - 1 kkk
array_tmp3_a1 array_tmp2
kkk - 1 2
--------------------------------, array_tmp3 :
kkk - 1 kkk
array_tmp3_a1 - ats(kkk, array_tmp3_a2, array_tmp3, 2)
kkk
---------------------------------------------------------,
array_tmp3_a2
1
array_tmp4 : array_tmp3 , order_d : 1,
kkk kkk
if 1 + order_d + kkk <= glob_max_terms
then (if not array_y_set_initial
1, order_d + kkk
then (temporary : array_tmp4 expt(glob_h, order_d)
kkk
factorial_3(kkk - 1, - 1 + order_d + kkk), array_y : temporary,
order_d + kkk
array_y_higher : temporary, term : - 1 + order_d + kkk,
1, order_d + kkk
adj2 : - 1 + order_d + kkk, adj3 : 2, while term >=
1 do (if adj3 <= 1 + order_d then (if adj2 > 0
temporary convfp(adj2)
then temporary : ---------------------- else temporary : temporary,
glob_h
array_y_higher : temporary), term : term - 1, adj2 : adj2 - 1,
adj3, term
adj3 : 1 + adj3))), kkk : 1 + kkk))
(%o12) atomall() := block([kkk, order_d, adj2, adj3, temporary, term, temp,
temp2], array_tmp1 : array_const_3D0 array_x ,
1 1 1
array_tmp2 : array_const_1D0 + array_tmp1 ,
1 1 1
array_tmp3_a1 : sinh(array_tmp2 ), array_tmp3_a2 : cosh(array_tmp2 ),
1 1 1 1
array_tmp3_a1
1
array_tmp3 : --------------, array_tmp4 : array_tmp3 + array_const_0D0 ,
1 array_tmp3_a2 1 1 1
1
if not array_y_set_initial then (if 1 <= glob_max_terms
1, 2
then (temporary : array_tmp4 expt(glob_h, 1) factorial_3(0, 1),
1
array_y : temporary, array_y_higher : temporary,
2 1, 2
temporary 1.0
temporary : -------------, array_y_higher : temporary, 0)), kkk : 2,
glob_h 2, 1
array_tmp1 : array_const_3D0 array_x , array_tmp2 : array_tmp1 ,
2 1 2 2 2
array_tmp3_a2 array_tmp2
1 2
array_tmp3_a1 : --------------------------,
2 1
array_tmp3_a1 array_tmp2
1 2
array_tmp3_a2 : --------------------------,
2 1
array_tmp3_a1 - ats(2, array_tmp3_a2, array_tmp3, 2)
2
array_tmp3 : -----------------------------------------------------,
2 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
2 2 1, 3
then (if 2 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(1, 2), array_y : temporary,
2 3
temporary 2.0
array_y_higher : temporary, temporary : -------------,
1, 3 glob_h
array_y_higher : temporary, 0)), kkk : 3,
2, 2
array_tmp3_a2 array_tmp2
2 2
array_tmp3_a1 : --------------------------,
3 2
array_tmp3_a1 array_tmp2
2 2
array_tmp3_a2 : --------------------------,
3 2
array_tmp3_a1 - ats(3, array_tmp3_a2, array_tmp3, 2)
3
array_tmp3 : -----------------------------------------------------,
3 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
3 3 1, 4
then (if 3 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(2, 3), array_y : temporary,
3 4
temporary 3.0
array_y_higher : temporary, temporary : -------------,
1, 4 glob_h
array_y_higher : temporary, 0)), kkk : 4,
2, 3
array_tmp3_a2 array_tmp2
3 2
array_tmp3_a1 : --------------------------,
4 3
array_tmp3_a1 array_tmp2
3 2
array_tmp3_a2 : --------------------------,
4 3
array_tmp3_a1 - ats(4, array_tmp3_a2, array_tmp3, 2)
4
array_tmp3 : -----------------------------------------------------,
4 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
4 4 1, 5
then (if 4 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(3, 4), array_y : temporary,
4 5
temporary 4.0
array_y_higher : temporary, temporary : -------------,
1, 5 glob_h
array_y_higher : temporary, 0)), kkk : 5,
2, 4
array_tmp3_a2 array_tmp2
4 2
array_tmp3_a1 : --------------------------,
5 4
array_tmp3_a1 array_tmp2
4 2
array_tmp3_a2 : --------------------------,
5 4
array_tmp3_a1 - ats(5, array_tmp3_a2, array_tmp3, 2)
5
array_tmp3 : -----------------------------------------------------,
5 array_tmp3_a2
1
array_tmp4 : array_tmp3 , if not array_y_set_initial
5 5 1, 6
then (if 5 <= glob_max_terms then (temporary :
array_tmp4 expt(glob_h, 1) factorial_3(4, 5), array_y : temporary,
5 6
temporary 5.0
array_y_higher : temporary, temporary : -------------,
1, 6 glob_h
array_y_higher : temporary, 0)), kkk : 6,
2, 5
while kkk <= glob_max_terms do (array_tmp3_a1 :
kkk
array_tmp3_a2 array_tmp2
kkk - 1 2
--------------------------------, array_tmp3_a2 :
kkk - 1 kkk
array_tmp3_a1 array_tmp2
kkk - 1 2
--------------------------------, array_tmp3 :
kkk - 1 kkk
array_tmp3_a1 - ats(kkk, array_tmp3_a2, array_tmp3, 2)
kkk
---------------------------------------------------------,
array_tmp3_a2
1
array_tmp4 : array_tmp3 , order_d : 1,
kkk kkk
if 1 + order_d + kkk <= glob_max_terms
then (if not array_y_set_initial
1, order_d + kkk
then (temporary : array_tmp4 expt(glob_h, order_d)
kkk
factorial_3(kkk - 1, - 1 + order_d + kkk), array_y : temporary,
order_d + kkk
array_y_higher : temporary, term : - 1 + order_d + kkk,
1, order_d + kkk
adj2 : - 1 + order_d + kkk, adj3 : 2, while term >=
1 do (if adj3 <= 1 + order_d then (if adj2 > 0
temporary convfp(adj2)
then temporary : ---------------------- else temporary : temporary,
glob_h
array_y_higher : temporary), term : term - 1, adj2 : adj2 - 1,
adj3, term
adj3 : 1 + adj3))), kkk : 1 + kkk))
log(x)
(%i13) log10(x) := ---------
log(10.0)
log(x)
(%o13) log10(x) := ---------
log(10.0)
(%i14) omniout_str(iolevel, str) := if glob_iolevel >= iolevel
then printf(true, "~a~%", string(str))
(%o14) omniout_str(iolevel, str) := if glob_iolevel >= iolevel
then printf(true, "~a~%", string(str))
(%i15) omniout_str_noeol(iolevel, str) :=
if glob_iolevel >= iolevel then printf(true, "~a", string(str))
(%o15) omniout_str_noeol(iolevel, str) :=
if glob_iolevel >= iolevel then printf(true, "~a", string(str))
(%i16) omniout_labstr(iolevel, label, str) :=
if glob_iolevel >= iolevel then printf(true, "~a = ~a~%", string(label),
string(str))
(%o16) omniout_labstr(iolevel, label, str) :=
if glob_iolevel >= iolevel then printf(true, "~a = ~a~%", string(label),
string(str))
(%i17) omniout_float(iolevel, prelabel, prelen, value, vallen, postlabel) :=
if glob_iolevel >= iolevel then (if vallen = 4
then printf(true, "~a = ~g ~s ~%", prelabel, value, postlabel)
else printf(true, "~a = ~g ~s ~%", prelabel, value, postlabel))
(%o17) omniout_float(iolevel, prelabel, prelen, value, vallen, postlabel) :=
if glob_iolevel >= iolevel then (if vallen = 4
then printf(true, "~a = ~g ~s ~%", prelabel, value, postlabel)
else printf(true, "~a = ~g ~s ~%", prelabel, value, postlabel))
(%i18) omniout_int(iolevel, prelabel, prelen, value, vallen, postlabel) :=
if glob_iolevel >= iolevel then (printf(true, "~a = ~d ~a~%", prelabel, value,
postlabel), newline())
(%o18) omniout_int(iolevel, prelabel, prelen, value, vallen, postlabel) :=
if glob_iolevel >= iolevel then (printf(true, "~a = ~d ~a~%", prelabel, value,
postlabel), newline())
(%i19) omniout_float_arr(iolevel, prelabel, elemnt, prelen, value, vallen,
postlabel) := if glob_iolevel >= iolevel
then (sprint(prelabel, "[", elemnt, "]=", value, postlabel), newline())
(%o19) omniout_float_arr(iolevel, prelabel, elemnt, prelen, value, vallen,
postlabel) := if glob_iolevel >= iolevel
then (sprint(prelabel, "[", elemnt, "]=", value, postlabel), newline())
(%i20) dump_series(iolevel, dump_label, series_name, arr_series, numb) :=
block([i], if glob_iolevel >= iolevel then (i : 1,
while i <= numb do (sprint(dump_label, series_name, "i = ", i, "series = ",
array_series ), newline(), i : 1 + i)))
i
(%o20) dump_series(iolevel, dump_label, series_name, arr_series, numb) :=
block([i], if glob_iolevel >= iolevel then (i : 1,
while i <= numb do (sprint(dump_label, series_name, "i = ", i, "series = ",
array_series ), newline(), i : 1 + i)))
i
(%i21) dump_series_2(iolevel, dump_label, series_name2, arr_series2, numb,
subnum, arr_x) := (array_series2, numb, subnum) :=
block([i, sub, ts_term], if glob_iolevel >= iolevel
then (sub : 1, while sub <= subnum do (i : 1,
while i <= num do (sprint(dump_label, series_name, "sub = ", sub, "i = ", i,
"series2 = ", array_series2 ), i : 1 + i), sub : 1 + sub)))
sub, i
(%o21) dump_series_2(iolevel, dump_label, series_name2, arr_series2, numb,
subnum, arr_x) := (array_series2, numb, subnum) :=
block([i, sub, ts_term], if glob_iolevel >= iolevel
then (sub : 1, while sub <= subnum do (i : 1,
while i <= num do (sprint(dump_label, series_name, "sub = ", sub, "i = ", i,
"series2 = ", array_series2 ), i : 1 + i), sub : 1 + sub)))
sub, i
(%i22) cs_info(iolevel, str) := if glob_iolevel >= iolevel
then sprint(concat("cs_info ", str, " glob_correct_start_flag = ",
glob_correct_start_flag, "glob_h := ", glob_h, "glob_reached_optimal_h := ",
glob_reached_optimal_h))
(%o22) cs_info(iolevel, str) := if glob_iolevel >= iolevel
then sprint(concat("cs_info ", str, " glob_correct_start_flag = ",
glob_correct_start_flag, "glob_h := ", glob_h, "glob_reached_optimal_h := ",
glob_reached_optimal_h))
(%i23) logitem_time(fd, secs_in) := block([days, days_int, hours, hours_int,
minutes, minutes_int, sec_int, seconds, secs, years, years_int],
secs : convfloat(secs_in), printf(fd, "
~%"),
secs
if secs >= 0 then (years_int : trunc(----------------),
glob_sec_in_year
sec_temp : mod(trunc(secs), trunc(glob_sec_in_year)),
sec_temp
days_int : trunc(---------------), sec_temp :
glob_sec_in_day
sec_temp
mod(sec_temp, trunc(glob_sec_in_day)), hours_int : trunc(----------------),
glob_sec_in_hour
sec_temp : mod(sec_temp, trunc(glob_sec_in_hour)),
sec_temp
minutes_int : trunc(------------------),
glob_sec_in_minute
sec_int : mod(sec_temp, trunc(glob_sec_in_minute)),
if years_int > 0 then printf(fd,
"= ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", years_int, days_int,
hours_int, minutes_int, sec_int) elseif days_int > 0
then printf(fd, "= ~d Days ~d Hours ~d Minutes ~d Seconds~%", days_int,
hours_int, minutes_int, sec_int) elseif hours_int > 0
then printf(fd, "= ~d Hours ~d Minutes ~d Seconds~%", hours_int, minutes_int,
sec_int) elseif minutes_int > 0 then printf(fd, "= ~d Minutes ~d Seconds~%",
minutes_int, sec_int) else printf(fd, "= ~d Seconds~%", sec_int))
else printf(fd, " Unknown~%"), printf(fd, " | ~%"))
(%o23) logitem_time(fd, secs_in) := block([days, days_int, hours, hours_int,
minutes, minutes_int, sec_int, seconds, secs, years, years_int],
secs : convfloat(secs_in), printf(fd, "~%"),
secs
if secs >= 0 then (years_int : trunc(----------------),
glob_sec_in_year
sec_temp : mod(trunc(secs), trunc(glob_sec_in_year)),
sec_temp
days_int : trunc(---------------), sec_temp :
glob_sec_in_day
sec_temp
mod(sec_temp, trunc(glob_sec_in_day)), hours_int : trunc(----------------),
glob_sec_in_hour
sec_temp : mod(sec_temp, trunc(glob_sec_in_hour)),
sec_temp
minutes_int : trunc(------------------),
glob_sec_in_minute
sec_int : mod(sec_temp, trunc(glob_sec_in_minute)),
if years_int > 0 then printf(fd,
"= ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", years_int, days_int,
hours_int, minutes_int, sec_int) elseif days_int > 0
then printf(fd, "= ~d Days ~d Hours ~d Minutes ~d Seconds~%", days_int,
hours_int, minutes_int, sec_int) elseif hours_int > 0
then printf(fd, "= ~d Hours ~d Minutes ~d Seconds~%", hours_int, minutes_int,
sec_int) elseif minutes_int > 0 then printf(fd, "= ~d Minutes ~d Seconds~%",
minutes_int, sec_int) else printf(fd, "= ~d Seconds~%", sec_int))
else printf(fd, " Unknown~%"), printf(fd, " | ~%"))
(%i24) omniout_timestr(secs_in) := block([days, days_int, hours, hours_int,
minutes, minutes_int, sec_int, seconds, secs, years, years_int],
secs : convfloat(secs_in), if secs >= 0
secs
then (years_int : trunc(----------------),
glob_sec_in_year
sec_temp : mod(trunc(secs), trunc(glob_sec_in_year)),
sec_temp
days_int : trunc(---------------), sec_temp :
glob_sec_in_day
sec_temp
mod(sec_temp, trunc(glob_sec_in_day)), hours_int : trunc(----------------),
glob_sec_in_hour
sec_temp : mod(sec_temp, trunc(glob_sec_in_hour)),
sec_temp
minutes_int : trunc(------------------),
glob_sec_in_minute
sec_int : mod(sec_temp, trunc(glob_sec_in_minute)),
if years_int > 0 then printf(true,
"= ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", years_int, days_int,
hours_int, minutes_int, sec_int) elseif days_int > 0
then printf(true, "= ~d Days ~d Hours ~d Minutes ~d Seconds~%", days_int,
hours_int, minutes_int, sec_int) elseif hours_int > 0
then printf(true, "= ~d Hours ~d Minutes ~d Seconds~%", hours_int,
minutes_int, sec_int) elseif minutes_int > 0
then printf(true, "= ~d Minutes ~d Seconds~%", minutes_int, sec_int)
else printf(true, "= ~d Seconds~%", sec_int)) else printf(true, " Unknown~%"))
(%o24) omniout_timestr(secs_in) := block([days, days_int, hours, hours_int,
minutes, minutes_int, sec_int, seconds, secs, years, years_int],
secs : convfloat(secs_in), if secs >= 0
secs
then (years_int : trunc(----------------),
glob_sec_in_year
sec_temp : mod(trunc(secs), trunc(glob_sec_in_year)),
sec_temp
days_int : trunc(---------------), sec_temp :
glob_sec_in_day
sec_temp
mod(sec_temp, trunc(glob_sec_in_day)), hours_int : trunc(----------------),
glob_sec_in_hour
sec_temp : mod(sec_temp, trunc(glob_sec_in_hour)),
sec_temp
minutes_int : trunc(------------------),
glob_sec_in_minute
sec_int : mod(sec_temp, trunc(glob_sec_in_minute)),
if years_int > 0 then printf(true,
"= ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", years_int, days_int,
hours_int, minutes_int, sec_int) elseif days_int > 0
then printf(true, "= ~d Days ~d Hours ~d Minutes ~d Seconds~%", days_int,
hours_int, minutes_int, sec_int) elseif hours_int > 0
then printf(true, "= ~d Hours ~d Minutes ~d Seconds~%", hours_int,
minutes_int, sec_int) elseif minutes_int > 0
then printf(true, "= ~d Minutes ~d Seconds~%", minutes_int, sec_int)
else printf(true, "= ~d Seconds~%", sec_int)) else printf(true, " Unknown~%"))
(%i25) ats(mmm_ats, arr_a, arr_b, jjj_ats) :=
block([iii_ats, lll_ats, ma_ats, ret_ats], ret_ats : 0.0,
if jjj_ats <= mmm_ats then (ma_ats : 1 + mmm_ats, iii_ats : jjj_ats,
while iii_ats <= mmm_ats do (lll_ats : ma_ats - iii_ats,
ret_ats : arr_a arr_b + ret_ats, iii_ats : 1 + iii_ats)),
iii_ats lll_ats
ret_ats)
(%o25) ats(mmm_ats, arr_a, arr_b, jjj_ats) :=
block([iii_ats, lll_ats, ma_ats, ret_ats], ret_ats : 0.0,
if jjj_ats <= mmm_ats then (ma_ats : 1 + mmm_ats, iii_ats : jjj_ats,
while iii_ats <= mmm_ats do (lll_ats : ma_ats - iii_ats,
ret_ats : arr_a arr_b + ret_ats, iii_ats : 1 + iii_ats)),
iii_ats lll_ats
ret_ats)
(%i26) att(mmm_att, arr_aa, arr_bb, jjj_att) :=
block([al_att, iii_att, lll_att, ma_att, ret_att], ret_att : 0.0,
if jjj_att <= mmm_att then (ma_att : 2 + mmm_att, iii_att : jjj_att,
while iii_att <= mmm_att do (lll_att : ma_att - iii_att,
al_att : lll_att - 1, if lll_att <= glob_max_terms
then ret_att : arr_aa arr_bb convfp(al_att) + ret_att,
iii_att lll_att
ret_att
iii_att : 1 + iii_att), ret_att : ---------------), ret_att)
convfp(mmm_att)
(%o26) att(mmm_att, arr_aa, arr_bb, jjj_att) :=
block([al_att, iii_att, lll_att, ma_att, ret_att], ret_att : 0.0,
if jjj_att <= mmm_att then (ma_att : 2 + mmm_att, iii_att : jjj_att,
while iii_att <= mmm_att do (lll_att : ma_att - iii_att,
al_att : lll_att - 1, if lll_att <= glob_max_terms
then ret_att : arr_aa arr_bb convfp(al_att) + ret_att,
iii_att lll_att
ret_att
iii_att : 1 + iii_att), ret_att : ---------------), ret_att)
convfp(mmm_att)
(%i27) display_pole() := if (array_pole # glob_large_float)
1
and (array_pole > 0.0) and (array_pole # glob_large_float)
1 2
and (array_pole > 0.0) and glob_display_flag
2
then (omniout_float(ALWAYS, "Radius of convergence ", 4,
array_pole , 4, " "), omniout_float(ALWAYS,
1
"Order of pole ", 4, array_pole , 4, " "))
2
(%o27) display_pole() := if (array_pole # glob_large_float)
1
and (array_pole > 0.0) and (array_pole # glob_large_float)
1 2
and (array_pole > 0.0) and glob_display_flag
2
then (omniout_float(ALWAYS, "Radius of convergence ", 4,
array_pole , 4, " "), omniout_float(ALWAYS,
1
"Order of pole ", 4, array_pole , 4, " "))
2
(%i28) logditto(file) := (printf(file, ""), printf(file, "ditto"),
printf(file, " | "))
(%o28) logditto(file) := (printf(file, ""), printf(file, "ditto"),
printf(file, " | "))
(%i29) logitem_integer(file, n) := (printf(file, ""),
printf(file, "~d", n), printf(file, " | "))
(%o29) logitem_integer(file, n) := (printf(file, ""),
printf(file, "~d", n), printf(file, " | "))
(%i30) logitem_str(file, str) := (printf(file, ""), printf(file, str),
printf(file, " | "))
(%o30) logitem_str(file, str) := (printf(file, ""), printf(file, str),
printf(file, " | "))
(%i31) logitem_good_digits(file, rel_error) :=
block([good_digits], printf(file, ""),
if rel_error # - 1.0 then (if rel_error > + 1.0E-34
then (good_digits : 1 - floor(log10(rel_error)),
printf(file, "~d", good_digits)) else (good_digits : 16,
printf(file, "~d", good_digits))) else printf(file, "Unknown"),
printf(file, " | "))
(%o31) logitem_good_digits(file, rel_error) :=
block([good_digits], printf(file, ""),
if rel_error # - 1.0 then (if rel_error > + 1.0E-34
then (good_digits : 1 - floor(log10(rel_error)),
printf(file, "~d", good_digits)) else (good_digits : 16,
printf(file, "~d", good_digits))) else printf(file, "Unknown"),
printf(file, " | "))
(%i32) log_revs(file, revs) := printf(file, revs)
(%o32) log_revs(file, revs) := printf(file, revs)
(%i33) logitem_float(file, x) := (printf(file, ""), printf(file, "~g", x),
printf(file, " | "))
(%o33) logitem_float(file, x) := (printf(file, ""), printf(file, "~g", x),
printf(file, " | "))
(%i34) logitem_pole(file, pole) := (printf(file, ""),
if pole = 0 then printf(file, "NA") elseif pole = 1 then printf(file, "Real")
elseif pole = 2 then printf(file, "Complex") else printf(file, "No Pole"),
printf(file, " | "))
(%o34) logitem_pole(file, pole) := (printf(file, ""),
if pole = 0 then printf(file, "NA") elseif pole = 1 then printf(file, "Real")
elseif pole = 2 then printf(file, "Complex") else printf(file, "No Pole"),
printf(file, " | "))
(%i35) logstart(file) := printf(file, "")
(%o35) logstart(file) := printf(file, "
")
(%i36) logend(file) := printf(file, "
~%")
(%o36) logend(file) := printf(file, "~%")
(%i37) chk_data() := block([errflag], errflag : false,
if (glob_max_terms < 15) or (glob_max_terms > 512)
then (omniout_str(ALWAYS, "Illegal max_terms = -- Using 30"),
glob_max_terms : 30), if glob_max_iter < 2
then (omniout_str(ALWAYS, "Illegal max_iter"), errflag : true),
if errflag then quit())
(%o37) chk_data() := block([errflag], errflag : false,
if (glob_max_terms < 15) or (glob_max_terms > 512)
then (omniout_str(ALWAYS, "Illegal max_terms = -- Using 30"),
glob_max_terms : 30), if glob_max_iter < 2
then (omniout_str(ALWAYS, "Illegal max_iter"), errflag : true),
if errflag then quit())
(%i38) comp_expect_sec(t_end2, t_start2, t2, clock_sec2) :=
block([ms2, rrr, sec_left, sub1, sub2], ms2 : clock_sec2,
sub1 : t_end2 - t_start2, sub2 : t2 - t_start2,
if sub1 = 0.0 then sec_left : 0.0 else (if sub2 > 0.0
sub1
then (rrr : ----, sec_left : rrr ms2 - ms2) else sec_left : 0.0), sec_left)
sub2
(%o38) comp_expect_sec(t_end2, t_start2, t2, clock_sec2) :=
block([ms2, rrr, sec_left, sub1, sub2], ms2 : clock_sec2,
sub1 : t_end2 - t_start2, sub2 : t2 - t_start2,
if sub1 = 0.0 then sec_left : 0.0 else (if sub2 > 0.0
sub1
then (rrr : ----, sec_left : rrr ms2 - ms2) else sec_left : 0.0), sec_left)
sub2
(%i39) comp_percent(t_end2, t_start2, t2) :=
block([rrr, sub1, sub2], sub1 : t_end2 - t_start2, sub2 : t2 - t_start2,
100.0 sub2
if sub2 > glob_small_float then rrr : ---------- else rrr : 0.0, rrr)
sub1
(%o39) comp_percent(t_end2, t_start2, t2) :=
block([rrr, sub1, sub2], sub1 : t_end2 - t_start2, sub2 : t2 - t_start2,
100.0 sub2
if sub2 > glob_small_float then rrr : ---------- else rrr : 0.0, rrr)
sub1
(%i40) factorial_2(nnn) := nnn!
(%o40) factorial_2(nnn) := nnn!
(%i41) factorial_1(nnn) := block([ret],
if nnn <= glob_max_terms then (if array_fact_1 = 0
nnn
then (ret : factorial_2(nnn), array_fact_1 : ret)
nnn
else ret : array_fact_1 ) else ret : factorial_2(nnn), ret)
nnn
(%o41) factorial_1(nnn) := block([ret],
if nnn <= glob_max_terms then (if array_fact_1 = 0
nnn
then (ret : factorial_2(nnn), array_fact_1 : ret)
nnn
else ret : array_fact_1 ) else ret : factorial_2(nnn), ret)
nnn
(%i42) factorial_3(mmm, nnn) := block([ret],
if (nnn <= glob_max_terms) and (mmm <= glob_max_terms)
factorial_1(mmm)
then (if array_fact_2 = 0 then (ret : ----------------,
mmm, nnn factorial_1(nnn)
array_fact_2 : ret) else ret : array_fact_2 )
mmm, nnn mmm, nnn
factorial_2(mmm)
else ret : ----------------, ret)
factorial_2(nnn)
(%o42) factorial_3(mmm, nnn) := block([ret],
if (nnn <= glob_max_terms) and (mmm <= glob_max_terms)
factorial_1(mmm)
then (if array_fact_2 = 0 then (ret : ----------------,
mmm, nnn factorial_1(nnn)
array_fact_2 : ret) else ret : array_fact_2 )
mmm, nnn mmm, nnn
factorial_2(mmm)
else ret : ----------------, ret)
factorial_2(nnn)
(%i43) convfp(mmm) := mmm
(%o43) convfp(mmm) := mmm
(%i44) convfloat(mmm) := mmm
(%o44) convfloat(mmm) := mmm
(%i45) elapsed_time_seconds() := block([t], t : elapsed_real_time(), t)
(%o45) elapsed_time_seconds() := block([t], t : elapsed_real_time(), t)
(%i46) Si(x) := 0.0
(%o46) Si(x) := 0.0
(%i47) Ci(x) := 0.0
(%o47) Ci(x) := 0.0
(%i48) ln(x) := log(x)
(%o48) ln(x) := log(x)
(%i49) arcsin(x) := asin(x)
(%o49) arcsin(x) := asin(x)
(%i50) arccos(x) := acos(x)
(%o50) arccos(x) := acos(x)
(%i51) arctan(x) := atan(x)
(%o51) arctan(x) := atan(x)
(%i52) omniabs(x) := abs(x)
(%o52) omniabs(x) := abs(x)
(%i53) expt(x, y) := (if (x = 0.0) and (y < 0.0)
y
then print("expt error x = ", x, "y = ", y), x )
(%o53) expt(x, y) := (if (x = 0.0) and (y < 0.0)
y
then print("expt error x = ", x, "y = ", y), x )
(%i54) estimated_needed_step_error(x_start, x_end, estimated_h,
estimated_answer) := block([desired_abs_gbl_error, range, estimated_steps,
step_error], omniout_float(ALWAYS, "glob_desired_digits_correct", 32,
glob_desired_digits_correct, 32, ""), desired_abs_gbl_error :
expt(10.0, - glob_desired_digits_correct) omniabs(estimated_answer),
omniout_float(ALWAYS, "desired_abs_gbl_error", 32, desired_abs_gbl_error, 32,
""), range : x_end - x_start, omniout_float(ALWAYS, "range", 32, range, 32,
range
""), estimated_steps : -----------, omniout_float(ALWAYS, "estimated_steps",
estimated_h
desired_abs_gbl_error
32, estimated_steps, 32, ""), step_error : omniabs(---------------------),
estimated_steps
omniout_float(ALWAYS, "step_error", 32, step_error, 32, ""), step_error)
(%o54) estimated_needed_step_error(x_start, x_end, estimated_h,
estimated_answer) := block([desired_abs_gbl_error, range, estimated_steps,
step_error], omniout_float(ALWAYS, "glob_desired_digits_correct", 32,
glob_desired_digits_correct, 32, ""), desired_abs_gbl_error :
expt(10.0, - glob_desired_digits_correct) omniabs(estimated_answer),
omniout_float(ALWAYS, "desired_abs_gbl_error", 32, desired_abs_gbl_error, 32,
""), range : x_end - x_start, omniout_float(ALWAYS, "range", 32, range, 32,
range
""), estimated_steps : -----------, omniout_float(ALWAYS, "estimated_steps",
estimated_h
desired_abs_gbl_error
32, estimated_steps, 32, ""), step_error : omniabs(---------------------),
estimated_steps
omniout_float(ALWAYS, "step_error", 32, step_error, 32, ""), step_error)
ln(cosh(1.0 + 3.0 x))
(%i55) exact_soln_y(x) := block(---------------------)
3.0
ln(cosh(1.0 + 3.0 x))
(%o55) exact_soln_y(x) := block(---------------------)
3.0
(%i56) main() := block([d1, d2, d3, d4, est_err_2, niii, done_once, term, ord,
order_diff, term_no, html_log_file, iiif, jjjf, rows, r_order, sub_iter,
calc_term, iii, temp_sum, current_iter, x_start, x_end, it, log10norm,
max_terms, opt_iter, tmp, subiter, est_needed_step_err, value3, min_value,
est_answer, best_h, found_h, repeat_it],
define_variable(glob_max_terms, 30, fixnum),
define_variable(glob_iolevel, 5, fixnum), define_variable(ALWAYS, 1, fixnum),
define_variable(INFO, 2, fixnum), define_variable(DEBUGL, 3, fixnum),
define_variable(DEBUGMASSIVE, 4, fixnum),
define_variable(MAX_UNCHANGED, 10, fixnum),
define_variable(glob_check_sign, 1.0, float),
define_variable(glob_desired_digits_correct, 8.0, float),
define_variable(glob_max_value3, 0.0, float),
define_variable(glob_ratio_of_radius, 0.01, float),
define_variable(glob_percent_done, 0.0, float),
define_variable(glob_subiter_method, 3, fixnum),
define_variable(glob_log10normmin, 0.1, float),
define_variable(glob_total_exp_sec, 0.1, float),
define_variable(glob_optimal_expect_sec, 0.1, float),
define_variable(glob_html_log, true, boolean),
define_variable(glob_good_digits, 0, fixnum),
define_variable(glob_max_opt_iter, 10, fixnum),
define_variable(glob_dump, false, boolean),
define_variable(glob_djd_debug, true, boolean),
define_variable(glob_display_flag, true, boolean),
define_variable(glob_djd_debug2, true, boolean),
define_variable(glob_sec_in_minute, 60, fixnum),
define_variable(glob_min_in_hour, 60, fixnum),
define_variable(glob_hours_in_day, 24, fixnum),
define_variable(glob_days_in_year, 365, fixnum),
define_variable(glob_sec_in_hour, 3600, fixnum),
define_variable(glob_sec_in_day, 86400, fixnum),
define_variable(glob_sec_in_year, 31536000, fixnum),
define_variable(glob_almost_1, 0.999, float),
define_variable(glob_clock_sec, 0.0, float),
define_variable(glob_clock_start_sec, 0.0, float),
define_variable(glob_not_yet_finished, true, boolean),
define_variable(glob_initial_pass, true, boolean),
define_variable(glob_not_yet_start_msg, true, boolean),
define_variable(glob_reached_optimal_h, false, boolean),
define_variable(glob_optimal_done, false, boolean),
define_variable(glob_disp_incr, 0.1, float),
define_variable(glob_h, 0.1, float), define_variable(glob_hmax, 1.0, float),
define_variable(glob_hmin, 1.0E-11, float),
define_variable(glob_hmin_init, 0.001, float),
define_variable(glob_large_float, 9.0E+100, float),
define_variable(glob_last_good_h, 0.1, float),
define_variable(glob_look_poles, false, boolean),
define_variable(glob_neg_h, false, boolean),
define_variable(glob_display_interval, 0.0, float),
define_variable(glob_next_display, 0.0, float),
define_variable(glob_dump_analytic, false, boolean),
define_variable(glob_log10_abserr, 1.0E-11, float),
define_variable(glob_log10_relerr, 1.0E-11, float),
define_variable(glob_abserr, 1.0E-11, float),
define_variable(glob_relerr, 1.0E-11, float),
define_variable(glob_max_hours, 0.0, float),
define_variable(glob_max_iter, 1000, fixnum),
define_variable(glob_max_rel_trunc_err, 1.0E-11, float),
define_variable(glob_max_trunc_err, 1.0E-11, float),
define_variable(glob_no_eqs, 0, fixnum),
define_variable(glob_optimal_clock_start_sec, 0.0, float),
define_variable(glob_optimal_start, 0.0, float),
define_variable(glob_small_float, 1.0E-51, float),
define_variable(glob_smallish_float, 1.0E-101, float),
define_variable(glob_unchanged_h_cnt, 0, fixnum),
define_variable(glob_warned, false, boolean),
define_variable(glob_warned2, false, boolean),
define_variable(glob_max_sec, 10000.0, float),
define_variable(glob_orig_start_sec, 0.0, float),
define_variable(glob_start, 0, fixnum),
define_variable(glob_curr_iter_when_opt, 0, fixnum),
define_variable(glob_current_iter, 0, fixnum),
define_variable(glob_iter, 0, fixnum),
define_variable(glob_normmax, 0.0, float),
define_variable(glob_log10abserr, 0.0, float),
define_variable(glob_log10relerr, 0.0, float),
define_variable(glob_max_minutes, 0.0, float), ALWAYS : 1, INFO : 2,
DEBUGL : 3, DEBUGMASSIVE : 4, glob_iolevel : INFO,
glob_orig_start_sec : elapsed_time_seconds(), MAX_UNCHANGED : 10,
glob_curr_iter_when_opt : 0, glob_display_flag : true, glob_no_eqs : 1,
glob_iter : - 1, opt_iter : - 1, glob_max_iter : 50000, glob_max_hours : 0.0,
glob_max_minutes : 15.0, omniout_str(ALWAYS,
"##############ECHO OF PROBLEM#################"),
omniout_str(ALWAYS, "##############temp/lin_tanhpostode.ode#################"),
omniout_str(ALWAYS, "diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"),
omniout_str(ALWAYS, "!"), omniout_str(ALWAYS,
"/* BEGIN FIRST INPUT BLOCK */"), omniout_str(ALWAYS, "Digits:32,"),
omniout_str(ALWAYS, "max_terms:30,"), omniout_str(ALWAYS, "!"),
omniout_str(ALWAYS, "/* END FIRST INPUT BLOCK */"),
omniout_str(ALWAYS, "/* BEGIN SECOND INPUT BLOCK */"),
omniout_str(ALWAYS, "x_start:1.1,"), omniout_str(ALWAYS, "x_end:2.0,"),
omniout_str(ALWAYS, "array_y_init[0 + 1] : exact_soln_y(x_start),"),
omniout_str(ALWAYS, "glob_h:0.00001,"),
omniout_str(ALWAYS, "glob_look_poles:true,"),
omniout_str(ALWAYS, "glob_max_iter:10,"),
omniout_str(ALWAYS, "/* END SECOND INPUT BLOCK */"),
omniout_str(ALWAYS, "/* BEGIN OVERRIDE BLOCK */"),
omniout_str(ALWAYS, "glob_desired_digits_correct:10,"),
omniout_str(ALWAYS, "glob_display_interval:0.001,"),
omniout_str(ALWAYS, "glob_look_poles:true,"),
omniout_str(ALWAYS, "glob_max_iter:10000000,"),
omniout_str(ALWAYS, "glob_max_minutes:3,"),
omniout_str(ALWAYS, "glob_subiter_method:3,"),
omniout_str(ALWAYS, "/* END OVERRIDE BLOCK */"), omniout_str(ALWAYS, "!"),
omniout_str(ALWAYS, "/* BEGIN USER DEF BLOCK */"),
omniout_str(ALWAYS, "exact_soln_y (x) := (block("),
omniout_str(ALWAYS, " (ln(cosh(3.0*x + 1.0))/3.0) "),
omniout_str(ALWAYS, "));"), omniout_str(ALWAYS, ""),
omniout_str(ALWAYS, "/* END USER DEF BLOCK */"),
omniout_str(ALWAYS, "#######END OF ECHO OF PROBLEM#################"),
glob_unchanged_h_cnt : 0, glob_warned : false, glob_warned2 : false,
glob_small_float : 1.0E-200, glob_smallish_float : 1.0E-64,
glob_large_float : 1.0E+100, glob_almost_1 : 0.99, glob_log10_abserr : - 8.0,
glob_log10_relerr : - 8.0, glob_hmax : 0.01, Digits : 32, max_terms : 30,
glob_max_terms : max_terms, glob_html_log : true,
array(array_y_init, 1 + max_terms), array(array_norms, 1 + max_terms),
array(array_fact_1, 1 + max_terms), array(array_pole, 1 + max_terms),
array(array_1st_rel_error, 1 + max_terms),
array(array_last_rel_error, 1 + max_terms),
array(array_type_pole, 1 + max_terms), array(array_y, 1 + max_terms),
array(array_x, 1 + max_terms), array(array_tmp0, 1 + max_terms),
array(array_tmp1, 1 + max_terms), array(array_tmp2, 1 + max_terms),
array(array_tmp3_g, 1 + max_terms), array(array_tmp3_a1, 1 + max_terms),
array(array_tmp3_a2, 1 + max_terms), array(array_tmp3, 1 + max_terms),
array(array_tmp4, 1 + max_terms), array(array_m1, 1 + max_terms),
array(array_y_higher, 1 + 2, 1 + max_terms),
array(array_y_higher_work, 1 + 2, 1 + max_terms),
array(array_y_higher_work2, 1 + 2, 1 + max_terms),
array(array_y_set_initial, 1 + 2, 1 + max_terms),
array(array_poles, 1 + 1, 1 + 3), array(array_real_pole, 1 + 1, 1 + 3),
array(array_complex_pole, 1 + 1, 1 + 3),
array(array_fact_2, 1 + max_terms, 1 + max_terms), term : 1,
while term <= max_terms do (array_y_init : 0.0, term : 1 + term),
term
term : 1, while term <= max_terms do (array_norms : 0.0,
term
term : 1 + term), term : 1, while term <=
max_terms do (array_fact_1 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_pole : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_1st_rel_error : 0.0, term : 1 + term),
term
term : 1, while term <= max_terms do (array_last_rel_error : 0.0,
term
term : 1 + term), term : 1, while term <=
max_terms do (array_type_pole : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_y : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_x : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp0 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp1 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp2 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp3_g : 0.0, term : 1 + term),
term
term : 1, while term <= max_terms do (array_tmp3_a1 : 0.0,
term
term : 1 + term), term : 1, while term <=
max_terms do (array_tmp3_a2 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp3 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp4 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_m1 : 0.0, term : 1 + term), ord : 1,
term
while ord <= 2 do (term : 1, while term <=
max_terms do (array_y_higher : 0.0, term : 1 + term), ord : 1 + ord),
ord, term
ord : 1, while ord <= 2 do (term : 1,
while term <= max_terms do (array_y_higher_work : 0.0,
ord, term
term : 1 + term), ord : 1 + ord), ord : 1,
while ord <= 2 do (term : 1, while term <=
max_terms do (array_y_higher_work2 : 0.0, term : 1 + term),
ord, term
ord : 1 + ord), ord : 1, while ord <= 2 do (term : 1,
while term <= max_terms do (array_y_set_initial : 0.0,
ord, term
term : 1 + term), ord : 1 + ord), ord : 1,
while ord <= 1 do (term : 1, while term <= 3 do (array_poles : 0.0,
ord, term
term : 1 + term), ord : 1 + ord), ord : 1,
while ord <= 1 do (term : 1, while term <=
3 do (array_real_pole : 0.0, term : 1 + term), ord : 1 + ord),
ord, term
ord : 1, while ord <= 1 do (term : 1,
while term <= 3 do (array_complex_pole : 0.0, term : 1 + term),
ord, term
ord : 1 + ord), ord : 1, while ord <= max_terms do (term : 1,
while term <= max_terms do (array_fact_2 : 0.0, term : 1 + term),
ord, term
ord : 1 + ord), array(array_y, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_y : 0.0, term : 1 + term),
term
array(array_x, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_x : 0.0, term : 1 + term),
term
array(array_tmp0, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp0 : 0.0, term : 1 + term),
term
array(array_tmp1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp1 : 0.0, term : 1 + term),
term
array(array_tmp2, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp2 : 0.0, term : 1 + term),
term
array(array_tmp3_g, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3_g : 0.0, term : 1 + term),
term
array(array_tmp3_a1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3_a1 : 0.0, term : 1 + term),
term
array(array_tmp3_a2, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3_a2 : 0.0, term : 1 + term),
term
array(array_tmp3, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3 : 0.0, term : 1 + term),
term
array(array_tmp4, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp4 : 0.0, term : 1 + term),
term
array(array_m1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_m1 : 0.0, term : 1 + term),
term
array(array_const_1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_const_1 : 0.0, term : 1 + term),
term
array_const_1 : 1, array(array_const_0D0, 1 + 1 + max_terms), term : 1,
1
while term <= 1 + max_terms do (array_const_0D0 : 0.0, term : 1 + term),
term
array_const_0D0 : 0.0, array(array_const_3D0, 1 + 1 + max_terms), term : 1,
1
while term <= 1 + max_terms do (array_const_3D0 : 0.0, term : 1 + term),
term
array_const_3D0 : 3.0, array(array_const_1D0, 1 + 1 + max_terms), term : 1,
1
while term <= 1 + max_terms do (array_const_1D0 : 0.0, term : 1 + term),
term
array_const_1D0 : 1.0, array(array_m1, 1 + 1 + max_terms), term : 1,
1
while term <= max_terms do (array_m1 : 0.0, term : 1 + term),
term
array_m1 : - 1.0, iiif : 0, while iiif <= glob_max_terms do (jjjf : 0,
1
while jjjf <= glob_max_terms do (array_fact_1 : 0,
iiif
array_fact_2 : 0, jjjf : 1 + jjjf), iiif : 1 + iiif), x_start : 1.1,
iiif, jjjf
x_end : 2.0, array_y_init : exact_soln_y(x_start), glob_h : 1.0E-5,
1 + 0
glob_look_poles : true, glob_max_iter : 10, 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,
glob_last_good_h : glob_h, glob_max_terms : max_terms,
glob_max_sec : convfloat(3600.0) convfloat(glob_max_hours)
+ convfloat(60.0) convfloat(glob_max_minutes),
glob_abserr : expt(10.0, glob_log10_abserr),
glob_relerr : expt(10.0, glob_log10_relerr),
if glob_h > 0.0 then (glob_neg_h : false,
glob_display_interval : omniabs(glob_display_interval))
else (glob_neg_h : true, glob_display_interval :
- omniabs(glob_display_interval)), chk_data(), array_y_set_initial : true,
1, 1
array_y_set_initial : false, array_y_set_initial : false,
1, 2 1, 3
array_y_set_initial : false, array_y_set_initial : false,
1, 4 1, 5
array_y_set_initial : false, array_y_set_initial : false,
1, 6 1, 7
array_y_set_initial : false, array_y_set_initial : false,
1, 8 1, 9
array_y_set_initial : false, array_y_set_initial : false,
1, 10 1, 11
array_y_set_initial : false, array_y_set_initial : false,
1, 12 1, 13
array_y_set_initial : false, array_y_set_initial : false,
1, 14 1, 15
array_y_set_initial : false, array_y_set_initial : false,
1, 16 1, 17
array_y_set_initial : false, array_y_set_initial : false,
1, 18 1, 19
array_y_set_initial : false, array_y_set_initial : false,
1, 20 1, 21
array_y_set_initial : false, array_y_set_initial : false,
1, 22 1, 23
array_y_set_initial : false, array_y_set_initial : false,
1, 24 1, 25
array_y_set_initial : false, array_y_set_initial : false,
1, 26 1, 27
array_y_set_initial : false, array_y_set_initial : false,
1, 28 1, 29
array_y_set_initial : false, omniout_str(ALWAYS, "START of Optimize"),
1, 30
glob_check_sign : check_sign(x_start, x_end),
glob_h : check_sign(x_start, x_end), if glob_display_interval < glob_h
then glob_h : glob_display_interval, found_h : - 1.0, best_h : 0.0,
min_value : glob_large_float, est_answer : est_size_answer(), opt_iter : 1,
while (opt_iter <= 20) and (found_h < 0.0) do (omniout_int(ALWAYS,
"opt_iter", 32, opt_iter, 4, ""), array_x : x_start, array_x : glob_h,
1 2
glob_next_display : x_start, order_diff : 1, term_no : 1,
while term_no <= order_diff do (array_y :
term_no
array_y_init expt(glob_h, term_no - 1)
term_no
---------------------------------------------, term_no : 1 + term_no),
factorial_1(term_no - 1)
rows : order_diff, r_order : 1, while r_order <= rows do (term_no : 1,
while term_no <= 1 - r_order + rows do (it : - 1 + r_order + term_no,
array_y_init expt(glob_h, term_no - 1)
it
array_y_higher : ----------------------------------------,
r_order, term_no factorial_1(term_no - 1)
term_no : 1 + term_no), r_order : 1 + r_order), atomall(),
est_needed_step_err : estimated_needed_step_error(x_start, x_end, glob_h,
est_answer), omniout_float(ALWAYS, "est_needed_step_err", 32,
est_needed_step_err, 16, ""), value3 : test_suggested_h(),
omniout_float(ALWAYS, "value3", 32, value3, 32, ""),
if (value3 < est_needed_step_err) and (found_h < 0.0)
then (best_h : glob_h, found_h : 1.0),
omniout_float(ALWAYS, "best_h", 32, best_h, 32, ""), opt_iter : 1 + opt_iter,
glob_h : glob_h 0.5), if found_h > 0.0 then glob_h : best_h
else omniout_str(ALWAYS, "No increment to obtain desired accuracy found"),
if glob_html_log then html_log_file : openw("html/entry.html"),
if found_h > 0.0 then (omniout_str(ALWAYS, "START of Soultion"),
array_x : x_start, array_x : glob_h, glob_next_display : x_start,
1 2
order_diff : 1, term_no : 1, while term_no <=
order_diff do (array_y : (array_y_init expt(glob_h, term_no - 1))
term_no term_no
/factorial_1(term_no - 1), term_no : 1 + term_no), rows : order_diff,
r_order : 1, while r_order <= rows do (term_no : 1,
while term_no <= 1 - r_order + rows do (it : - 1 + r_order + term_no,
array_y_init expt(glob_h, term_no - 1)
it
array_y_higher : ----------------------------------------,
r_order, term_no factorial_1(term_no - 1)
term_no : 1 + term_no), r_order : 1 + r_order), current_iter : 1,
glob_clock_start_sec : elapsed_time_seconds(),
glob_log10normmin : - glob_large_float,
if omniabs(array_y_higher ) > glob_small_float
1, 1
then (tmp : omniabs(array_y_higher ), log10norm : log10(tmp),
1, 1
if log10norm < glob_log10normmin then glob_log10normmin : log10norm),
display_alot(current_iter), glob_clock_sec : elapsed_time_seconds(),
glob_current_iter : 0, glob_iter : 0, omniout_str(DEBUGL, " "),
glob_reached_optimal_h : true, glob_optimal_clock_start_sec :
elapsed_time_seconds(), while (glob_current_iter < glob_max_iter)
and (glob_check_sign array_x < glob_check_sign x_end)
1
and (convfloat(glob_clock_sec) - convfloat(glob_orig_start_sec) < convfloat(glob_max_sec)) do (if reached_interval
() then (omniout_str(INFO, " "), omniout_str(INFO, "TOP MAIN SOLVE Loop")),
glob_iter : 1 + glob_iter, glob_clock_sec : elapsed_time_seconds(),
glob_current_iter : 1 + glob_current_iter, atomall(),
display_alot(current_iter), if glob_look_poles then check_for_pole(),
if reached_interval() then glob_next_display :
glob_display_interval + glob_next_display, array_x : glob_h + array_x ,
1 1
array_x : glob_h, order_diff : 2, ord : 2, calc_term : 1,
2
iii : glob_max_terms, while iii >= calc_term do (array_y_higher_work :
2, iii
array_y_higher
2, iii
---------------------------
expt(glob_h, calc_term - 1)
-------------------------------------, iii : iii - 1), temp_sum : 0.0,
factorial_3(iii - calc_term, iii - 1)
ord : 2, calc_term : 1, iii : glob_max_terms,
while iii >= calc_term do (temp_sum : array_y_higher_work + temp_sum,
ord, iii
iii : iii - 1), array_y_higher_work2 :
ord, calc_term
temp_sum expt(glob_h, calc_term - 1)
------------------------------------, ord : 1, calc_term : 2,
factorial_1(calc_term - 1)
iii : glob_max_terms, while iii >= calc_term do (array_y_higher_work :
1, iii
array_y_higher
1, iii
---------------------------
expt(glob_h, calc_term - 1)
-------------------------------------, iii : iii - 1), temp_sum : 0.0,
factorial_3(iii - calc_term, iii - 1)
ord : 1, calc_term : 2, iii : glob_max_terms,
while iii >= calc_term do (temp_sum : array_y_higher_work + temp_sum,
ord, iii
iii : iii - 1), array_y_higher_work2 :
ord, calc_term
temp_sum expt(glob_h, calc_term - 1)
------------------------------------, ord : 1, calc_term : 1,
factorial_1(calc_term - 1)
iii : glob_max_terms, while iii >= calc_term do (array_y_higher_work :
1, iii
array_y_higher
1, iii
---------------------------
expt(glob_h, calc_term - 1)
-------------------------------------, iii : iii - 1), temp_sum : 0.0,
factorial_3(iii - calc_term, iii - 1)
ord : 1, calc_term : 1, iii : glob_max_terms,
while iii >= calc_term do (temp_sum : array_y_higher_work + temp_sum,
ord, iii
iii : iii - 1), array_y_higher_work2 :
ord, calc_term
temp_sum expt(glob_h, calc_term - 1)
------------------------------------, term_no : glob_max_terms,
factorial_1(calc_term - 1)
while term_no >= 1 do (array_y : array_y_higher_work2 ,
term_no 1, term_no
ord : 1, while ord <= order_diff do (array_y_higher :
ord, term_no
array_y_higher_work2 , ord : 1 + ord), term_no : term_no - 1)),
ord, term_no
omniout_str(ALWAYS, "Finished!"), if glob_iter >= glob_max_iter
then omniout_str(ALWAYS,
"Maximum Iterations Reached before Solution Completed!"),
if elapsed_time_seconds() - convfloat(glob_orig_start_sec) >=
convfloat(glob_max_sec) then omniout_str(ALWAYS,
"Maximum Time Reached before Solution Completed!"),
glob_clock_sec : elapsed_time_seconds(),
omniout_str(INFO, "diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"),
omniout_int(INFO, "Iterations ", 32, glob_iter, 4, " "),
prog_report(x_start, x_end), if glob_html_log
then (logstart(html_log_file), logitem_str(html_log_file,
"2013-01-13T00:13:13-06:00"), logitem_str(html_log_file, "Maxima"),
logitem_str(html_log_file, "lin_tanh"),
logitem_str(html_log_file, "diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"),
logitem_float(html_log_file, x_start), logitem_float(html_log_file, x_end),
logitem_float(html_log_file, array_x ), logitem_float(html_log_file, glob_h),
1
logitem_str(html_log_file, "16"), logitem_good_digits(html_log_file,
array_last_rel_error ), logitem_integer(html_log_file, glob_max_terms),
1
logitem_float(html_log_file, array_1st_rel_error ),
1
logitem_float(html_log_file, array_last_rel_error ),
1
logitem_integer(html_log_file, glob_iter),
logitem_pole(html_log_file, array_type_pole ),
1
if (array_type_pole = 1) or (array_type_pole = 2)
1 1
then (logitem_float(html_log_file, array_pole ),
1
logitem_float(html_log_file, array_pole ), 0)
2
else (logitem_str(html_log_file, "NA"), logitem_str(html_log_file, "NA"), 0),
logitem_time(html_log_file, convfloat(glob_clock_sec)),
if glob_percent_done < 100.0 then (logitem_time(html_log_file,
convfloat(glob_total_exp_sec)), 0) else (logitem_str(html_log_file, "Done"),
0), log_revs(html_log_file, " 156 | "), logitem_str(html_log_file, "lin_tanh diffeq.max"),
logitem_str(html_log_file,
"lin_tanh maxima results"),
logitem_str(html_log_file, "Languages compared - single equations"),
logend(html_log_file)), if glob_html_log then close(html_log_file)))
(%o56) main() := block([d1, d2, d3, d4, est_err_2, niii, done_once, term, ord,
order_diff, term_no, html_log_file, iiif, jjjf, rows, r_order, sub_iter,
calc_term, iii, temp_sum, current_iter, x_start, x_end, it, log10norm,
max_terms, opt_iter, tmp, subiter, est_needed_step_err, value3, min_value,
est_answer, best_h, found_h, repeat_it],
define_variable(glob_max_terms, 30, fixnum),
define_variable(glob_iolevel, 5, fixnum), define_variable(ALWAYS, 1, fixnum),
define_variable(INFO, 2, fixnum), define_variable(DEBUGL, 3, fixnum),
define_variable(DEBUGMASSIVE, 4, fixnum),
define_variable(MAX_UNCHANGED, 10, fixnum),
define_variable(glob_check_sign, 1.0, float),
define_variable(glob_desired_digits_correct, 8.0, float),
define_variable(glob_max_value3, 0.0, float),
define_variable(glob_ratio_of_radius, 0.01, float),
define_variable(glob_percent_done, 0.0, float),
define_variable(glob_subiter_method, 3, fixnum),
define_variable(glob_log10normmin, 0.1, float),
define_variable(glob_total_exp_sec, 0.1, float),
define_variable(glob_optimal_expect_sec, 0.1, float),
define_variable(glob_html_log, true, boolean),
define_variable(glob_good_digits, 0, fixnum),
define_variable(glob_max_opt_iter, 10, fixnum),
define_variable(glob_dump, false, boolean),
define_variable(glob_djd_debug, true, boolean),
define_variable(glob_display_flag, true, boolean),
define_variable(glob_djd_debug2, true, boolean),
define_variable(glob_sec_in_minute, 60, fixnum),
define_variable(glob_min_in_hour, 60, fixnum),
define_variable(glob_hours_in_day, 24, fixnum),
define_variable(glob_days_in_year, 365, fixnum),
define_variable(glob_sec_in_hour, 3600, fixnum),
define_variable(glob_sec_in_day, 86400, fixnum),
define_variable(glob_sec_in_year, 31536000, fixnum),
define_variable(glob_almost_1, 0.999, float),
define_variable(glob_clock_sec, 0.0, float),
define_variable(glob_clock_start_sec, 0.0, float),
define_variable(glob_not_yet_finished, true, boolean),
define_variable(glob_initial_pass, true, boolean),
define_variable(glob_not_yet_start_msg, true, boolean),
define_variable(glob_reached_optimal_h, false, boolean),
define_variable(glob_optimal_done, false, boolean),
define_variable(glob_disp_incr, 0.1, float),
define_variable(glob_h, 0.1, float), define_variable(glob_hmax, 1.0, float),
define_variable(glob_hmin, 1.0E-11, float),
define_variable(glob_hmin_init, 0.001, float),
define_variable(glob_large_float, 9.0E+100, float),
define_variable(glob_last_good_h, 0.1, float),
define_variable(glob_look_poles, false, boolean),
define_variable(glob_neg_h, false, boolean),
define_variable(glob_display_interval, 0.0, float),
define_variable(glob_next_display, 0.0, float),
define_variable(glob_dump_analytic, false, boolean),
define_variable(glob_log10_abserr, 1.0E-11, float),
define_variable(glob_log10_relerr, 1.0E-11, float),
define_variable(glob_abserr, 1.0E-11, float),
define_variable(glob_relerr, 1.0E-11, float),
define_variable(glob_max_hours, 0.0, float),
define_variable(glob_max_iter, 1000, fixnum),
define_variable(glob_max_rel_trunc_err, 1.0E-11, float),
define_variable(glob_max_trunc_err, 1.0E-11, float),
define_variable(glob_no_eqs, 0, fixnum),
define_variable(glob_optimal_clock_start_sec, 0.0, float),
define_variable(glob_optimal_start, 0.0, float),
define_variable(glob_small_float, 1.0E-51, float),
define_variable(glob_smallish_float, 1.0E-101, float),
define_variable(glob_unchanged_h_cnt, 0, fixnum),
define_variable(glob_warned, false, boolean),
define_variable(glob_warned2, false, boolean),
define_variable(glob_max_sec, 10000.0, float),
define_variable(glob_orig_start_sec, 0.0, float),
define_variable(glob_start, 0, fixnum),
define_variable(glob_curr_iter_when_opt, 0, fixnum),
define_variable(glob_current_iter, 0, fixnum),
define_variable(glob_iter, 0, fixnum),
define_variable(glob_normmax, 0.0, float),
define_variable(glob_log10abserr, 0.0, float),
define_variable(glob_log10relerr, 0.0, float),
define_variable(glob_max_minutes, 0.0, float), ALWAYS : 1, INFO : 2,
DEBUGL : 3, DEBUGMASSIVE : 4, glob_iolevel : INFO,
glob_orig_start_sec : elapsed_time_seconds(), MAX_UNCHANGED : 10,
glob_curr_iter_when_opt : 0, glob_display_flag : true, glob_no_eqs : 1,
glob_iter : - 1, opt_iter : - 1, glob_max_iter : 50000, glob_max_hours : 0.0,
glob_max_minutes : 15.0, omniout_str(ALWAYS,
"##############ECHO OF PROBLEM#################"),
omniout_str(ALWAYS, "##############temp/lin_tanhpostode.ode#################"),
omniout_str(ALWAYS, "diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"),
omniout_str(ALWAYS, "!"), omniout_str(ALWAYS,
"/* BEGIN FIRST INPUT BLOCK */"), omniout_str(ALWAYS, "Digits:32,"),
omniout_str(ALWAYS, "max_terms:30,"), omniout_str(ALWAYS, "!"),
omniout_str(ALWAYS, "/* END FIRST INPUT BLOCK */"),
omniout_str(ALWAYS, "/* BEGIN SECOND INPUT BLOCK */"),
omniout_str(ALWAYS, "x_start:1.1,"), omniout_str(ALWAYS, "x_end:2.0,"),
omniout_str(ALWAYS, "array_y_init[0 + 1] : exact_soln_y(x_start),"),
omniout_str(ALWAYS, "glob_h:0.00001,"),
omniout_str(ALWAYS, "glob_look_poles:true,"),
omniout_str(ALWAYS, "glob_max_iter:10,"),
omniout_str(ALWAYS, "/* END SECOND INPUT BLOCK */"),
omniout_str(ALWAYS, "/* BEGIN OVERRIDE BLOCK */"),
omniout_str(ALWAYS, "glob_desired_digits_correct:10,"),
omniout_str(ALWAYS, "glob_display_interval:0.001,"),
omniout_str(ALWAYS, "glob_look_poles:true,"),
omniout_str(ALWAYS, "glob_max_iter:10000000,"),
omniout_str(ALWAYS, "glob_max_minutes:3,"),
omniout_str(ALWAYS, "glob_subiter_method:3,"),
omniout_str(ALWAYS, "/* END OVERRIDE BLOCK */"), omniout_str(ALWAYS, "!"),
omniout_str(ALWAYS, "/* BEGIN USER DEF BLOCK */"),
omniout_str(ALWAYS, "exact_soln_y (x) := (block("),
omniout_str(ALWAYS, " (ln(cosh(3.0*x + 1.0))/3.0) "),
omniout_str(ALWAYS, "));"), omniout_str(ALWAYS, ""),
omniout_str(ALWAYS, "/* END USER DEF BLOCK */"),
omniout_str(ALWAYS, "#######END OF ECHO OF PROBLEM#################"),
glob_unchanged_h_cnt : 0, glob_warned : false, glob_warned2 : false,
glob_small_float : 1.0E-200, glob_smallish_float : 1.0E-64,
glob_large_float : 1.0E+100, glob_almost_1 : 0.99, glob_log10_abserr : - 8.0,
glob_log10_relerr : - 8.0, glob_hmax : 0.01, Digits : 32, max_terms : 30,
glob_max_terms : max_terms, glob_html_log : true,
array(array_y_init, 1 + max_terms), array(array_norms, 1 + max_terms),
array(array_fact_1, 1 + max_terms), array(array_pole, 1 + max_terms),
array(array_1st_rel_error, 1 + max_terms),
array(array_last_rel_error, 1 + max_terms),
array(array_type_pole, 1 + max_terms), array(array_y, 1 + max_terms),
array(array_x, 1 + max_terms), array(array_tmp0, 1 + max_terms),
array(array_tmp1, 1 + max_terms), array(array_tmp2, 1 + max_terms),
array(array_tmp3_g, 1 + max_terms), array(array_tmp3_a1, 1 + max_terms),
array(array_tmp3_a2, 1 + max_terms), array(array_tmp3, 1 + max_terms),
array(array_tmp4, 1 + max_terms), array(array_m1, 1 + max_terms),
array(array_y_higher, 1 + 2, 1 + max_terms),
array(array_y_higher_work, 1 + 2, 1 + max_terms),
array(array_y_higher_work2, 1 + 2, 1 + max_terms),
array(array_y_set_initial, 1 + 2, 1 + max_terms),
array(array_poles, 1 + 1, 1 + 3), array(array_real_pole, 1 + 1, 1 + 3),
array(array_complex_pole, 1 + 1, 1 + 3),
array(array_fact_2, 1 + max_terms, 1 + max_terms), term : 1,
while term <= max_terms do (array_y_init : 0.0, term : 1 + term),
term
term : 1, while term <= max_terms do (array_norms : 0.0,
term
term : 1 + term), term : 1, while term <=
max_terms do (array_fact_1 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_pole : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_1st_rel_error : 0.0, term : 1 + term),
term
term : 1, while term <= max_terms do (array_last_rel_error : 0.0,
term
term : 1 + term), term : 1, while term <=
max_terms do (array_type_pole : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_y : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_x : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp0 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp1 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp2 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp3_g : 0.0, term : 1 + term),
term
term : 1, while term <= max_terms do (array_tmp3_a1 : 0.0,
term
term : 1 + term), term : 1, while term <=
max_terms do (array_tmp3_a2 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp3 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_tmp4 : 0.0, term : 1 + term), term : 1,
term
while term <= max_terms do (array_m1 : 0.0, term : 1 + term), ord : 1,
term
while ord <= 2 do (term : 1, while term <=
max_terms do (array_y_higher : 0.0, term : 1 + term), ord : 1 + ord),
ord, term
ord : 1, while ord <= 2 do (term : 1,
while term <= max_terms do (array_y_higher_work : 0.0,
ord, term
term : 1 + term), ord : 1 + ord), ord : 1,
while ord <= 2 do (term : 1, while term <=
max_terms do (array_y_higher_work2 : 0.0, term : 1 + term),
ord, term
ord : 1 + ord), ord : 1, while ord <= 2 do (term : 1,
while term <= max_terms do (array_y_set_initial : 0.0,
ord, term
term : 1 + term), ord : 1 + ord), ord : 1,
while ord <= 1 do (term : 1, while term <= 3 do (array_poles : 0.0,
ord, term
term : 1 + term), ord : 1 + ord), ord : 1,
while ord <= 1 do (term : 1, while term <=
3 do (array_real_pole : 0.0, term : 1 + term), ord : 1 + ord),
ord, term
ord : 1, while ord <= 1 do (term : 1,
while term <= 3 do (array_complex_pole : 0.0, term : 1 + term),
ord, term
ord : 1 + ord), ord : 1, while ord <= max_terms do (term : 1,
while term <= max_terms do (array_fact_2 : 0.0, term : 1 + term),
ord, term
ord : 1 + ord), array(array_y, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_y : 0.0, term : 1 + term),
term
array(array_x, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_x : 0.0, term : 1 + term),
term
array(array_tmp0, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp0 : 0.0, term : 1 + term),
term
array(array_tmp1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp1 : 0.0, term : 1 + term),
term
array(array_tmp2, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp2 : 0.0, term : 1 + term),
term
array(array_tmp3_g, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3_g : 0.0, term : 1 + term),
term
array(array_tmp3_a1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3_a1 : 0.0, term : 1 + term),
term
array(array_tmp3_a2, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3_a2 : 0.0, term : 1 + term),
term
array(array_tmp3, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp3 : 0.0, term : 1 + term),
term
array(array_tmp4, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_tmp4 : 0.0, term : 1 + term),
term
array(array_m1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_m1 : 0.0, term : 1 + term),
term
array(array_const_1, 1 + 1 + max_terms), term : 1,
while term <= 1 + max_terms do (array_const_1 : 0.0, term : 1 + term),
term
array_const_1 : 1, array(array_const_0D0, 1 + 1 + max_terms), term : 1,
1
while term <= 1 + max_terms do (array_const_0D0 : 0.0, term : 1 + term),
term
array_const_0D0 : 0.0, array(array_const_3D0, 1 + 1 + max_terms), term : 1,
1
while term <= 1 + max_terms do (array_const_3D0 : 0.0, term : 1 + term),
term
array_const_3D0 : 3.0, array(array_const_1D0, 1 + 1 + max_terms), term : 1,
1
while term <= 1 + max_terms do (array_const_1D0 : 0.0, term : 1 + term),
term
array_const_1D0 : 1.0, array(array_m1, 1 + 1 + max_terms), term : 1,
1
while term <= max_terms do (array_m1 : 0.0, term : 1 + term),
term
array_m1 : - 1.0, iiif : 0, while iiif <= glob_max_terms do (jjjf : 0,
1
while jjjf <= glob_max_terms do (array_fact_1 : 0,
iiif
array_fact_2 : 0, jjjf : 1 + jjjf), iiif : 1 + iiif), x_start : 1.1,
iiif, jjjf
x_end : 2.0, array_y_init : exact_soln_y(x_start), glob_h : 1.0E-5,
1 + 0
glob_look_poles : true, glob_max_iter : 10, 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,
glob_last_good_h : glob_h, glob_max_terms : max_terms,
glob_max_sec : convfloat(3600.0) convfloat(glob_max_hours)
+ convfloat(60.0) convfloat(glob_max_minutes),
glob_abserr : expt(10.0, glob_log10_abserr),
glob_relerr : expt(10.0, glob_log10_relerr),
if glob_h > 0.0 then (glob_neg_h : false,
glob_display_interval : omniabs(glob_display_interval))
else (glob_neg_h : true, glob_display_interval :
- omniabs(glob_display_interval)), chk_data(), array_y_set_initial : true,
1, 1
array_y_set_initial : false, array_y_set_initial : false,
1, 2 1, 3
array_y_set_initial : false, array_y_set_initial : false,
1, 4 1, 5
array_y_set_initial : false, array_y_set_initial : false,
1, 6 1, 7
array_y_set_initial : false, array_y_set_initial : false,
1, 8 1, 9
array_y_set_initial : false, array_y_set_initial : false,
1, 10 1, 11
array_y_set_initial : false, array_y_set_initial : false,
1, 12 1, 13
array_y_set_initial : false, array_y_set_initial : false,
1, 14 1, 15
array_y_set_initial : false, array_y_set_initial : false,
1, 16 1, 17
array_y_set_initial : false, array_y_set_initial : false,
1, 18 1, 19
array_y_set_initial : false, array_y_set_initial : false,
1, 20 1, 21
array_y_set_initial : false, array_y_set_initial : false,
1, 22 1, 23
array_y_set_initial : false, array_y_set_initial : false,
1, 24 1, 25
array_y_set_initial : false, array_y_set_initial : false,
1, 26 1, 27
array_y_set_initial : false, array_y_set_initial : false,
1, 28 1, 29
array_y_set_initial : false, omniout_str(ALWAYS, "START of Optimize"),
1, 30
glob_check_sign : check_sign(x_start, x_end),
glob_h : check_sign(x_start, x_end), if glob_display_interval < glob_h
then glob_h : glob_display_interval, found_h : - 1.0, best_h : 0.0,
min_value : glob_large_float, est_answer : est_size_answer(), opt_iter : 1,
while (opt_iter <= 20) and (found_h < 0.0) do (omniout_int(ALWAYS,
"opt_iter", 32, opt_iter, 4, ""), array_x : x_start, array_x : glob_h,
1 2
glob_next_display : x_start, order_diff : 1, term_no : 1,
while term_no <= order_diff do (array_y :
term_no
array_y_init expt(glob_h, term_no - 1)
term_no
---------------------------------------------, term_no : 1 + term_no),
factorial_1(term_no - 1)
rows : order_diff, r_order : 1, while r_order <= rows do (term_no : 1,
while term_no <= 1 - r_order + rows do (it : - 1 + r_order + term_no,
array_y_init expt(glob_h, term_no - 1)
it
array_y_higher : ----------------------------------------,
r_order, term_no factorial_1(term_no - 1)
term_no : 1 + term_no), r_order : 1 + r_order), atomall(),
est_needed_step_err : estimated_needed_step_error(x_start, x_end, glob_h,
est_answer), omniout_float(ALWAYS, "est_needed_step_err", 32,
est_needed_step_err, 16, ""), value3 : test_suggested_h(),
omniout_float(ALWAYS, "value3", 32, value3, 32, ""),
if (value3 < est_needed_step_err) and (found_h < 0.0)
then (best_h : glob_h, found_h : 1.0),
omniout_float(ALWAYS, "best_h", 32, best_h, 32, ""), opt_iter : 1 + opt_iter,
glob_h : glob_h 0.5), if found_h > 0.0 then glob_h : best_h
else omniout_str(ALWAYS, "No increment to obtain desired accuracy found"),
if glob_html_log then html_log_file : openw("html/entry.html"),
if found_h > 0.0 then (omniout_str(ALWAYS, "START of Soultion"),
array_x : x_start, array_x : glob_h, glob_next_display : x_start,
1 2
order_diff : 1, term_no : 1, while term_no <=
order_diff do (array_y : (array_y_init expt(glob_h, term_no - 1))
term_no term_no
/factorial_1(term_no - 1), term_no : 1 + term_no), rows : order_diff,
r_order : 1, while r_order <= rows do (term_no : 1,
while term_no <= 1 - r_order + rows do (it : - 1 + r_order + term_no,
array_y_init expt(glob_h, term_no - 1)
it
array_y_higher : ----------------------------------------,
r_order, term_no factorial_1(term_no - 1)
term_no : 1 + term_no), r_order : 1 + r_order), current_iter : 1,
glob_clock_start_sec : elapsed_time_seconds(),
glob_log10normmin : - glob_large_float,
if omniabs(array_y_higher ) > glob_small_float
1, 1
then (tmp : omniabs(array_y_higher ), log10norm : log10(tmp),
1, 1
if log10norm < glob_log10normmin then glob_log10normmin : log10norm),
display_alot(current_iter), glob_clock_sec : elapsed_time_seconds(),
glob_current_iter : 0, glob_iter : 0, omniout_str(DEBUGL, " "),
glob_reached_optimal_h : true, glob_optimal_clock_start_sec :
elapsed_time_seconds(), while (glob_current_iter < glob_max_iter)
and (glob_check_sign array_x < glob_check_sign x_end)
1
and (convfloat(glob_clock_sec) - convfloat(glob_orig_start_sec) < convfloat(glob_max_sec)) do (if reached_interval
() then (omniout_str(INFO, " "), omniout_str(INFO, "TOP MAIN SOLVE Loop")),
glob_iter : 1 + glob_iter, glob_clock_sec : elapsed_time_seconds(),
glob_current_iter : 1 + glob_current_iter, atomall(),
display_alot(current_iter), if glob_look_poles then check_for_pole(),
if reached_interval() then glob_next_display :
glob_display_interval + glob_next_display, array_x : glob_h + array_x ,
1 1
array_x : glob_h, order_diff : 2, ord : 2, calc_term : 1,
2
iii : glob_max_terms, while iii >= calc_term do (array_y_higher_work :
2, iii
array_y_higher
2, iii
---------------------------
expt(glob_h, calc_term - 1)
-------------------------------------, iii : iii - 1), temp_sum : 0.0,
factorial_3(iii - calc_term, iii - 1)
ord : 2, calc_term : 1, iii : glob_max_terms,
while iii >= calc_term do (temp_sum : array_y_higher_work + temp_sum,
ord, iii
iii : iii - 1), array_y_higher_work2 :
ord, calc_term
temp_sum expt(glob_h, calc_term - 1)
------------------------------------, ord : 1, calc_term : 2,
factorial_1(calc_term - 1)
iii : glob_max_terms, while iii >= calc_term do (array_y_higher_work :
1, iii
array_y_higher
1, iii
---------------------------
expt(glob_h, calc_term - 1)
-------------------------------------, iii : iii - 1), temp_sum : 0.0,
factorial_3(iii - calc_term, iii - 1)
ord : 1, calc_term : 2, iii : glob_max_terms,
while iii >= calc_term do (temp_sum : array_y_higher_work + temp_sum,
ord, iii
iii : iii - 1), array_y_higher_work2 :
ord, calc_term
temp_sum expt(glob_h, calc_term - 1)
------------------------------------, ord : 1, calc_term : 1,
factorial_1(calc_term - 1)
iii : glob_max_terms, while iii >= calc_term do (array_y_higher_work :
1, iii
array_y_higher
1, iii
---------------------------
expt(glob_h, calc_term - 1)
-------------------------------------, iii : iii - 1), temp_sum : 0.0,
factorial_3(iii - calc_term, iii - 1)
ord : 1, calc_term : 1, iii : glob_max_terms,
while iii >= calc_term do (temp_sum : array_y_higher_work + temp_sum,
ord, iii
iii : iii - 1), array_y_higher_work2 :
ord, calc_term
temp_sum expt(glob_h, calc_term - 1)
------------------------------------, term_no : glob_max_terms,
factorial_1(calc_term - 1)
while term_no >= 1 do (array_y : array_y_higher_work2 ,
term_no 1, term_no
ord : 1, while ord <= order_diff do (array_y_higher :
ord, term_no
array_y_higher_work2 , ord : 1 + ord), term_no : term_no - 1)),
ord, term_no
omniout_str(ALWAYS, "Finished!"), if glob_iter >= glob_max_iter
then omniout_str(ALWAYS,
"Maximum Iterations Reached before Solution Completed!"),
if elapsed_time_seconds() - convfloat(glob_orig_start_sec) >=
convfloat(glob_max_sec) then omniout_str(ALWAYS,
"Maximum Time Reached before Solution Completed!"),
glob_clock_sec : elapsed_time_seconds(),
omniout_str(INFO, "diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"),
omniout_int(INFO, "Iterations ", 32, glob_iter, 4, " "),
prog_report(x_start, x_end), if glob_html_log
then (logstart(html_log_file), logitem_str(html_log_file,
"2013-01-13T00:13:13-06:00"), logitem_str(html_log_file, "Maxima"),
logitem_str(html_log_file, "lin_tanh"),
logitem_str(html_log_file, "diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"),
logitem_float(html_log_file, x_start), logitem_float(html_log_file, x_end),
logitem_float(html_log_file, array_x ), logitem_float(html_log_file, glob_h),
1
logitem_str(html_log_file, "16"), logitem_good_digits(html_log_file,
array_last_rel_error ), logitem_integer(html_log_file, glob_max_terms),
1
logitem_float(html_log_file, array_1st_rel_error ),
1
logitem_float(html_log_file, array_last_rel_error ),
1
logitem_integer(html_log_file, glob_iter),
logitem_pole(html_log_file, array_type_pole ),
1
if (array_type_pole = 1) or (array_type_pole = 2)
1 1
then (logitem_float(html_log_file, array_pole ),
1
logitem_float(html_log_file, array_pole ), 0)
2
else (logitem_str(html_log_file, "NA"), logitem_str(html_log_file, "NA"), 0),
logitem_time(html_log_file, convfloat(glob_clock_sec)),
if glob_percent_done < 100.0 then (logitem_time(html_log_file,
convfloat(glob_total_exp_sec)), 0) else (logitem_str(html_log_file, "Done"),
0), log_revs(html_log_file, " 156 | "), logitem_str(html_log_file, "lin_tanh diffeq.max"),
logitem_str(html_log_file,
"lin_tanh maxima results"),
logitem_str(html_log_file, "Languages compared - single equations"),
logend(html_log_file)), if glob_html_log then close(html_log_file)))
(%i57) main()
"##############ECHO OF PROBLEM#################"
"##############temp/lin_tanhpostode.ode#################"
"diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"
"!"
"/* BEGIN FIRST INPUT BLOCK */"
"Digits:32,"
"max_terms:30,"
"!"
"/* END FIRST INPUT BLOCK */"
"/* BEGIN SECOND INPUT BLOCK */"
"x_start:1.1,"
"x_end:2.0,"
"array_y_init[0 + 1] : exact_soln_y(x_start),"
"glob_h:0.00001,"
"glob_look_poles:true,"
"glob_max_iter:10,"
"/* 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 */"
"exact_soln_y (x) := (block("
" (ln(cosh(3.0*x + 1.0))/3.0) "
"));"
""
"/* END USER DEF BLOCK */"
"#######END OF ECHO OF PROBLEM#################"
"START of Optimize"
min_size = 0.0 ""
min_size = 1. ""
opt_iter = 1
glob_desired_digits_correct = 10. ""
desired_abs_gbl_error = 1.0000000000E-10 ""
range = 0.8999999999999999 ""
estimated_steps = 899.9999999999999 ""
step_error = 1.11111111111111120000000000000E-13 ""
est_needed_step_err = 1.11111111111111120000000000000E-13 ""
hn_div_ho = 0.5 ""
hn_div_ho_2 = 0.25 ""
hn_div_ho_3 = 0.125 ""
value3 = 4.1135837558932980000000000000000000000000000000000000000000000000000000000000000000000000000000000000E-85 ""
max_value3 = 4.1135837558932980000000000000000000000000000000000000000000000000000000000000000000000000000000000000E-85 ""
value3 = 4.1135837558932980000000000000000000000000000000000000000000000000000000000000000000000000000000000000E-85 ""
best_h = 1.000E-3 ""
"START of Soultion"
x[1] = 1.1 " "
y[1] (analytic) = 1.2023456360961104 " "
y[1] (numeric) = 1.2023456360961104 " "
absolute error = 0.0 " "
relative error = 0.0 "%"
Correct digits = 16
h = 1.000E-3 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1 " "
y[1] (analytic) = 1.2023456360961104 " "
y[1] (numeric) = 1.2023456360961104 " "
absolute error = 0.0 " "
relative error = 0.0 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5274316090036404 " "
Order of pole = 3.16433092264194470E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.101 " "
y[1] (analytic) = 1.203345269054324 " "
y[1] (numeric) = 1.2033452690543243 " "
absolute error = 2.2204460492503130000000000000000E-16 " "
relative error = 1.845227721712240400000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5284837627226338 " "
Order of pole = 3.318824382308349600E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1019999999999999 " "
y[1] (analytic) = 1.2043449042077943 " "
y[1] (numeric) = 1.2043449042077947 " "
absolute error = 4.4408920985006260000000000000000E-16 " "
relative error = 3.68739227690077660000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5295392324768 " "
Order of pole = 3.477937366910488500E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1029999999999998 " "
y[1] (analytic) = 1.205344541543394 " "
y[1] (numeric) = 1.2053445415433945 " "
absolute error = 4.4408920985006260000000000000000E-16 " "
relative error = 3.684334184493212600000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5305980549378 " "
Order of pole = 3.641724304361204600E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1039999999999996 " "
y[1] (analytic) = 1.2063441810480742 " "
y[1] (numeric) = 1.2063441810480748 " "
absolute error = 6.6613381477509390000000000000000E-16 " "
relative error = 5.52192173046630500000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5316602665679084 " "
Order of pole = 3.81023934306998770E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1049999999999995 " "
y[1] (analytic) = 1.2073438227088638 " "
y[1] (numeric) = 1.2073438227088646 " "
absolute error = 8.8817841970012520000000000000000E-16 " "
relative error = 7.35646634367465200000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5327259036303695 " "
Order of pole = 3.98353636683452800E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1059999999999994 " "
y[1] (analytic) = 1.2083434665128696 " "
y[1] (numeric) = 1.2083434665128707 " "
absolute error = 1.1102230246251565000000000000000E-15 " "
relative error = 9.1879755664928900000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5337950021541171 " "
Order of pole = 4.16166894597829200E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1069999999999993 " "
y[1] (analytic) = 1.2093431124472758 " "
y[1] (numeric) = 1.2093431124472769 " "
absolute error = 1.1102230246251565000000000000000E-15 " "
relative error = 9.18038076372274700000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5348675979421174 " "
Order of pole = 4.34469035060658370E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1079999999999992 " "
y[1] (analytic) = 1.2103427604993422 " "
y[1] (numeric) = 1.2103427604993433 " "
absolute error = 1.1102230246251565000000000000000E-15 " "
relative error = 9.1727984902980710000000000000E-14 "%"
Correct digits = 16
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.535943726560557 " "
Order of pole = 4.532653534114139400E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.108999999999999 " "
y[1] (analytic) = 1.2113424106564052 " "
y[1] (numeric) = 1.2113424106564068 " "
absolute error = 1.5543122344752192000000000000000E-15 " "
relative error = 1.28313202014694150000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5370234233187965 " "
Order of pole = 4.725611106932703400E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.109999999999999 " "
y[1] (analytic) = 1.2123420629058776 " "
y[1] (numeric) = 1.2123420629058794 " "
absolute error = 1.7763568394002505000000000000000E-15 " "
relative error = 1.46522742528827130000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5381067232719867 " "
Order of pole = 4.923615341032494500E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1109999999999989 " "
y[1] (analytic) = 1.2133417172352468 " "
y[1] (numeric) = 1.2133417172352485 " "
absolute error = 1.7763568394002505000000000000000E-15 " "
relative error = 1.4640202460423970000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.539193661206581 " "
Order of pole = 5.12671814851017900E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1119999999999988 " "
y[1] (analytic) = 1.214341373632075 " "
y[1] (numeric) = 1.2143413736320765 " "
absolute error = 1.5543122344752192000000000000000E-15 " "
relative error = 1.27996317034500500000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5402842716327803 " "
Order of pole = 5.33497107390861200E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1129999999999987 " "
y[1] (analytic) = 1.2153410320839981 " "
y[1] (numeric) = 1.2153410320840001 " "
absolute error = 1.9984014443252818000000000000000E-15 " "
relative error = 1.6443133174715050000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5413785887641425 " "
Order of pole = 5.54842526383350100E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1139999999999985 " "
y[1] (analytic) = 1.2163406925787275 " "
y[1] (numeric) = 1.2163406925787297 " "
absolute error = 2.220446049250313000000000000000E-15 " "
relative error = 1.825513248712260000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5424766465338642 " "
Order of pole = 5.76713149172185500E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1149999999999984 " "
y[1] (analytic) = 1.2173403551040471 " "
y[1] (numeric) = 1.2173403551040494 " "
absolute error = 2.220446049250313000000000000000E-15 " "
relative error = 1.82401416328675760000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5435784785526727 " "
Order of pole = 5.9911400993149400E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1159999999999983 " "
y[1] (analytic) = 1.2183400196478138 " "
y[1] (numeric) = 1.218340019647816 " "
absolute error = 2.220446049250313000000000000000E-15 " "
relative error = 1.82251753487682240000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5446841181275703 " "
Order of pole = 6.220501023072345000E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1169999999999982 " "
y[1] (analytic) = 1.2193396861979569 " "
y[1] (numeric) = 1.2193396861979593 " "
absolute error = 2.4424906541753444000000000000000E-15 " "
relative error = 2.00312569321131100000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5457935982395477 " "
Order of pole = 6.45526376330316500E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.117999999999998 " "
y[1] (analytic) = 1.2203393547424783 " "
y[1] (numeric) = 1.220339354742481 " "
absolute error = 2.6645352591003757000000000000000E-15 " "
relative error = 2.18343795006320870000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5469069515283547 " "
Order of pole = 6.69547736461382200E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.118999999999998 " "
y[1] (analytic) = 1.2213390252694523 " "
y[1] (numeric) = 1.2213390252694547 " "
absolute error = 2.4424906541753444000000000000000E-15 " "
relative error = 1.99984656482787900000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.548024210300398 " "
Order of pole = 6.94119042594536500E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1199999999999979 " "
y[1] (analytic) = 1.2223386977670228 " "
y[1] (numeric) = 1.2223386977670252 " "
absolute error = 2.4424906541753444000000000000000E-15 " "
relative error = 1.99821101846591630000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.54914540650461 " "
Order of pole = 7.1924510674712390E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1209999999999978 " "
y[1] (analytic) = 1.2233383722234055 " "
y[1] (numeric) = 1.2233383722234081 " "
absolute error = 2.6645352591003757000000000000000E-15 " "
relative error = 2.1780852457506170000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.550270571727666 " "
Order of pole = 7.44930692631076100E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1219999999999977 " "
y[1] (analytic) = 1.2243380486268867 " "
y[1] (numeric) = 1.2243380486268896 " "
absolute error = 2.886579864025407000000000000000E-15 " "
relative error = 2.35766573395538050000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.55139973718793 " "
Order of pole = 7.71180514617029200E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1229999999999976 " "
y[1] (analytic) = 1.2253377269658228 " "
y[1] (numeric) = 1.2253377269658257 " "
absolute error = 2.886579864025407000000000000000E-15 " "
relative error = 2.3557422582369570000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.552532933724637 " "
Order of pole = 7.97999236429642900E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1239999999999974 " "
y[1] (analytic) = 1.2263374072286393 " "
y[1] (numeric) = 1.2263374072286422 " "
absolute error = 2.886579864025407000000000000000E-15 " "
relative error = 2.35382191476055270000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5536701917814064 " "
Order of pole = 8.25391468817144900E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1249999999999973 " "
y[1] (analytic) = 1.227337089403831 " "
y[1] (numeric) = 1.227337089403834 " "
absolute error = 3.1086244689504383000000000000000E-15 " "
relative error = 2.53282044174223340000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5548115414109112 " "
Order of pole = 8.5336177025123300E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1259999999999972 " "
y[1] (analytic) = 1.2283367734799617 " "
y[1] (numeric) = 1.228336773479965 " "
absolute error = 3.3306690738754696000000000000000E-15 " "
relative error = 2.7115276085396817000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5559570122534743 " "
Order of pole = 8.81914643967771600E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1269999999999971 " "
y[1] (analytic) = 1.229336459445664 " "
y[1] (numeric) = 1.2293364594456673 " "
absolute error = 3.3306690738754696000000000000000E-15 " "
relative error = 2.7093226173226370000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5571066335370913 " "
Order of pole = 9.1105453818386200E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.127999999999997 " "
y[1] (analytic) = 1.2303361472896375 " "
y[1] (numeric) = 1.2303361472896408 " "
absolute error = 3.3306690738754696000000000000000E-15 " "
relative error = 2.7071212052191990000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5582604340633848 " "
Order of pole = 9.407858438840400E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.128999999999997 " "
y[1] (analytic) = 1.2313358370006495 " "
y[1] (numeric) = 1.231335837000653 " "
absolute error = 3.552713678800501000000000000000E-15 " "
relative error = 2.88525158778321740000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5594184421909505 " "
Order of pole = 9.71112892830721100E-2 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1299999999999968 " "
y[1] (analytic) = 1.2323355285675353 " "
y[1] (numeric) = 1.232335528567539 " "
absolute error = 3.774758283725532000000000000000E-15 " "
relative error = 3.06309296147073250000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5605806858460713 " "
Order of pole = 0.1002039958918779 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1309999999999967 " "
y[1] (analytic) = 1.2333352219791964 " "
y[1] (numeric) = 1.2333352219792002 " "
absolute error = 3.774758283725532000000000000000E-15 " "
relative error = 3.06061013782447900000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.561747192496718 " "
Order of pole = 0.10335712545387388 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1319999999999966 " "
y[1] (analytic) = 1.2343349172246005 " "
y[1] (numeric) = 1.2343349172246043 " "
absolute error = 3.774758283725532000000000000000E-15 " "
relative error = 3.05813133133515160000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5629179891492544 " "
Order of pole = 0.10657109303564738 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1329999999999965 " "
y[1] (analytic) = 1.2353346142927812 " "
y[1] (numeric) = 1.2353346142927852 " "
absolute error = 3.9968028886505635000000000000000E-15 " "
relative error = 3.2354010341875670000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5640931023395797 " "
Order of pole = 0.10984630739761769 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1339999999999963 " "
y[1] (analytic) = 1.236334313172838 " "
y[1] (numeric) = 1.2363343131728424 " "
absolute error = 4.440892098500626000000000000000E-15 " "
relative error = 3.59198321294168900000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5652725581284679 " "
Order of pole = 0.11318317093347652 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1349999999999962 " "
y[1] (analytic) = 1.2373340138539357 " "
y[1] (numeric) = 1.2373340138539402 " "
absolute error = 4.440892098500626000000000000000E-15 " "
relative error = 3.58908107978745200000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5664563820752926 " "
Order of pole = 0.11658207930511644 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1359999999999961 " "
y[1] (analytic) = 1.2383337163253036 " "
y[1] (numeric) = 1.238333716325308 " "
absolute error = 4.440892098500626000000000000000E-15 " "
relative error = 3.5861836272041125000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5676445992575168 " "
Order of pole = 0.12004342172343385 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.136999999999996 " "
y[1] (analytic) = 1.2393334205762352 " "
y[1] (numeric) = 1.2393334205762399 " "
absolute error = 4.6629367034256575000000000000000E-15 " "
relative error = 3.7624553861039250000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5688372342429104 " "
Order of pole = 0.12356758055689099 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.137999999999996 " "
y[1] (analytic) = 1.2403331265960882 " "
y[1] (numeric) = 1.240333126596093 " "
absolute error = 4.884981308350689000000000000000E-15 " "
relative error = 3.9384429905188467000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5700343110747907 " "
Order of pole = 0.12715493113855736 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1389999999999958 " "
y[1] (analytic) = 1.2413328343742842 " "
y[1] (numeric) = 1.241332834374289 " "
absolute error = 4.884981308350689000000000000000E-15 " "
relative error = 3.9352711642506820000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.571235853290182 " "
Order of pole = 0.13080584200569234 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1399999999999957 " "
y[1] (analytic) = 1.242332543900308 " "
y[1] (numeric) = 1.2423325439003126 " "
absolute error = 4.6629367034256575000000000000000E-15 " "
relative error = 3.75337241732905840000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5724418838789669 " "
Order of pole = 0.13452067433972204 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1409999999999956 " "
y[1] (analytic) = 1.243332255163706 " "
y[1] (numeric) = 1.243332255163711 " "
absolute error = 4.884981308350689000000000000000E-15 " "
relative error = 3.92894279711861640000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5736524253043438 " "
Order of pole = 0.13829978225833983 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1419999999999955 " "
y[1] (analytic) = 1.244331968154089 " "
y[1] (numeric) = 1.2443319681540943 " "
absolute error = 5.329070518200751000000000000000E-15 " "
relative error = 4.28267588922125850000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5748674994629126 " "
Order of pole = 0.14214351225650645 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1429999999999954 " "
y[1] (analytic) = 1.2453316828611292 " "
y[1] (numeric) = 1.2453316828611345 " "
absolute error = 5.329070518200751000000000000000E-15 " "
relative error = 4.2792378862130120000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5760871277031536 " "
Order of pole = 0.14605220345101522 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1439999999999952 " "
y[1] (analytic) = 1.2463313992745602 " "
y[1] (numeric) = 1.2463313992745655 " "
absolute error = 5.329070518200751000000000000000E-15 " "
relative error = 4.2758053927732150000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5773113308244793 " "
Order of pole = 0.1500261876062563 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1449999999999951 " "
y[1] (analytic) = 1.2473311173841768 " "
y[1] (numeric) = 1.2473311173841823 " "
absolute error = 5.551115123125783000000000000000E-15 " "
relative error = 4.45039416219105250000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5785401290208598 " "
Order of pole = 0.1540657882961174 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.145999999999995 " "
y[1] (analytic) = 1.2483308371798356 " "
y[1] (numeric) = 1.2483308371798414 " "
absolute error = 5.773159728050814000000000000000E-15 " "
relative error = 4.624703288667639000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5797735419246206 " "
Order of pole = 0.1581713215549172 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.146999999999995 " "
y[1] (analytic) = 1.2493305586514536 " "
y[1] (numeric) = 1.2493305586514594 " "
absolute error = 5.773159728050814000000000000000E-15 " "
relative error = 4.6210025745967910000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5810115885757097 " "
Order of pole = 0.16234309543805026 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1479999999999948 " "
y[1] (analytic) = 1.2503302817890078 " "
y[1] (numeric) = 1.2503302817890136 " "
absolute error = 5.773159728050814000000000000000E-15 " "
relative error = 4.6173077723035016000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5822542874115313 " "
Order of pole = 0.16658140987013859 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1489999999999947 " "
y[1] (analytic) = 1.251330006582535 " "
y[1] (numeric) = 1.251330006582541 " "
absolute error = 5.995204332975845000000000000000E-15 " "
relative error = 4.7910657471957740000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5835016562570863 " "
Order of pole = 0.17088655651821583 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1499999999999946 " "
y[1] (analytic) = 1.2523297330221321 " "
y[1] (numeric) = 1.2523297330221383 " "
absolute error = 6.217248937900877000000000000000E-15 " "
relative error = 4.9645462963634673000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5847537123333257 " "
Order of pole = 0.1752588189089863 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1509999999999945 " "
y[1] (analytic) = 1.253329461097955 " "
y[1] (numeric) = 1.2533294610979613 " "
absolute error = 6.217248937900877000000000000000E-15 " "
relative error = 4.9605862870680273000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.586010472232364 " "
Order of pole = 0.17969847205891654 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1519999999999944 " "
y[1] (analytic) = 1.2543291908002183 " "
y[1] (numeric) = 1.2543291908002245 " "
absolute error = 6.217248937900877000000000000000E-15 " "
relative error = 4.9566325837753070000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5872719519165297 " "
Order of pole = 0.18420578249215325 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1529999999999943 " "
y[1] (analytic) = 1.2553289221191948 " "
y[1] (numeric) = 1.2553289221192014 " "
absolute error = 6.661338147750939000000000000000E-15 " "
relative error = 5.3064483980067490000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.588538166708287 " "
Order of pole = 0.18878100806523435 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1539999999999941 " "
y[1] (analytic) = 1.2563286550452164 " "
y[1] (numeric) = 1.256328655045223 " "
absolute error = 6.661338147750939000000000000000E-15 " "
relative error = 5.3022257519957560000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5898091312866967 " "
Order of pole = 0.19342439795089028 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.154999999999994 " "
y[1] (analytic) = 1.257328389568672 " "
y[1] (numeric) = 1.2573283895686789 " "
absolute error = 6.8833827526759700000000000000000E-15 " "
relative error = 5.474610141458209000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5910848596731622 " "
Order of pole = 0.1981361923973992 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.155999999999994 " "
y[1] (analytic) = 1.2583281256800085 " "
y[1] (numeric) = 1.2583281256800154 " "
absolute error = 6.8833827526759700000000000000000E-15 " "
relative error = 5.4702605879973850000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5923653652225238 " "
Order of pole = 0.202916622627761 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1569999999999938 " "
y[1] (analytic) = 1.2593278633697291 " "
y[1] (numeric) = 1.2593278633697362 " "
absolute error = 7.105427357601002000000000000000E-15 " "
relative error = 5.6422378669429170000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.593650660631642 " "
Order of pole = 0.2077659109517498 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1579999999999937 " "
y[1] (analytic) = 1.2603276026283947 " "
y[1] (numeric) = 1.260327602628402 " "
absolute error = 7.327471962526033000000000000000E-15 " "
relative error = 5.8139423013863210000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5949407579247021 " "
Order of pole = 0.2126842705498042 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1589999999999936 " "
y[1] (analytic) = 1.2613273434466226 " "
y[1] (numeric) = 1.26132734344663 " "
absolute error = 7.327471962526033000000000000000E-15 " "
relative error = 5.8093341118757090000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5962356684113983 " "
Order of pole = 0.21767190487801358 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1599999999999935 " "
y[1] (analytic) = 1.2623270858150855 " "
y[1] (numeric) = 1.2623270858150926 " "
absolute error = 7.105427357601002000000000000000E-15 " "
relative error = 5.6288322079479280000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5975354027373159 " "
Order of pole = 0.22272900838598275 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1609999999999934 " "
y[1] (analytic) = 1.2633268297245117 " "
y[1] (numeric) = 1.263326829724519 " "
absolute error = 7.327471962526033000000000000000E-15 " "
relative error = 5.8001395918456860000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.5988399708369951 " "
Order of pole = 0.22785576584555045 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1619999999999933 " "
y[1] (analytic) = 1.264326575165686 " "
y[1] (numeric) = 1.2643265751656936 " "
absolute error = 7.549516567451064000000000000000E-15 " "
relative error = 5.9711760519324090000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6001493819526549 " "
Order of pole = 0.2330523526027939 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1629999999999932 " "
y[1] (analytic) = 1.2653263221294477 " "
y[1] (numeric) = 1.2653263221294555 " "
absolute error = 7.771561172376096000000000000000E-15 " "
relative error = 6.1419422298092650000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.601463644590359 " "
Order of pole = 0.23831893396578074 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.163999999999993 " "
y[1] (analytic) = 1.2663260706066912 " "
y[1] (numeric) = 1.266326070606699 " "
absolute error = 7.771561172376096000000000000000E-15 " "
relative error = 6.1370932438062940000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6027827665522463 " "
Order of pole = 0.24365566564109464 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.164999999999993 " "
y[1] (analytic) = 1.2673258205883646 " "
y[1] (numeric) = 1.2673258205883724 " "
absolute error = 7.771561172376096000000000000000E-15 " "
relative error = 6.1322519009106080000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6041067549078405 " "
Order of pole = 0.2490626933402691 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1659999999999928 " "
y[1] (analytic) = 1.2683255720654707 " "
y[1] (numeric) = 1.2683255720654787 " "
absolute error = 7.993605777301127000000000000000E-15 " "
relative error = 6.3024872740549780000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6054356159924859 " "
Order of pole = 0.2545401527341138 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1669999999999927 " "
y[1] (analytic) = 1.2693253250290661 " "
y[1] (numeric) = 1.2693253250290741 " "
absolute error = 7.993605777301127000000000000000E-15 " "
relative error = 6.2975232745144220000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.606769355399838 " "
Order of pole = 0.2600881693434616 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1679999999999926 " "
y[1] (analytic) = 1.270325079470261 " "
y[1] (numeric) = 1.270325079470269 " "
absolute error = 7.993605777301127000000000000000E-15 " "
relative error = 6.2925670810455430000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6081079779704828 " "
Order of pole = 0.2657068583724769 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1689999999999925 " "
y[1] (analytic) = 1.2713248353802178 " "
y[1] (numeric) = 1.271324835380226 " "
absolute error = 8.215650382226158000000000000000E-15 " "
relative error = 6.4622747496072360000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6094514877956991 " "
Order of pole = 0.27139632475709696 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1699999999999924 " "
y[1] (analytic) = 1.2723245927501527 " "
y[1] (numeric) = 1.2723245927501614 " "
absolute error = 8.659739592076221000000000000000E-15 " "
relative error = 6.8062345422075320000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6107998882032881 " "
Order of pole = 0.2771566629578004 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1709999999999923 " "
y[1] (analytic) = 1.2733243515713348 " "
y[1] (numeric) = 1.2733243515713435 " "
absolute error = 8.659739592076221000000000000000E-15 " "
relative error = 6.8008905832906940000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6121531817620254 " "
Order of pole = 0.28298795700101564 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1719999999999922 " "
y[1] (analytic) = 1.2743241118350843 " "
y[1] (numeric) = 1.274324111835093 " "
absolute error = 8.659739592076221000000000000000E-15 " "
relative error = 6.7955550017850670000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6135113702459498 " "
Order of pole = 0.28889027997091254 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.172999999999992 " "
y[1] (analytic) = 1.2753238735327737 " "
y[1] (numeric) = 1.2753238735327825 " "
absolute error = 8.881784197001252000000000000000E-15 " "
relative error = 6.9643361826183260000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6148744546695766 " "
Order of pole = 0.2948636945065406 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.173999999999992 " "
y[1] (analytic) = 1.276323636655827 " "
y[1] (numeric) = 1.2763236366558361 " "
absolute error = 9.103828801926284000000000000000E-15 " "
relative error = 7.1328529382874850000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6162424352429188 " "
Order of pole = 0.30090825214708694 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1749999999999918 " "
y[1] (analytic) = 1.2773234011957204 " "
y[1] (numeric) = 1.2773234011957293 " "
absolute error = 8.881784197001252000000000000000E-15 " "
relative error = 6.9534341801668150000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6176153114006937 " "
Order of pole = 0.30702399373076616 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1759999999999917 " "
y[1] (analytic) = 1.2783231671439794 " "
y[1] (numeric) = 1.2783231671439885 " "
absolute error = 9.103828801926284000000000000000E-15 " "
relative error = 7.1216958558812580000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6189930817627558 " "
Order of pole = 0.3132109488159802 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1769999999999916 " "
y[1] (analytic) = 1.2793229344921815 " "
y[1] (numeric) = 1.2793229344921908 " "
absolute error = 9.325873406851315000000000000000E-15 " "
relative error = 7.2896945371757570000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6203757441509845 " "
Order of pole = 0.3194691359348152 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1779999999999915 " "
y[1] (analytic) = 1.280322703231954 " "
y[1] (numeric) = 1.2803227032319635 " "
absolute error = 9.547918011776346000000000000000E-15 " "
relative error = 7.4574308396424380000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6217632955743095 " "
Order of pole = 0.325798562341582 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1789999999999914 " "
y[1] (analytic) = 1.2813224733549748 " "
y[1] (numeric) = 1.2813224733549842 " "
absolute error = 9.325873406851315000000000000000E-15 " "
relative error = 7.2783187689143850000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6231557322250503 " "
Order of pole = 0.33219922397166 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1799999999999913 " "
y[1] (analytic) = 1.2823222448529707 " "
y[1] (numeric) = 1.28232224485298 " "
absolute error = 9.325873406851315000000000000000E-15 " "
relative error = 7.2726441768313910000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6245530494757816 " "
Order of pole = 0.3386711053465188 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1809999999999912 " "
y[1] (analytic) = 1.2833220177177185 " "
y[1] (numeric) = 1.283322017717728 " "
absolute error = 9.547918011776346000000000000000E-15 " "
relative error = 7.4400017142669500000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6259552418723844 " "
Order of pole = 0.34521417950535316 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.181999999999991 " "
y[1] (analytic) = 1.284321791941044 " "
y[1] (numeric) = 1.284321791941054 " "
absolute error = 9.992007221626409000000000000000E-15 " "
relative error = 7.7799872931573580000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6273623031219095 " "
Order of pole = 0.3518284077855274 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.182999999999991 " "
y[1] (analytic) = 1.2853215675148224 " "
y[1] (numeric) = 1.2853215675148324 " "
absolute error = 9.992007221626409000000000000000E-15 " "
relative error = 7.7739357015116620000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6287742260909996 " "
Order of pole = 0.3585137398039393 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1839999999999908 " "
y[1] (analytic) = 1.286321344430977 " "
y[1] (numeric) = 1.286321344430987 " "
absolute error = 9.992007221626409000000000000000E-15 " "
relative error = 7.767893508792330000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.630191002807852 " "
Order of pole = 0.3652701134433318 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1849999999999907 " "
y[1] (analytic) = 1.2873211226814785 " "
y[1] (numeric) = 1.287321122681489 " "
absolute error = 1.043609643147647100000000000000E-14 " "
relative error = 8.1068322795311360000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6316126244622657 " "
Order of pole = 0.3720974548799525 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1859999999999906 " "
y[1] (analytic) = 1.2883209022583475 " "
y[1] (numeric) = 1.2883209022583582 " "
absolute error = 1.065814103640150300000000000000E-14 " "
relative error = 8.2728930484000030000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.633039081354331 " "
Order of pole = 0.3789956777905896 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1869999999999905 " "
y[1] (analytic) = 1.2893206831536517 " "
y[1] (numeric) = 1.289320683153662 " "
absolute error = 1.043609643147647100000000000000E-14 " "
relative error = 8.0942596887145210000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6344703629596868 " "
Order of pole = 0.3859646842956863 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1879999999999904 " "
y[1] (analytic) = 1.290320465359505 " "
y[1] (numeric) = 1.2903204653595155 " "
absolute error = 1.043609643147647100000000000000E-14 " "
relative error = 8.0879879934081320000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6359064578853955 " "
Order of pole = 0.3930043642942245 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1889999999999903 " "
y[1] (analytic) = 1.2913202488680695 " "
y[1] (numeric) = 1.2913202488680802 " "
absolute error = 1.065814103640150300000000000000E-14 " "
relative error = 8.253677618502531000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6373473538507377 " "
Order of pole = 0.40011459517417336 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1899999999999902 " "
y[1] (analytic) = 1.2923200336715541 " "
y[1] (numeric) = 1.2923200336715652 " "
absolute error = 1.110223024625156500000000000000E-14 " "
relative error = 8.5909294578599870000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.638793037717229 " "
Order of pole = 0.40729524224870595 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.19099999999999 " "
y[1] (analytic) = 1.2933198197622147 " "
y[1] (numeric) = 1.2933198197622258 " "
absolute error = 1.110223024625156500000000000000E-14 " "
relative error = 8.5842883381256640000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6402434954664187 " "
Order of pole = 0.4145461583800376 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.19199999999999 " "
y[1] (analytic) = 1.2943196071323526 " "
y[1] (numeric) = 1.2943196071323637 " "
absolute error = 1.110223024625156500000000000000E-14 " "
relative error = 8.5776574696641300000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6416987121907674 " "
Order of pole = 0.42186718384832567 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1929999999999898 " "
y[1] (analytic) = 1.2953193957743154 " "
y[1] (numeric) = 1.2953193957743265 " "
absolute error = 1.110223024625156500000000000000E-14 " "
relative error = 8.5710368288084490000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6431586721125855 " "
Order of pole = 0.4292581465961316 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1939999999999897 " "
y[1] (analytic) = 1.2963191856804965 " "
y[1] (numeric) = 1.296319185680508 " "
absolute error = 1.154631945610162800000000000000E-14 " "
relative error = 8.9070034476427520000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6446233585396046 " "
Order of pole = 0.43671886158249507 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1949999999999896 " "
y[1] (analytic) = 1.2973189768433355 " "
y[1] (numeric) = 1.297318976843347 " "
absolute error = 1.154631945610162800000000000000E-14 " "
relative error = 8.9001391810334730000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6460927539022325 " "
Order of pole = 0.4442491312853125 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1959999999999895 " "
y[1] (analytic) = 1.2983187692553164 " "
y[1] (numeric) = 1.2983187692553277 " "
absolute error = 1.132427485117659700000000000000E-14 " "
relative error = 8.7222607570188030000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6475668397334036 " "
Order of pole = 0.4518487453930735 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1969999999999894 " "
y[1] (analytic) = 1.2993185629089679 " "
y[1] (numeric) = 1.2993185629089794 " "
absolute error = 1.154631945610162800000000000000E-14 " "
relative error = 8.8864423134625680000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6490455966478552 " "
Order of pole = 0.4595174804825639 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1979999999999893 " "
y[1] (analytic) = 1.3003183577968638 " "
y[1] (numeric) = 1.3003183577968758 " "
absolute error = 1.19904086659516900000000000000E-14 " "
relative error = 9.2211331125610680000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6505290043578456 " "
Order of pole = 0.4672551002220633 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.1989999999999892 " "
y[1] (analytic) = 1.3013181539116232 " "
y[1] (numeric) = 1.301318153911635 " "
absolute error = 1.17683640610266600000000000000E-14 " "
relative error = 9.0434180339774840000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6520170416935829 " "
Order of pole = 0.47506135564884566 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.199999999999989 " "
y[1] (analytic) = 1.3023179512459075 " "
y[1] (numeric) = 1.3023179512459195 " "
absolute error = 1.19904086659516900000000000000E-14 " "
relative error = 9.206974882347780000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.653509686533083 " "
Order of pole = 0.48293598412363714 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.200999999999989 " "
y[1] (analytic) = 1.3033177497924233 " "
y[1] (numeric) = 1.3033177497924355 " "
absolute error = 1.221245327087672200000000000000E-14 " "
relative error = 9.3702807874915960000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6550069158697756 " "
Order of pole = 0.49087871028414654 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2019999999999889 " "
y[1] (analytic) = 1.3043175495439203 " "
y[1] (numeric) = 1.3043175495439328 " "
absolute error = 1.243449787580175300000000000000E-14 " "
relative error = 9.5333363260731350000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.656508705752891 " "
Order of pole = 0.49888924516909583 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2029999999999887 " "
y[1] (analytic) = 1.305317350493192 " "
y[1] (numeric) = 1.3053173504932045 " "
absolute error = 1.243449787580175300000000000000E-14 " "
relative error = 9.5260343173279570000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.658015031314746 " "
Order of pole = 0.5069672865799735 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2039999999999886 " "
y[1] (analytic) = 1.3063171526330748 " "
y[1] (numeric) = 1.306317152633087 " "
absolute error = 1.221245327087672200000000000000E-14 " "
relative error = 9.3487659151230780000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6595258667788118 " "
Order of pole = 0.5151125191451147 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2049999999999885 " "
y[1] (analytic) = 1.307316955956447 " "
y[1] (numeric) = 1.3073169559564595 " "
absolute error = 1.243449787580175300000000000000E-14 " "
relative error = 9.5114637801852280000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6610411854280271 " "
Order of pole = 0.5233246139001491 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2059999999999884 " "
y[1] (analytic) = 1.308316760456231 " "
y[1] (numeric) = 1.3083167604562436 " "
absolute error = 1.265654248072678500000000000000E-14 " "
relative error = 9.6739129722019670000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6625609596049646 " "
Order of pole = 0.5316032281926653 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2069999999999883 " "
y[1] (analytic) = 1.309316566125391 " "
y[1] (numeric) = 1.3093165661254034 " "
absolute error = 1.243449787580175300000000000000E-14 " "
relative error = 9.4969377135421670000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6640851607408358 " "
Order of pole = 0.5399480061126276 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2079999999999882 " "
y[1] (analytic) = 1.3103163729569325 " "
y[1] (numeric) = 1.310316372956945 " "
absolute error = 1.243449787580175300000000000000E-14 " "
relative error = 9.4896912932114060000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.665613759309192 " "
Order of pole = 0.5483585777894628 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.208999999999988 " "
y[1] (analytic) = 1.3113161809439038 " "
y[1] (numeric) = 1.3113161809439162 " "
absolute error = 1.243449787580175300000000000000E-14 " "
relative error = 9.4824559145233970000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6671467247058285 " "
Order of pole = 0.5568345576192915 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.209999999999988 " "
y[1] (analytic) = 1.3123159900793941 " "
y[1] (numeric) = 1.3123159900794068 " "
absolute error = 1.265654248072678500000000000000E-14 " "
relative error = 9.6444321157445270000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6686840260033298 " "
Order of pole = 0.5653755551636337 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2109999999999879 " "
y[1] (analytic) = 1.3133158003565348 " "
y[1] (numeric) = 1.3133158003565475 " "
absolute error = 1.265654248072678500000000000000E-14 " "
relative error = 9.6370899347215850000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6702256303988372 " "
Order of pole = 0.5739811526374723 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2119999999999878 " "
y[1] (analytic) = 1.3143156117684973 " "
y[1] (numeric) = 1.3143156117685102 " "
absolute error = 1.287858708565181600000000000000E-14 " "
relative error = 9.7987020547696590000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6717715047628012 " "
Order of pole = 0.5826509273204721 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2129999999999876 " "
y[1] (analytic) = 1.3153154243084944 " "
y[1] (numeric) = 1.3153154243085075 " "
absolute error = 1.310063169057684700000000000000E-14 " "
relative error = 9.9600684736623470000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6733216148768577 " "
Order of pole = 0.5913844404748545 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2139999999999875 " "
y[1] (analytic) = 1.3163152379697793 " "
y[1] (numeric) = 1.3163152379697927 " "
absolute error = 1.332267629550187800000000000000E-14 " "
relative error = 1.0121189750906574000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6748759255603594 " "
Order of pole = 0.6001812391367078 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2149999999999874 " "
y[1] (analytic) = 1.3173150527456463 " "
y[1] (numeric) = 1.3173150527456596 " "
absolute error = 1.332267629550187800000000000000E-14 " "
relative error = 1.0113507978014646000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6764344006997611 " "
Order of pole = 0.609040856562105 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2159999999999873 " "
y[1] (analytic) = 1.3183148686294288 " "
y[1] (numeric) = 1.318314868629442 " "
absolute error = 1.310063169057684700000000000000E-14 " "
relative error = 9.9374072175919340000000000000E-13 "%"
Correct digits = 15
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6779970032209983 " "
Order of pole = 0.6179628117224407 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2169999999999872 " "
y[1] (analytic) = 1.3193146856145 " "
y[1] (numeric) = 1.3193146856145133 " "
absolute error = 1.332267629550187800000000000000E-14 " "
relative error = 1.0098179335657548000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.679563695102282 " "
Order of pole = 0.626946609515528 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.217999999999987 " "
y[1] (analytic) = 1.3203145036942734 " "
y[1] (numeric) = 1.320314503694287 " "
absolute error = 1.35447209004269100000000000000E-14 " "
relative error = 1.0258707953694698000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6811344373866617 " "
Order of pole = 0.6359917408819467 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.218999999999987 " "
y[1] (analytic) = 1.321314322862202 " "
y[1] (numeric) = 1.3213143228622155 " "
absolute error = 1.35447209004269100000000000000E-14 " "
relative error = 1.0250945339854209000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6827091901314348 " "
Order of pole = 0.6450976820471279 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2199999999999869 " "
y[1] (analytic) = 1.3223141431117773 " "
y[1] (numeric) = 1.3223141431117909 " "
absolute error = 1.35447209004269100000000000000E-14 " "
relative error = 1.0243194456464308000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6842879124685646 " "
Order of pole = 0.6542638953844708 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2209999999999868 " "
y[1] (analytic) = 1.32331396443653 " "
y[1] (numeric) = 1.3233139644365437 " "
absolute error = 1.376676550535194000000000000000E-14 " "
relative error = 1.040324962580884000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.68587056257195 " "
Order of pole = 0.6634898288584878 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2219999999999867 " "
y[1] (analytic) = 1.3243137868300294 " "
y[1] (numeric) = 1.3243137868300434 " "
absolute error = 1.398881011027697200000000000000E-14 " "
relative error = 1.056306311192423000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6874570976496452 " "
Order of pole = 0.6727749159454053 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2229999999999865 " "
y[1] (analytic) = 1.3253136102858836 " "
y[1] (numeric) = 1.3253136102858978 " "
absolute error = 1.421085471520200400000000000000E-14 " "
relative error = 1.0722635461456234000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6890474739746866 " "
Order of pole = 0.682118575965216 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2239999999999864 " "
y[1] (analytic) = 1.3263134347977388 " "
y[1] (numeric) = 1.3263134347977528 " "
absolute error = 1.398881011027697200000000000000E-14 " "
relative error = 1.0547137458809086000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6906416468535992 " "
Order of pole = 0.6915202136666991 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2249999999999863 " "
y[1] (analytic) = 1.3273132603592785 " "
y[1] (numeric) = 1.3273132603592928 " "
absolute error = 1.421085471520200400000000000000E-14 " "
relative error = 1.0706481385830044000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6922395706515956 " "
Order of pole = 0.7009792195050171 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2259999999999862 " "
y[1] (analytic) = 1.3283130869642248 " "
y[1] (numeric) = 1.3283130869642394 " "
absolute error = 1.465494392505206600000000000000E-14 " "
relative error = 1.1032748279658232000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6938411987687598 " "
Order of pole = 0.7104949692853122 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.226999999999986 " "
y[1] (analytic) = 1.3293129146063374 " "
y[1] (numeric) = 1.329312914606352 " "
absolute error = 1.465494392505206600000000000000E-14 " "
relative error = 1.1024450123086317000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6954464836646932 " "
Order of pole = 0.7200668244823127 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.227999999999986 " "
y[1] (analytic) = 1.3303127432794128 " "
y[1] (numeric) = 1.3303127432794275 " "
absolute error = 1.465494392505206600000000000000E-14 " "
relative error = 1.1016164431324257000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6970553768316134 " "
Order of pole = 0.72969413180701 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2289999999999859 " "
y[1] (analytic) = 1.3313125729772846 " "
y[1] (numeric) = 1.3313125729772994 " "
absolute error = 1.487698852997709800000000000000E-14 " "
relative error = 1.1174677406303542000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.6986678288311303 " "
Order of pole = 0.7393762236958263 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2299999999999858 " "
y[1] (analytic) = 1.3323124036938234 " "
y[1] (numeric) = 1.3323124036938385 " "
absolute error = 1.50990331349021300000000000000E-14 " "
relative error = 1.1332952461479909000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7002837892635498 " "
Order of pole = 0.749112417851201 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2309999999999857 " "
y[1] (analytic) = 1.3333122354229372 " "
y[1] (numeric) = 1.3333122354229523 " "
absolute error = 1.50990331349021300000000000000E-14 " "
relative error = 1.1324454042914109000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.701903206788809 " "
Order of pole = 0.7589020175015424 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2319999999999856 " "
y[1] (analytic) = 1.3343120681585692 " "
y[1] (numeric) = 1.3343120681585843 " "
absolute error = 1.50990331349021300000000000000E-14 " "
relative error = 1.131596835194611000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.703526029107744 " "
Order of pole = 0.7687443110786969 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2329999999999854 " "
y[1] (analytic) = 1.3353119018946993 " "
y[1] (numeric) = 1.3353119018947146 " "
absolute error = 1.53210777398271600000000000000E-14 " "
relative error = 1.1473782056527614000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7051522030057364 " "
Order of pole = 0.7786385728329197 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2339999999999853 " "
y[1] (analytic) = 1.3363117366253434 " "
y[1] (numeric) = 1.3363117366253592 " "
absolute error = 1.576516694967722300000000000000E-14 " "
relative error = 1.1797521878756978000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7067816742901893 " "
Order of pole = 0.7885840618620801 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2349999999999852 " "
y[1] (analytic) = 1.337311572344554 " "
y[1] (numeric) = 1.3373115723445697 " "
absolute error = 1.576516694967722300000000000000E-14 " "
relative error = 1.1788701508084595000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7084143878503522 " "
Order of pole = 0.798580022969956 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.235999999999985 " "
y[1] (analytic) = 1.338311409046418 " "
y[1] (numeric) = 1.3383114090464336 " "
absolute error = 1.55431223447521920000000000000E-14 " "
relative error = 1.1613980303602935000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7100502876296857 " "
Order of pole = 0.8086256862244152 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.236999999999985 " "
y[1] (analytic) = 1.3393112467250576 " "
y[1] (numeric) = 1.3393112467250736 " "
absolute error = 1.598721155460225400000000000000E-14 " "
relative error = 1.193689039324868100000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7116893166346578 " "
Order of pole = 0.8187202670455438 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.237999999999985 " "
y[1] (analytic) = 1.3403110853746314 " "
y[1] (numeric) = 1.3403110853746474 " "
absolute error = 1.598721155460225400000000000000E-14 " "
relative error = 1.192798577065686000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7133314169455554 " "
Order of pole = 0.8288629662823492 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2389999999999848 " "
y[1] (analytic) = 1.3413109249893316 " "
y[1] (numeric) = 1.3413109249893478 " "
absolute error = 1.620925615952728500000000000000E-14 " "
relative error = 1.2084637392822406000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7149765297058737 " "
Order of pole = 0.8390529701011413 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2399999999999847 " "
y[1] (analytic) = 1.3423107655633864 " "
y[1] (numeric) = 1.3423107655634023 " "
absolute error = 1.598721155460225400000000000000E-14 " "
relative error = 1.191021629621826100000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7166245951395762 " "
Order of pole = 0.8492894501254238 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2409999999999846 " "
y[1] (analytic) = 1.343310607091057 " "
y[1] (numeric) = 1.3433106070910732 " "
absolute error = 1.620925615952728500000000000000E-14 " "
relative error = 1.2066647932326297000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7182755525364108 " "
Order of pole = 0.8595715632451686 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2419999999999844 " "
y[1] (analytic) = 1.34431044956664 " "
y[1] (numeric) = 1.3443104495666565 " "
absolute error = 1.643130076445231700000000000000E-14 " "
relative error = 1.2222846865282651000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7199293402684732 " "
Order of pole = 0.8698984517654598 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2429999999999843 " "
y[1] (analytic) = 1.3453102929844662 " "
y[1] (numeric) = 1.3453102929844827 " "
absolute error = 1.643130076445231700000000000000E-14 " "
relative error = 1.2213762765466361000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7215858957988623 " "
Order of pole = 0.8802692435180468 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2439999999999842 " "
y[1] (analytic) = 1.3463101373389 " "
y[1] (numeric) = 1.3463101373389161 " "
absolute error = 1.620925615952728500000000000000E-14 " "
relative error = 1.2039763877561156000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.723245155671108 " "
Order of pole = 0.890683051692843 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2449999999999841 " "
y[1] (analytic) = 1.3473099826243384 " "
y[1] (numeric) = 1.347309982624355 " "
absolute error = 1.665334536937734800000000000000E-14 " "
relative error = 1.2360440866725687000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.72490705552708 " "
Order of pole = 0.9011389750145913 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.245999999999984 " "
y[1] (analytic) = 1.348309828835214 " "
y[1] (numeric) = 1.3483098288352309 " "
absolute error = 1.68753899743023800000000000000E-14 " "
relative error = 1.2515958582665523000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7265715300897215 " "
Order of pole = 0.9116360975041573 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.246999999999984 " "
y[1] (analytic) = 1.3493096759659915 " "
y[1] (numeric) = 1.3493096759660084 " "
absolute error = 1.68753899743023800000000000000E-14 " "
relative error = 1.2506684177018912000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7282385131888793 " "
Order of pole = 0.9221734887639279 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2479999999999838 " "
y[1] (analytic) = 1.3503095240111689 " "
y[1] (numeric) = 1.3503095240111855 " "
absolute error = 1.665334536937734800000000000000E-14 " "
relative error = 1.2332983714658004000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7299079377679711 " "
Order of pole = 0.9327502040769797 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2489999999999837 " "
y[1] (analytic) = 1.3513093729652763 " "
y[1] (numeric) = 1.3513093729652934 " "
absolute error = 1.70974345792274100000000000000E-14 " "
relative error = 1.2652494625793403000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7315797358756806 " "
Order of pole = 0.9433652842413736 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2499999999999836 " "
y[1] (analytic) = 1.3523092228228781 " "
y[1] (numeric) = 1.3523092228228952 " "
absolute error = 1.70974345792274100000000000000E-14 " "
relative error = 1.2643139816452165000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7332538386685323 " "
Order of pole = 0.9540177555735596 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2509999999999835 " "
y[1] (analytic) = 1.3533090735785702 " "
y[1] (numeric) = 1.3533090735785873 " "
absolute error = 1.70974345792274100000000000000E-14 " "
relative error = 1.263379882174031000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7349301764346508 " "
Order of pole = 0.9647066301911558 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2519999999999833 " "
y[1] (analytic) = 1.3543089252269809 " "
y[1] (numeric) = 1.354308925226998 " "
absolute error = 1.70974345792274100000000000000E-14 " "
relative error = 1.2624471611129565000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7366086785709338 " "
Order of pole = 0.9754309056705548 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2529999999999832 " "
y[1] (analytic) = 1.3553087777627706 " "
y[1] (numeric) = 1.3553087777627877 " "
absolute error = 1.70974345792274100000000000000E-14 " "
relative error = 1.2615158154181227000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7382892736205486 " "
Order of pole = 0.9861895655637554 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2539999999999831 " "
y[1] (analytic) = 1.3563086311806316 " "
y[1] (numeric) = 1.3563086311806491 " "
absolute error = 1.754152378907747300000000000000E-14 " "
relative error = 1.2933283314586025000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7399718892524276 " "
Order of pole = 0.9969815790277305 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.254999999999983 " "
y[1] (analytic) = 1.3573084854752888 " "
y[1] (numeric) = 1.3573084854753064 " "
absolute error = 1.754152378907747300000000000000E-14 " "
relative error = 1.292375607814384000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7416564522999396 " "
Order of pole = 1.0078059014013547 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.255999999999983 " "
y[1] (analytic) = 1.3583083406414975 " "
y[1] (numeric) = 1.3583083406415148 " "
absolute error = 1.731947918415244200000000000000E-14 " "
relative error = 1.275077143086146000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7433428887292803 " "
Order of pole = 1.018661473654804 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2569999999999828 " "
y[1] (analytic) = 1.3593081966740443 " "
y[1] (numeric) = 1.359308196674062 " "
absolute error = 1.776356839400250500000000000000E-14 " "
relative error = 1.3068094812836714000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7450311236747655 " "
Order of pole = 1.0295472229219413 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2579999999999827 " "
y[1] (analytic) = 1.360308053567748 " "
y[1] (numeric) = 1.360308053567766 " "
absolute error = 1.798561299892753600000000000000E-14 " "
relative error = 1.322172058877088000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7467210814141827 " "
Order of pole = 1.0404620620608416 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2589999999999826 " "
y[1] (analytic) = 1.3613079113174582 " "
y[1] (numeric) = 1.3613079113174762 " "
absolute error = 1.798561299892753600000000000000E-14 " "
relative error = 1.3212009457523290000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7484126854265418 " "
Order of pole = 1.0514048905110656 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2599999999999825 " "
y[1] (analytic) = 1.3623077699180548 " "
y[1] (numeric) = 1.3623077699180726 " "
absolute error = 1.776356839400250500000000000000E-14 " "
relative error = 1.3039321059639125000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7501058583426516 " "
Order of pole = 1.0623745934679363 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2609999999999824 " "
y[1] (analytic) = 1.3633076293644477 " "
y[1] (numeric) = 1.3633076293644657 " "
absolute error = 1.798561299892753600000000000000E-14 " "
relative error = 1.3192629903576600000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7518005220029484 " "
Order of pole = 1.073370042789282 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2619999999999822 " "
y[1] (analytic) = 1.3643074896515788 " "
y[1] (numeric) = 1.364307489651597 " "
absolute error = 1.820765760385256700000000000000E-14 " "
relative error = 1.3345714028516031000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.753496597428777 " "
Order of pole = 1.0843900964438369 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2629999999999821 " "
y[1] (analytic) = 1.3653073507744191 " "
y[1] (numeric) = 1.3653073507744375 " "
absolute error = 1.8429702208777599000000000000E-14 " "
relative error = 1.3498573927932084000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7551940048526096 " "
Order of pole = 1.0954335989998487 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.263999999999982 " "
y[1] (analytic) = 1.3663072127279705 " "
y[1] (numeric) = 1.366307212727989 " "
absolute error = 1.8429702208777599000000000000E-14 " "
relative error = 1.3488695687978428000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7568926637062015 " "
Order of pole = 1.1064993813695985 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.264999999999982 " "
y[1] (analytic) = 1.3673070755072638 " "
y[1] (numeric) = 1.3673070755072825 " "
absolute error = 1.86517468137026300000000000000E-14 " "
relative error = 1.3641227451984720000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7585924926467331 " "
Order of pole = 1.1175862611708158 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2659999999999818 " "
y[1] (analytic) = 1.3683069391073603 " "
y[1] (numeric) = 1.3683069391073792 " "
absolute error = 1.88737914186276600000000000000E-14 " "
relative error = 1.379353628867826800000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7602934095716847 " "
Order of pole = 1.128693042956609 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2669999999999817 " "
y[1] (analytic) = 1.3693068035233509 " "
y[1] (numeric) = 1.3693068035233698 " "
absolute error = 1.88737914186276600000000000000E-14 " "
relative error = 1.3783464282850040000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7619953316018024 " "
Order of pole = 1.1398185178662317 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2679999999999816 " "
y[1] (analytic) = 1.3703066687503556 " "
y[1] (numeric) = 1.370306668750374 " "
absolute error = 1.8429702208777599000000000000E-14 " "
relative error = 1.3449326803309272000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7636981751081928 " "
Order of pole = 1.1509614640270982 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2689999999999815 " "
y[1] (analytic) = 1.3713065347835223 " "
y[1] (numeric) = 1.3713065347835411 " "
absolute error = 1.88737914186276600000000000000E-14 " "
relative error = 1.3763364309794618000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7654018557085394 " "
Order of pole = 1.1621206465029612 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2699999999999814 " "
y[1] (analytic) = 1.3723064016180304 " "
y[1] (numeric) = 1.3723064016180495 " "
absolute error = 1.909583602355269200000000000000E-14 " "
relative error = 1.391514023474463800000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7671062883328659 " "
Order of pole = 1.1732948181515983 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2709999999999813 " "
y[1] (analytic) = 1.3733062692490867 " "
y[1] (numeric) = 1.3733062692491058 " "
absolute error = 1.909583602355269200000000000000E-14 " "
relative error = 1.3905008992636542000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7688113871283118 " "
Order of pole = 1.184482718274026 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2719999999999811 " "
y[1] (analytic) = 1.3743061376719268 " "
y[1] (numeric) = 1.3743061376719459 " "
absolute error = 1.909583602355269200000000000000E-14 " "
relative error = 1.3894892484363797000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7705170655737694 " "
Order of pole = 1.1956830742251743 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.272999999999981 " "
y[1] (analytic) = 1.3753060068818144 " "
y[1] (numeric) = 1.3753060068818337 " "
absolute error = 1.931788062847772400000000000000E-14 " "
relative error = 1.4046241732250200000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7722232364168093 " "
Order of pole = 1.2068946004906174 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.273999999999981 " "
y[1] (analytic) = 1.3763058768740422 " "
y[1] (numeric) = 1.376305876874062 " "
absolute error = 1.976196983832778600000000000000E-14 " "
relative error = 1.43587048274563000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7739298117381062 " "
Order of pole = 1.2181159996001032 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2749999999999808 " "
y[1] (analytic) = 1.3773057476439312 " "
y[1] (numeric) = 1.377305747643951 " "
absolute error = 1.976196983832778600000000000000E-14 " "
relative error = 1.43482809624031000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.775636702932559 " "
Order of pole = 1.2293459617822293 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2759999999999807 " "
y[1] (analytic) = 1.37830561918683 " "
y[1] (numeric) = 1.3783056191868495 " "
absolute error = 1.953992523340275500000000000000E-14 " "
relative error = 1.4176772525189937000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7773438207132288 " "
Order of pole = 1.2405831650855834 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2769999999999806 " "
y[1] (analytic) = 1.379305491498114 " "
y[1] (numeric) = 1.379305491498134 " "
absolute error = 1.998401444325281800000000000000E-14 " "
relative error = 1.4488461451384096000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7790510751423056 " "
Order of pole = 1.2518262757054686 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2779999999999805 " "
y[1] (analytic) = 1.3803053645731884 " "
y[1] (numeric) = 1.3803053645732084 " "
absolute error = 1.998401444325281800000000000000E-14 " "
relative error = 1.4477966221215247000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7807583758149885 " "
Order of pole = 1.2630739507997006 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2789999999999804 " "
y[1] (analytic) = 1.3813052384074844 " "
y[1] (numeric) = 1.3813052384075044 " "
absolute error = 1.998401444325281800000000000000E-14 " "
relative error = 1.4467486177271371000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7824656310471847 " "
Order of pole = 1.2743248262783986 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2799999999999803 " "
y[1] (analytic) = 1.382305112996461 " "
y[1] (numeric) = 1.382305112996481 " "
absolute error = 1.998401444325281800000000000000E-14 " "
relative error = 1.4457021286662913000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7841727494748685 " "
Order of pole = 1.2855775406965027 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2809999999999802 " "
y[1] (analytic) = 1.3833049883356037 " "
y[1] (numeric) = 1.3833049883356239 " "
absolute error = 2.02060590481778500000000000000E-14 " "
relative error = 1.4607088977890434000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7858796385205675 " "
Order of pole = 1.2968307122664875 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.28199999999998 " "
y[1] (analytic) = 1.384304864420426 " "
y[1] (numeric) = 1.3843048644204463 " "
absolute error = 2.02060590481778500000000000000E-14 " "
relative error = 1.4596538354748628000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7875862051448526 " "
Order of pole = 1.3080829501869182 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.28299999999998 " "
y[1] (analytic) = 1.385304741246468 " "
y[1] (numeric) = 1.3853047412464883 " "
absolute error = 2.02060590481778500000000000000E-14 " "
relative error = 1.4586002954120306000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7892923557367133 " "
Order of pole = 1.31933285282798 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2839999999999798 " "
y[1] (analytic) = 1.386304618809296 " "
y[1] (numeric) = 1.3863046188093162 " "
absolute error = 2.02060590481778500000000000000E-14 " "
relative error = 1.4575482743131113000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.790997996132686 " "
Order of pole = 1.3305790081146558 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2849999999999797 " "
y[1] (analytic) = 1.3873044971045025 " "
y[1] (numeric) = 1.387304497104523 " "
absolute error = 2.04281036531028800000000000000E-14 " "
relative error = 1.4725032388880144000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7927030315825785 " "
Order of pole = 1.3418199929371646 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2859999999999796 " "
y[1] (analytic) = 1.3883043761277074 " "
y[1] (numeric) = 1.3883043761277283 " "
absolute error = 2.087219286295294300000000000000E-14 " "
relative error = 1.5034306036814618000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7944073668534348 " "
Order of pole = 1.3530543747455592 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2869999999999795 " "
y[1] (analytic) = 1.389304255874557 " "
y[1] (numeric) = 1.3893042558745776 " "
absolute error = 2.065014825802791200000000000000E-14 " "
relative error = 1.4863661556286525000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.796110906124419 " "
Order of pole = 1.3642807098581695 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2879999999999794 " "
y[1] (analytic) = 1.390304136340722 " "
y[1] (numeric) = 1.3903041363407427 " "
absolute error = 2.065014825802791200000000000000E-14 " "
relative error = 1.4852971891732310000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7978135531181463 " "
Order of pole = 1.3754975455234284 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2889999999999793 " "
y[1] (analytic) = 1.3913040175218996 " "
y[1] (numeric) = 1.391304017521921 " "
absolute error = 2.131628207280300600000000000000E-14 " "
relative error = 1.5321081377145868000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.7995152110195165 " "
Order of pole = 1.3867034185927771 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2899999999999792 " "
y[1] (analytic) = 1.3923038994138144 " "
y[1] (numeric) = 1.3923038994138355 " "
absolute error = 2.109423746787797400000000000000E-14 " "
relative error = 1.5150598570297072000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8012157825683546 " "
Order of pole = 1.3978968569757395 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.290999999999979 " "
y[1] (analytic) = 1.3933037820122143 " "
y[1] (numeric) = 1.3933037820122356 " "
absolute error = 2.131628207280300600000000000000E-14 " "
relative error = 1.5299091517585600000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8029151700166661 " "
Order of pole = 1.4090763789507452 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.291999999999979 " "
y[1] (analytic) = 1.394303665312874 " "
y[1] (numeric) = 1.3943036653128953 " "
absolute error = 2.131628207280300600000000000000E-14 " "
relative error = 1.5288120230265442000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8046132751791897 " "
Order of pole = 1.420240493904343 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2929999999999788 " "
y[1] (analytic) = 1.3953035493115922 " "
y[1] (numeric) = 1.3953035493116142 " "
absolute error = 2.1982415887578100000000000000E-14 " "
relative error = 1.5754576055098313000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.806309999443223 " "
Order of pole = 1.4313877024448285 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2939999999999787 " "
y[1] (analytic) = 1.3963034340041949 " "
y[1] (numeric) = 1.3963034340042169 " "
absolute error = 2.1982415887578100000000000000E-14 " "
relative error = 1.5743294295666724000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8080052437620615 " "
Order of pole = 1.4425164963738144 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2949999999999786 " "
y[1] (analytic) = 1.3973033193865307 " "
y[1] (numeric) = 1.397303319386553 " "
absolute error = 2.22044604925031300000000000000E-14 " "
relative error = 1.589093805506146800000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8096989087028446 " "
Order of pole = 1.4536253592819115 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2959999999999785 " "
y[1] (analytic) = 1.3983032054544748 " "
y[1] (numeric) = 1.398303205454497 " "
absolute error = 2.22044604925031300000000000000E-14 " "
relative error = 1.5879574906135085000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8113908944527877 " "
Order of pole = 1.4647127667548752 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2969999999999784 " "
y[1] (analytic) = 1.3993030922039258 " "
y[1] (numeric) = 1.3993030922039482 " "
absolute error = 2.242650509742816200000000000000E-14 " "
relative error = 1.6026910268672415000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8130811008152132 " "
Order of pole = 1.4757771862154794 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2979999999999783 " "
y[1] (analytic) = 1.4003029796308073 " "
y[1] (numeric) = 1.40030297963083 " "
absolute error = 2.264854970235319300000000000000E-14 " "
relative error = 1.6174035213668210000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8147694272547914 " "
Order of pole = 1.4868170776356724 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.2989999999999782 " "
y[1] (analytic) = 1.4013028677310677 " "
y[1] (numeric) = 1.4013028677310906 " "
absolute error = 2.287059430727822500000000000000E-14 " "
relative error = 1.6320950191381078000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8164557729232371 " "
Order of pole = 1.4978308939005949 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.299999999999978 " "
y[1] (analytic) = 1.4023027565006796 " "
y[1] (numeric) = 1.4023027565007022 " "
absolute error = 2.264854970235319300000000000000E-14 " "
relative error = 1.61509699651954000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8181400366132836 " "
Order of pole = 1.5088170801578773 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.300999999999978 " "
y[1] (analytic) = 1.4033026459356381 " "
y[1] (numeric) = 1.4033026459356612 " "
absolute error = 2.309263891220325600000000000000E-14 " "
relative error = 1.6455922020161565000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.819822116866598 " "
Order of pole = 1.519774075357331 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3019999999999778 " "
y[1] (analytic) = 1.4043025360319648 " "
y[1] (numeric) = 1.4043025360319878 " "
absolute error = 2.309263891220325600000000000000E-14 " "
relative error = 1.6444205090916122000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.821501911918115 " "
Order of pole = 1.530700311497025 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3029999999999777 " "
y[1] (analytic) = 1.4053024267857033 " "
y[1] (numeric) = 1.4053024267857264 " "
absolute error = 2.309263891220325600000000000000E-14 " "
relative error = 1.6432504827464223000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8231793197748438 " "
Order of pole = 1.541594214755408 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3039999999999776 " "
y[1] (analytic) = 1.4063023181929213 " "
y[1] (numeric) = 1.4063023181929444 " "
absolute error = 2.309263891220325600000000000000E-14 " "
relative error = 1.6420821194319704000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.824854238173095 " "
Order of pole = 1.5524542048597123 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3049999999999775 " "
y[1] (analytic) = 1.40730221024971 " "
y[1] (numeric) = 1.4073022102497332 " "
absolute error = 2.331468351712828700000000000000E-14 " "
relative error = 1.656693448452082000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8265265646276059 " "
Order of pole = 1.563278695852457 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3059999999999774 " "
y[1] (analytic) = 1.4083021029521843 " "
y[1] (numeric) = 1.4083021029522076 " "
absolute error = 2.331468351712828700000000000000E-14 " "
relative error = 1.6555171982101260000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8281961964655677 " "
Order of pole = 1.5740660965639428 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3069999999999773 " "
y[1] (analytic) = 1.4093019962964821 " "
y[1] (numeric) = 1.4093019962965054 " "
absolute error = 2.331468351712828700000000000000E-14 " "
relative error = 1.6543426163020533000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8298630308536432 " "
Order of pole = 1.5848148110536897 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3079999999999772 " "
y[1] (analytic) = 1.4103018902787643 " "
y[1] (numeric) = 1.4103018902787878 " "
absolute error = 2.35367281220533200000000000000E-14 " "
relative error = 1.6689141725110349000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8315269647335395 " "
Order of pole = 1.595523237635378 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.308999999999977 " "
y[1] (analytic) = 1.4113017848952147 " "
y[1] (numeric) = 1.4113017848952385 " "
absolute error = 2.37587727269783500000000000000E-14 " "
relative error = 1.6834650803437035000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.83318789496808 " "
Order of pole = 1.6061897711205173 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.309999999999977 " "
y[1] (analytic) = 1.4123016801420405 " "
y[1] (numeric) = 1.4123016801420643 " "
absolute error = 2.37587727269783500000000000000E-14 " "
relative error = 1.6822732041633515000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8348457182482814 " "
Order of pole = 1.6168128013206218 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3109999999999769 " "
y[1] (analytic) = 1.4133015760154712 " "
y[1] (numeric) = 1.4133015760154948 " "
absolute error = 2.35367281220533200000000000000E-14 " "
relative error = 1.6653719575131687000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8365003311851322 " "
Order of pole = 1.6273907145829902 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3119999999999767 " "
y[1] (analytic) = 1.414301472511758 " "
y[1] (numeric) = 1.414301472511782 " "
absolute error = 2.39808173319033800000000000000E-14 " "
relative error = 1.6955944540815723000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8381516303153074 " "
Order of pole = 1.637921893780275 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3129999999999766 " "
y[1] (analytic) = 1.4153013696271761 " "
y[1] (numeric) = 1.4153013696272 " "
absolute error = 2.39808173319033800000000000000E-14 " "
relative error = 1.6943965325364232000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8397995120963277 " "
Order of pole = 1.6484047183043256 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3139999999999765 " "
y[1] (analytic) = 1.416301267358022 " "
y[1] (numeric) = 1.4163012673580462 " "
absolute error = 2.420286193682841300000000000000E-14 " "
relative error = 1.7088780822724672000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8414438729619595 " "
Order of pole = 1.658837564878267 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3149999999999764 " "
y[1] (analytic) = 1.4173011657006145 " "
y[1] (numeric) = 1.417301165700639 " "
absolute error = 2.442490654175344400000000000000E-14 " "
relative error = 1.7233391979664026000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8430846093072135 " "
Order of pole = 1.669218807373472 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3159999999999763 " "
y[1] (analytic) = 1.4183010646512948 " "
y[1] (numeric) = 1.4183010646513192 " "
absolute error = 2.442490654175344400000000000000E-14 " "
relative error = 1.722124247841454000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.844721617542753 " "
Order of pole = 1.679546817630177 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3169999999999762 " "
y[1] (analytic) = 1.419300964206425 " "
y[1] (numeric) = 1.4193009642064498 " "
absolute error = 2.486899575160350700000000000000E-14 " "
relative error = 1.7522002999207803000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8463547941175944 " "
Order of pole = 1.6898199658257127 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.317999999999976 " "
y[1] (analytic) = 1.42030086436239 " "
y[1] (numeric) = 1.420300864362415 " "
absolute error = 2.509104035652854000000000000000E-14 " "
relative error = 1.7666003722241314000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8479840354905142 " "
Order of pole = 1.7000366200505326 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.318999999999976 " "
y[1] (analytic) = 1.4213007651155964 " "
y[1] (numeric) = 1.4213007651156213 " "
absolute error = 2.486899575160350700000000000000E-14 " "
relative error = 1.7497349162111284000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8496092382153662 " "
Order of pole = 1.7101951476249084 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3199999999999759 " "
y[1] (analytic) = 1.4223006664624716 " "
y[1] (numeric) = 1.4223006664624962 " "
absolute error = 2.464695114667847500000000000000E-14 " "
relative error = 1.7328931728605643000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8512302989212097 " "
Order of pole = 1.7202939148233725 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3209999999999757 " "
y[1] (analytic) = 1.4233005683994637 " "
y[1] (numeric) = 1.4233005683994886 " "
absolute error = 2.486899575160350700000000000000E-14 " "
relative error = 1.747276457534848000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8528471143709542 " "
Order of pole = 1.73033128772148 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3219999999999756 " "
y[1] (analytic) = 1.4243004709230436 " "
y[1] (numeric) = 1.4243004709230689 " "
absolute error = 2.53130849614535700000000000000E-14 " "
relative error = 1.7772292769831752000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.854459581438884 " "
Order of pole = 1.7403056320164954 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3229999999999755 " "
y[1] (analytic) = 1.4253003740297032 " "
y[1] (numeric) = 1.4253003740297285 " "
absolute error = 2.53130849614535700000000000000E-14 " "
relative error = 1.7759824821968400000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8560675971659872 " "
Order of pole = 1.7502153137150387 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3239999999999754 " "
y[1] (analytic) = 1.4263002777159546 " "
y[1] (numeric) = 1.4263002777159797 " "
absolute error = 2.509104035652854000000000000000E-14 " "
relative error = 1.759169562576877000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8576710587771428 " "
Order of pole = 1.7600586995709566 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3249999999999753 " "
y[1] (analytic) = 1.42730018197833 " "
y[1] (numeric) = 1.4273001819783555 " "
absolute error = 2.5535129566378600000000000000E-14 " "
relative error = 1.7890510972250606000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8592698636979579 " "
Order of pole = 1.7698341572328218 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3259999999999752 " "
y[1] (analytic) = 1.4283000868133842 " "
y[1] (numeric) = 1.42830008681341 " "
absolute error = 2.575717417130363000000000000000E-14 " "
relative error = 1.8033447178995346000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8608639095794486 " "
Order of pole = 1.7795400557646932 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.326999999999975 " "
y[1] (analytic) = 1.429299992217692 " "
y[1] (numeric) = 1.4292999922177176 " "
absolute error = 2.5535129566378600000000000000E-14 " "
relative error = 1.7865479399295642000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.862453094320917 " "
Order of pole = 1.789174765930042 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.327999999999975 " "
y[1] (analytic) = 1.430299898187848 " "
y[1] (numeric) = 1.4302998981878736 " "
absolute error = 2.5535129566378600000000000000E-14 " "
relative error = 1.785298985110111000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8640373161240873 " "
Order of pole = 1.7987366610716435 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3289999999999749 " "
y[1] (analytic) = 1.4312998047204673 " "
y[1] (numeric) = 1.4312998047204932 " "
absolute error = 2.597921877622866300000000000000E-14 " "
relative error = 1.8150787620139724000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8656164734388216 " "
Order of pole = 1.80822411634799 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3299999999999748 " "
y[1] (analytic) = 1.4322997118121856 " "
y[1] (numeric) = 1.432299711812212 " "
absolute error = 2.642330798607872600000000000000E-14 " "
relative error = 1.844816958920365800000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8671904650729427 " "
Order of pole = 1.8176355103705326 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3309999999999746 " "
y[1] (analytic) = 1.4332996194596594 " "
y[1] (numeric) = 1.4332996194596859 " "
absolute error = 2.642330798607872600000000000000E-14 " "
relative error = 1.8435299659145982000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8687591901696006 " "
Order of pole = 1.8269692249267386 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3319999999999745 " "
y[1] (analytic) = 1.434299527659564 " "
y[1] (numeric) = 1.4342995276595905 " "
absolute error = 2.642330798607872600000000000000E-14 " "
relative error = 1.8422447666280198000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8703225482250874 " "
Order of pole = 1.836223645290314 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3329999999999744 " "
y[1] (analytic) = 1.435299436408595 " "
y[1] (numeric) = 1.4352994364086213 " "
absolute error = 2.642330798607872600000000000000E-14 " "
relative error = 1.8409613573175437000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8718804391621495 " "
Order of pole = 1.8453971613435982 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3339999999999743 " "
y[1] (analytic) = 1.4362993457034676 " "
y[1] (numeric) = 1.436299345703494 " "
absolute error = 2.642330798607872600000000000000E-14 " "
relative error = 1.8396797342504653000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8734327632951082 " "
Order of pole = 1.8544881670843445 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3349999999999742 " "
y[1] (analytic) = 1.437299255540917 " "
y[1] (numeric) = 1.4372992555409436 " "
absolute error = 2.664535259100375700000000000000E-14 " "
relative error = 1.853848632306984200000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.874979421385681 " "
Order of pole = 1.8634950615216148 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.335999999999974 " "
y[1] (analytic) = 1.4382991659176982 " "
y[1] (numeric) = 1.4382991659177247 " "
absolute error = 2.642330798607872600000000000000E-14 " "
relative error = 1.8371218319673774000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8765203146917924 " "
Order of pole = 1.8724162494406524 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.336999999999974 " "
y[1] (analytic) = 1.4392990768305844 " "
y[1] (numeric) = 1.4392990768306113 " "
absolute error = 2.686739719592879000000000000000E-14 " "
relative error = 1.86670009231801020000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8780553449106014 " "
Order of pole = 1.8812501405935294 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3379999999999739 " "
y[1] (analytic) = 1.4402989882763697 " "
y[1] (numeric) = 1.4402989882763966 " "
absolute error = 2.686739719592879000000000000000E-14 " "
relative error = 1.8654041566800975000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8795844142992684 " "
Order of pole = 1.8899951515533449 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3389999999999738 " "
y[1] (analytic) = 1.441298900251866 " "
y[1] (numeric) = 1.441298900251893 " "
absolute error = 2.70894418008538200000000000000E-14 " "
relative error = 1.8795158864077370000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8811074255998388 " "
Order of pole = 1.8986497046070756 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3399999999999737 " "
y[1] (analytic) = 1.4422988127539051 " "
y[1] (numeric) = 1.4422988127539322 " "
absolute error = 2.70894418008538200000000000000E-14 " "
relative error = 1.878212861392406000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.882624282220699 " "
Order of pole = 1.9072122305252694 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3409999999999735 " "
y[1] (analytic) = 1.4432987257793373 " "
y[1] (numeric) = 1.4432987257793646 " "
absolute error = 2.73114864057788500000000000000E-14 " "
relative error = 1.8922961628079787000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8841348880560118 " "
Order of pole = 1.9156811659620718 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3419999999999734 " "
y[1] (analytic) = 1.444298639325032 " "
y[1] (numeric) = 1.4442986393250594 " "
absolute error = 2.73114864057788500000000000000E-14 " "
relative error = 1.8909860926367972000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8856391476579568 " "
Order of pole = 1.9240549559986473 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3429999999999733 " "
y[1] (analytic) = 1.4452985533878773 " "
y[1] (numeric) = 1.4452985533879048 " "
absolute error = 2.753353101070388000000000000000E-14 " "
relative error = 1.9050410689309435000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8871369662358033 " "
Order of pole = 1.9323320542198275 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3439999999999732 " "
y[1] (analytic) = 1.4462984679647797 " "
y[1] (numeric) = 1.4462984679648074 " "
absolute error = 2.775557561562891400000000000000E-14 " "
relative error = 1.919076610423735800000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8886282496308897 " "
Order of pole = 1.9405109223528125 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.344999999999973 " "
y[1] (analytic) = 1.4472983830526644 " "
y[1] (numeric) = 1.4472983830526924 " "
absolute error = 2.797762022055394500000000000000E-14 " "
relative error = 1.9330927573858755000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8901129044048757 " "
Order of pole = 1.9485900317045726 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.345999999999973 " "
y[1] (analytic) = 1.4482982986484754 " "
y[1] (numeric) = 1.4482982986485031 " "
absolute error = 2.775557561562891400000000000000E-14 " "
relative error = 1.9164267224182957000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.891590837791488 " "
Order of pole = 1.956567862382311 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3469999999999729 " "
y[1] (analytic) = 1.4492982147491738 " "
y[1] (numeric) = 1.4492982147492017 " "
absolute error = 2.797762022055394500000000000000E-14 " "
relative error = 1.93042535592966000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8930619578408483 " "
Order of pole = 1.9644429056350745 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3479999999999728 " "
y[1] (analytic) = 1.45029813135174 " "
y[1] (numeric) = 1.450298131351768 " "
absolute error = 2.797762022055394500000000000000E-14 " "
relative error = 1.9290944127796403000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8945261732721212 " "
Order of pole = 1.972213661587249 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3489999999999727 " "
y[1] (analytic) = 1.4512980484531717 " "
y[1] (numeric) = 1.4512980484532 " "
absolute error = 2.84217094304040100000000000000E-14 " "
relative error = 1.958364752209000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8959833936562414 " "
Order of pole = 1.9798786420321974 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3499999999999726 " "
y[1] (analytic) = 1.4522979660504856 " "
y[1] (numeric) = 1.452297966050514 " "
absolute error = 2.84217094304040100000000000000E-14 " "
relative error = 1.9570164039888216000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.897433529382746 " "
Order of pole = 1.9874363701072824 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3509999999999724 " "
y[1] (analytic) = 1.4532978841407145 " "
y[1] (numeric) = 1.4532978841407433 " "
absolute error = 2.88657986402540700000000000000E-14 " "
relative error = 1.9862272528746874000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.8988764916390668 " "
Order of pole = 1.9948853798753259 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3519999999999723 " "
y[1] (analytic) = 1.4542978027209112 " "
y[1] (numeric) = 1.4542978027209401 " "
absolute error = 2.88657986402540700000000000000E-14 " "
relative error = 1.9848616003027544000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9003121925365247 " "
Order of pole = 2.002224218331655 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3529999999999722 " "
y[1] (analytic) = 1.455297721788144 " "
y[1] (numeric) = 1.4552977217881733 " "
absolute error = 2.93098878501041300000000000000E-14 " "
relative error = 2.014013174849932000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9017405449737954 " "
Order of pole = 2.0094514434622504 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3539999999999721 " "
y[1] (analytic) = 1.4562976413395006 " "
y[1] (numeric) = 1.4562976413395297 " "
absolute error = 2.9087843245179100000000000000E-14 " "
relative error = 1.9973831186339178000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9031614628898212 " "
Order of pole = 2.016565627987042 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.354999999999972 " "
y[1] (analytic) = 1.4572975613720842 " "
y[1] (numeric) = 1.4572975613721135 " "
absolute error = 2.93098878501041300000000000000E-14 " "
relative error = 2.011249358196146000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9045748610994306 " "
Order of pole = 2.0235653570198124 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.355999999999972 " "
y[1] (analytic) = 1.458297481883017 " "
y[1] (numeric) = 1.4582974818830463 " "
absolute error = 2.93098878501041300000000000000E-14 " "
relative error = 2.0098702915030706000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9059806553609155 " "
Order of pole = 2.030449229084816 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3569999999999718 " "
y[1] (analytic) = 1.459297402869437 " "
y[1] (numeric) = 1.4592974028694667 " "
absolute error = 2.975397705995419500000000000000E-14 " "
relative error = 2.0389248278965297000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9073787624837164 " "
Order of pole = 2.0372158578602395 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3579999999999717 " "
y[1] (analytic) = 1.4602973243285013 " "
y[1] (numeric) = 1.4602973243285307 " "
absolute error = 2.93098878501041300000000000000E-14 " "
relative error = 2.0071178219532726000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9087691002046954 " "
Order of pole = 2.0438638702685665 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3589999999999716 " "
y[1] (analytic) = 1.461297246257381 " "
y[1] (numeric) = 1.461297246257411 " "
absolute error = 2.997602166487922700000000000000E-14 " "
relative error = 2.051329511613921900000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.91015158736094 " "
Order of pole = 2.0503919091927294 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3599999999999715 " "
y[1] (analytic) = 1.4622971686532678 " "
y[1] (numeric) = 1.4622971686532975 " "
absolute error = 2.975397705995419500000000000000E-14 " "
relative error = 2.0347421644368446000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9115261438571538 " "
Order of pole = 2.056798633023611 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3609999999999713 " "
y[1] (analytic) = 1.4632970915133665 " "
y[1] (numeric) = 1.4632970915133967 " "
absolute error = 3.01980662698042600000000000000E-14 " "
relative error = 2.063700286492943800000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9128926906421557 " "
Order of pole = 2.063082715361457 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3619999999999712 " "
y[1] (analytic) = 1.4642970148349017 " "
y[1] (numeric) = 1.464297014834932 " "
absolute error = 3.01980662698042600000000000000E-14 " "
relative error = 2.062291049142722000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.91425114978881 " "
Order of pole = 2.0692428462845562 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3629999999999711 " "
y[1] (analytic) = 1.4652969386151122 " "
y[1] (numeric) = 1.4652969386151429 " "
absolute error = 3.06421554796543200000000000000E-14 " "
relative error = 2.0911908482259556000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9156014445063683 " "
Order of pole = 2.07527773247611 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.363999999999971 " "
y[1] (analytic) = 1.4662968628512558 " "
y[1] (numeric) = 1.466296862851286 " "
absolute error = 3.01980662698042600000000000000E-14 " "
relative error = 2.0594783385871304000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9169434991912966 " "
Order of pole = 2.081186098240206 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.364999999999971 " "
y[1] (analytic) = 1.4672967875406036 " "
y[1] (numeric) = 1.4672967875406342 " "
absolute error = 3.06421554796543200000000000000E-14 " "
relative error = 2.08834066426431000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.918277239361803 " "
Order of pole = 2.0869666843701147 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3659999999999708 " "
y[1] (analytic) = 1.468296712680446 " "
y[1] (numeric) = 1.4682967126804767 " "
absolute error = 3.06421554796543200000000000000E-14 " "
relative error = 2.0869184828259676000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9196025917767512 " "
Order of pole = 2.0926182501119257 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3669999999999707 " "
y[1] (analytic) = 1.4692966382680879 " "
y[1] (numeric) = 1.4692966382681185 " "
absolute error = 3.06421554796543200000000000000E-14 " "
relative error = 2.0854982364741076000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9209194844076902 " "
Order of pole = 2.09813957274514 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3679999999999706 " "
y[1] (analytic) = 1.4702965643008505 " "
y[1] (numeric) = 1.4702965643008814 " "
absolute error = 3.08642000845793500000000000000E-14 " "
relative error = 2.099181949680728000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9222278465078155 " "
Order of pole = 2.103529448685336 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3689999999999705 " "
y[1] (analytic) = 1.4712964907760717 " "
y[1] (numeric) = 1.4712964907761026 " "
absolute error = 3.08642000845793500000000000000E-14 " "
relative error = 2.097755298002462200000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9235276085470392 " "
Order of pole = 2.1087866925276053 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3699999999999704 " "
y[1] (analytic) = 1.4722964176911046 " "
y[1] (numeric) = 1.4722964176911357 " "
absolute error = 3.10862446895043830000000000000E-14 " "
relative error = 2.11141209854023000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9248187023455283 " "
Order of pole = 2.1139101391335977 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3709999999999702 " "
y[1] (analytic) = 1.4732963450433185 " "
y[1] (numeric) = 1.4732963450433498 " "
absolute error = 3.130828929442941400000000000000E-14 " "
relative error = 2.1250503606936508000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9261010610016303 " "
Order of pole = 2.118898642596573 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3719999999999701 " "
y[1] (analytic) = 1.4742962728300986 " "
y[1] (numeric) = 1.4742962728301299 " "
absolute error = 3.130828929442941400000000000000E-14 " "
relative error = 2.1236090649764164000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9273746189924466 " "
Order of pole = 2.1237510777788238 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.37299999999997 " "
y[1] (analytic) = 1.475296201048845 " "
y[1] (numeric) = 1.4752962010488764 " "
absolute error = 3.130828929442941400000000000000E-14 " "
relative error = 2.1221697224036193000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9286393120950456 " "
Order of pole = 2.128466339164042 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.37399999999997 " "
y[1] (analytic) = 1.4762961296969743 " "
y[1] (numeric) = 1.4762961296970056 " "
absolute error = 3.130828929442941400000000000000E-14 " "
relative error = 2.120732329011509000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9298950775505068 " "
Order of pole = 2.1330433434559666 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3749999999999698 " "
y[1] (analytic) = 1.4772960587719177 " "
y[1] (numeric) = 1.4772960587719492 " "
absolute error = 3.153033389935444600000000000000E-14 " "
relative error = 2.1343273551792824000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9311418539193466 " "
Order of pole = 2.137481027332221 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3759999999999697 " "
y[1] (analytic) = 1.4782959882711226 " "
y[1] (numeric) = 1.4782959882711542 " "
absolute error = 3.153033389935444600000000000000E-14 " "
relative error = 2.1328836815845917000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9323795812595783 " "
Order of pole = 2.1417783502316645 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3769999999999696 " "
y[1] (analytic) = 1.4792959181920509 " "
y[1] (numeric) = 1.4792959181920826 " "
absolute error = 3.17523785042794770000000000000E-14 " "
relative error = 2.1464521137249024000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.933608201044439 " "
Order of pole = 2.145934293149974 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3779999999999695 " "
y[1] (analytic) = 1.4802958485321802 " "
y[1] (numeric) = 1.480295848532212 " "
absolute error = 3.17523785042794770000000000000E-14 " "
relative error = 2.1450021991052834000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9348276561898266 " "
Order of pole = 2.1499478590999317 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3789999999999694 " "
y[1] (analytic) = 1.4812957792890027 " "
y[1] (numeric) = 1.4812957792890347 " "
absolute error = 3.19744231092045100000000000000E-14 " "
relative error = 2.1585441311762665000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9360378911416687 " "
Order of pole = 2.153818074451383 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3799999999999693 " "
y[1] (analytic) = 1.4822957104600265 " "
y[1] (numeric) = 1.4822957104600585 " "
absolute error = 3.19744231092045100000000000000E-14 " "
relative error = 2.157088014461118200000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9372388518347803 " "
Order of pole = 2.157543988390074 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3809999999999691 " "
y[1] (analytic) = 1.4832956420427734 " "
y[1] (numeric) = 1.4832956420428056 " "
absolute error = 3.21964677141295400000000000000E-14 " "
relative error = 2.1706035399516868000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.938430485742844 " "
Order of pole = 2.161124673657284 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.381999999999969 " "
y[1] (analytic) = 1.484295574034781 " "
y[1] (numeric) = 1.4842955740348132 " "
absolute error = 3.21964677141295400000000000000E-14 " "
relative error = 2.1691412598239743000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9396127418435285 " "
Order of pole = 2.1645592260968556 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.382999999999969 " "
y[1] (analytic) = 1.4852955064336006 " "
y[1] (numeric) = 1.4852955064336333 " "
absolute error = 3.2640556923979600000000000000E-14 " "
relative error = 2.197579995535978200000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.9407855707351795 " "
Order of pole = 2.1678467664256615 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3839999999999688 " "
y[1] (analytic) = 1.4862954392368 " "
y[1] (numeric) = 1.4862954392368324 " "
absolute error = 3.24185123190545700000000000000E-14 " "
relative error = 2.1811620666548773000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.941948924548343 " "
Order of pole = 2.1709864390063487 " "
" "
"TOP MAIN SOLVE Loop"
x[1] = 1.3849999999999687 " "
y[1] (analytic) = 1.4872953724419589 " "
y[1] (numeric) = 1.4872953724419917 " "
absolute error = 3.286260152890463400000000000000E-14 " "
relative error = 2.209554479749925000000000000E-12 "%"
Correct digits = 14
h = 1.000E-3 " "
"Complex estimate of poles used"
Radius of convergence = 1.94310275704797 " "
Order of pole = 2.173977413305604 " "
"Finished!"
"Maximum Time Reached before Solution Completed!"
"diff ( y , x , 1 ) = tanh (3.0 * x + 1.0 ) ;"
Iterations = 286
"Total Elapsed Time "= 0 Years 0 Days 0 Hours 3 Minutes 0 Seconds
"Elapsed Time(since restart) "= 0 Years 0 Days 0 Hours 2 Minutes 59 Seconds
"Expected Time Remaining "= 0 Years 0 Days 0 Hours 6 Minutes 25 Seconds
"Optimized Time Remaining "= 0 Years 0 Days 0 Hours 6 Minutes 23 Seconds
"Expected Total Time "= 0 Years 0 Days 0 Hours 9 Minutes 24 Seconds
"Time to Timeout " Unknown
Percent Done = 31.888888888885383 "%"
(%o57) true
(%o57) diffeq.max