(%i1) batch(diffeq.max) read and interpret file: /home/dennis/mastersource/mine/omnisode/diffeq.max (%i2) load(stringproc) (%o2) /usr/local/share/maxima/5.26.0/share/contrib/stringproc/stringproc.mac (%i3) display_alot(iter) := if iter >= 0 then (ind_var : array_x , omniout_float(ALWAYS, 1 "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 : abs(numeric_val - analytic_val_y), omniout_float(ALWAYS, "y[1] (numeric) ", 33, numeric_val, abserr 100.0 20, " "), if abs(analytic_val_y) # 0.0 then relerr : ------------------- abs(analytic_val_y) else relerr : - 1.0, 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_float(ALWAYS, "h ", 4, glob_h, 20, " ")) (%o3) display_alot(iter) := if iter >= 0 then (ind_var : array_x , omniout_float(ALWAYS, 1 "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 : abs(numeric_val - analytic_val_y), omniout_float(ALWAYS, "y[1] (numeric) ", 33, numeric_val, abserr 100.0 20, " "), if abs(analytic_val_y) # 0.0 then relerr : ------------------- abs(analytic_val_y) else relerr : - 1.0, 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_float(ALWAYS, "h ", 4, glob_h, 20, " ")) (%i4) adjust_for_pole(h_param) := block(hnew : h_param, glob_normmax : glob_small_float, if !array_y_higher ! > glob_small_float ! 1, 1! then (tmp : !array_y_higher !, if tmp < glob_normmax ! 1, 1! then glob_normmax : tmp), if glob_look_poles and (!array_pole ! > glob_small_float) and (array_pole # glob_large_float) ! 1! 1 array_pole 1 then (sz2 : -----------, if sz2 < hnew 10.0 then (omniout_float(INFO, "glob_h adjusted to ", 20, h_param, 12, "due to singularity."), omniout_str(INFO, "Reached Optimal"), newline(), 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) 1 (%o4) adjust_for_pole(h_param) := block(hnew : h_param, glob_normmax : glob_small_float, if !array_y_higher ! > glob_small_float ! 1, 1! then (tmp : !array_y_higher !, if tmp < glob_normmax ! 1, 1! then glob_normmax : tmp), if glob_look_poles and (!array_pole ! > glob_small_float) and (array_pole # glob_large_float) ! 1! 1 array_pole 1 then (sz2 : -----------, if sz2 < hnew 10.0 then (omniout_float(INFO, "glob_h adjusted to ", 20, h_param, 12, "due to singularity."), omniout_str(INFO, "Reached Optimal"), newline(), 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) 1 (%i5) prog_report(x_start, x_end) := (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)), 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, "Time to Timeout "), omniout_timestr(convfloat(left_sec)), omniout_float(INFO, "Percent Done ", 33, percent_done, 4, "%")) (%o5) prog_report(x_start, x_end) := (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)), 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, "Time to Timeout "), omniout_timestr(convfloat(left_sec)), omniout_float(INFO, "Percent Done ", 33, percent_done, 4, "%")) (%i6) check_for_pole() := (n : glob_max_terms, m : - 1 - 1 + n, while (m >= 10) and ((!array_y_higher ! < glob_small_float) ! 1, m! or (!array_y_higher ! < glob_small_float) ! 1, m - 1! or (!array_y_higher ! < glob_small_float)) do m : m - 1, ! 1, m - 2! array_y_higher array_y_higher 1, m 1, m - 1 if m > 10 then (rm0 : ----------------------, rm1 : ----------------------, array_y_higher array_y_higher 1, m - 1 1, m - 2 hdrc : convfloat(m - 1) rm0 - convfloat(m - 2) rm1, glob_h if abs(hdrc) > glob_small_float then (rcs : ------, hdrc convfloat(m - 1) rm0 ord_no : 2.0 - convfloat(m) + --------------------, array_real_pole : rcs, hdrc 1, 1 array_real_pole : ord_no) else (array_real_pole : glob_large_float, 1, 2 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 !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 (array_complex_pole : glob_large_float, 1, 1 array_complex_pole : glob_large_float) 1, 2 elseif (!array_y_higher ! >= glob_large_float) ! 1, m! or (!array_y_higher ! >= glob_large_float) ! 1, m - 1! or (!array_y_higher ! >= glob_large_float) ! 1, m - 2! or (!array_y_higher ! >= glob_large_float) ! 1, m - 3! or (!array_y_higher ! >= glob_large_float) ! 1, m - 4! or (!array_y_higher ! >= glob_large_float) ! 1, m - 5! then (array_complex_pole : glob_large_float, 1, 1 array_complex_pole : glob_large_float) 1, 2 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 (abs(nr1 dr2 - nr2 dr1) <= glob_small_float) rm4 rm3 rm2 or (abs(dr1) <= glob_small_float) then (array_complex_pole : 1, 1 glob_large_float, array_complex_pole : glob_large_float) 1, 2 else (if abs(nr1 dr2 - nr2 dr1) > glob_small_float dr1 dr2 - ds2 dr1 + ds1 dr2 then (rcs : ---------------------------, nr1 dr2 - nr2 dr1 rcs nr1 - ds1 convfloat(m) ord_no : ------------- - ------------, 2.0 dr1 2.0 if abs(rcs) > glob_small_float then (if rcs > 0.0 then rad_c : sqrt(rcs) 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, 1, 1 array_complex_pole : ord_no), found : false, 1, 2 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 omniout_str(ALWAYS, "Complex estimate of poles used")), if (not found) and ((array_real_pole # glob_large_float) and (array_real_pole # glob_large_float) 1, 1 1, 2 and (array_real_pole > 0.0) and (array_real_pole > 0.0) 1, 1 1, 2 and ((array_complex_pole = glob_large_float) or (array_complex_pole = glob_large_float) or (array_complex_pole <= 0.0) or (array_complex_pole <= 0.0))) 1, 1 1, 2 1, 1 1, 2 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 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 glob_display_flag 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 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 omniout_str(ALWAYS, "Complex estimate of poles used")), if not found then (array_poles : glob_large_float, array_poles : glob_large_float, 1, 1 1, 2 array_type_pole : 3, if glob_display_flag 1 then omniout_str(ALWAYS, "NO POLE")), array_pole : glob_large_float, 1 array_pole : glob_large_float, if array_pole > array_poles 2 1 1, 1 then (array_pole : array_poles , array_pole : array_poles ), 1 1, 1 2 1, 2 display_pole()) (%o6) check_for_pole() := (n : glob_max_terms, m : - 1 - 1 + n, while (m >= 10) and ((!array_y_higher ! < glob_small_float) ! 1, m! or (!array_y_higher ! < glob_small_float) ! 1, m - 1! or (!array_y_higher ! < glob_small_float)) do m : m - 1, ! 1, m - 2! array_y_higher array_y_higher 1, m 1, m - 1 if m > 10 then (rm0 : ----------------------, rm1 : ----------------------, array_y_higher array_y_higher 1, m - 1 1, m - 2 hdrc : convfloat(m - 1) rm0 - convfloat(m - 2) rm1, glob_h if abs(hdrc) > glob_small_float then (rcs : ------, hdrc convfloat(m - 1) rm0 ord_no : 2.0 - convfloat(m) + --------------------, array_real_pole : rcs, hdrc 1, 1 array_real_pole : ord_no) else (array_real_pole : glob_large_float, 1, 2 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 !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 (array_complex_pole : glob_large_float, 1, 1 array_complex_pole : glob_large_float) 1, 2 elseif (!array_y_higher ! >= glob_large_float) ! 1, m! or (!array_y_higher ! >= glob_large_float) ! 1, m - 1! or (!array_y_higher ! >= glob_large_float) ! 1, m - 2! or (!array_y_higher ! >= glob_large_float) ! 1, m - 3! or (!array_y_higher ! >= glob_large_float) ! 1, m - 4! or (!array_y_higher ! >= glob_large_float) ! 1, m - 5! then (array_complex_pole : glob_large_float, 1, 1 array_complex_pole : glob_large_float) 1, 2 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 (abs(nr1 dr2 - nr2 dr1) <= glob_small_float) rm4 rm3 rm2 or (abs(dr1) <= glob_small_float) then (array_complex_pole : 1, 1 glob_large_float, array_complex_pole : glob_large_float) 1, 2 else (if abs(nr1 dr2 - nr2 dr1) > glob_small_float dr1 dr2 - ds2 dr1 + ds1 dr2 then (rcs : ---------------------------, nr1 dr2 - nr2 dr1 rcs nr1 - ds1 convfloat(m) ord_no : ------------- - ------------, 2.0 dr1 2.0 if abs(rcs) > glob_small_float then (if rcs > 0.0 then rad_c : sqrt(rcs) 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, 1, 1 array_complex_pole : ord_no), found : false, 1, 2 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 omniout_str(ALWAYS, "Complex estimate of poles used")), if (not found) and ((array_real_pole # glob_large_float) and (array_real_pole # glob_large_float) 1, 1 1, 2 and (array_real_pole > 0.0) and (array_real_pole > 0.0) 1, 1 1, 2 and ((array_complex_pole = glob_large_float) or (array_complex_pole = glob_large_float) or (array_complex_pole <= 0.0) or (array_complex_pole <= 0.0))) 1, 1 1, 2 1, 1 1, 2 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 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 glob_display_flag 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 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 omniout_str(ALWAYS, "Complex estimate of poles used")), if not found then (array_poles : glob_large_float, array_poles : glob_large_float, 1, 1 1, 2 array_type_pole : 3, if glob_display_flag 1 then omniout_str(ALWAYS, "NO POLE")), array_pole : glob_large_float, 1 array_pole : glob_large_float, if array_pole > array_poles 2 1 1, 1 then (array_pole : array_poles , array_pole : array_poles ), 1 1, 1 2 1, 2 display_pole()) (%i7) get_norms() := if not glob_initial_pass then (set_z(array_norms, 1 + glob_max_terms), iii : 1, while iii <= glob_max_terms do (if !array_y ! > array_norms ! iii! iii then array_norms : !array_y !, iii : 1 + iii)) iii ! iii! (%o7) get_norms() := if not glob_initial_pass then (set_z(array_norms, 1 + glob_max_terms), iii : 1, while iii <= glob_max_terms do (if !array_y ! > array_norms ! iii! iii then array_norms : !array_y !, iii : 1 + iii)) iii ! iii! (%i8) atomall() := (array_tmp1 : array_y array_y , 1 1 1 array_tmp2 : array_tmp1 + array_const_0D0 , 1 1 1 if not array_y_set_initial then (if 1 <= glob_max_terms 1, 2 1 then (temporary : array_tmp2 glob_h factorial_3(0, 1), 1 array_y : temporary, array_y_higher : temporary, 2 1, 2 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 2, glob_h 2, 1 array_tmp1 : ats(2, array_y, array_y, 1), 2 array_tmp2 : array_tmp1 + array_const_0D0 , 2 2 2 if not array_y_set_initial then (if 2 <= glob_max_terms 1, 3 1 then (temporary : array_tmp2 glob_h factorial_3(1, 2), 2 array_y : temporary, array_y_higher : temporary, 3 1, 3 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 3, glob_h 2, 2 array_tmp1 : ats(3, array_y, array_y, 1), 3 array_tmp2 : array_tmp1 + array_const_0D0 , 3 3 3 if not array_y_set_initial then (if 3 <= glob_max_terms 1, 4 1 then (temporary : array_tmp2 glob_h factorial_3(2, 3), 3 array_y : temporary, array_y_higher : temporary, 4 1, 4 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 4, glob_h 2, 3 array_tmp1 : ats(4, array_y, array_y, 1), 4 array_tmp2 : array_tmp1 + array_const_0D0 , 4 4 4 if not array_y_set_initial then (if 4 <= glob_max_terms 1, 5 1 then (temporary : array_tmp2 glob_h factorial_3(3, 4), 4 array_y : temporary, array_y_higher : temporary, 5 1, 5 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 5, glob_h 2, 4 array_tmp1 : ats(5, array_y, array_y, 1), 5 array_tmp2 : array_tmp1 + array_const_0D0 , 5 5 5 if not array_y_set_initial then (if 5 <= glob_max_terms 1, 6 1 then (temporary : array_tmp2 glob_h factorial_3(4, 5), 5 array_y : temporary, array_y_higher : temporary, 6 1, 6 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 6, glob_h 2, 5 while kkk <= glob_max_terms do (array_tmp1 : kkk ats(kkk, array_y, array_y, 1), array_tmp2 : kkk array_tmp1 + array_const_0D0 , order_d : 1, kkk kkk if 1 + order_d + kkk <= glob_max_terms then (if not array_y_set_initial 1, order_d + kkk order_d array_tmp2 glob_h kkk then (temporary : -----------------------------------------, factorial_3(kkk - 1, - 1 + order_d + kkk) array_y : temporary, array_y_higher : temporary, order_d + kkk 1, order_d + kkk term : - 1 + order_d + kkk, adj2 : 2, while (adj2 <= 1 + order_d) temporary convfp(adj2) and (term >= 1) do (temporary : ----------------------, glob_h array_y_higher : temporary, adj2 : 1 + adj2, term : term - 1))), adj2, term kkk : 1 + kkk)) (%o8) atomall() := (array_tmp1 : array_y array_y , 1 1 1 array_tmp2 : array_tmp1 + array_const_0D0 , 1 1 1 if not array_y_set_initial then (if 1 <= glob_max_terms 1, 2 1 then (temporary : array_tmp2 glob_h factorial_3(0, 1), 1 array_y : temporary, array_y_higher : temporary, 2 1, 2 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 2, glob_h 2, 1 array_tmp1 : ats(2, array_y, array_y, 1), 2 array_tmp2 : array_tmp1 + array_const_0D0 , 2 2 2 if not array_y_set_initial then (if 2 <= glob_max_terms 1, 3 1 then (temporary : array_tmp2 glob_h factorial_3(1, 2), 2 array_y : temporary, array_y_higher : temporary, 3 1, 3 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 3, glob_h 2, 2 array_tmp1 : ats(3, array_y, array_y, 1), 3 array_tmp2 : array_tmp1 + array_const_0D0 , 3 3 3 if not array_y_set_initial then (if 3 <= glob_max_terms 1, 4 1 then (temporary : array_tmp2 glob_h factorial_3(2, 3), 3 array_y : temporary, array_y_higher : temporary, 4 1, 4 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 4, glob_h 2, 3 array_tmp1 : ats(4, array_y, array_y, 1), 4 array_tmp2 : array_tmp1 + array_const_0D0 , 4 4 4 if not array_y_set_initial then (if 4 <= glob_max_terms 1, 5 1 then (temporary : array_tmp2 glob_h factorial_3(3, 4), 4 array_y : temporary, array_y_higher : temporary, 5 1, 5 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 5, glob_h 2, 4 array_tmp1 : ats(5, array_y, array_y, 1), 5 array_tmp2 : array_tmp1 + array_const_0D0 , 5 5 5 if not array_y_set_initial then (if 5 <= glob_max_terms 1, 6 1 then (temporary : array_tmp2 glob_h factorial_3(4, 5), 5 array_y : temporary, array_y_higher : temporary, 6 1, 6 temporary 2.0 temporary : -------------, array_y_higher : temporary)), kkk : 6, glob_h 2, 5 while kkk <= glob_max_terms do (array_tmp1 : kkk ats(kkk, array_y, array_y, 1), array_tmp2 : kkk array_tmp1 + array_const_0D0 , order_d : 1, kkk kkk if 1 + order_d + kkk <= glob_max_terms then (if not array_y_set_initial 1, order_d + kkk order_d array_tmp2 glob_h kkk then (temporary : -----------------------------------------, factorial_3(kkk - 1, - 1 + order_d + kkk) array_y : temporary, array_y_higher : temporary, order_d + kkk 1, order_d + kkk term : - 1 + order_d + kkk, adj2 : 2, while (adj2 <= 1 + order_d) temporary convfp(adj2) and (term >= 1) do (temporary : ----------------------, glob_h array_y_higher : temporary, adj2 : 1 + adj2, term : term - 1))), adj2, term kkk : 1 + kkk)) log(x) (%i9) log10(x) := --------- log(10.0) log(x) (%o9) log10(x) := --------- log(10.0) (%i10) omniout_str(iolevel, str) := if glob_iolevel >= iolevel then printf(true, "~a~%", string(str)) (%o10) omniout_str(iolevel, str) := if glob_iolevel >= iolevel then printf(true, "~a~%", string(str)) (%i11) omniout_str_noeol(iolevel, str) := if glob_iolevel >= iolevel then printf(true, "~a", string(str)) (%o11) omniout_str_noeol(iolevel, str) := if glob_iolevel >= iolevel then printf(true, "~a", string(str)) (%i12) omniout_labstr(iolevel, label, str) := if glob_iolevel >= iolevel then printf(true, "~a = ~a~%", string(label), string(str)) (%o12) omniout_labstr(iolevel, label, str) := if glob_iolevel >= iolevel then printf(true, "~a = ~a~%", string(label), string(str)) (%i13) 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)) (%o13) 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)) (%i14) omniout_int(iolevel, prelabel, prelen, value, vallen, postlabel) := if glob_iolevel >= iolevel then (printf(true, "~a = ~d ~a~%", prelabel, value, postlabel), newline()) (%o14) omniout_int(iolevel, prelabel, prelen, value, vallen, postlabel) := if glob_iolevel >= iolevel then (printf(true, "~a = ~d ~a~%", prelabel, value, postlabel), newline()) (%i15) omniout_float_arr(iolevel, prelabel, elemnt, prelen, value, vallen, postlabel) := if glob_iolevel >= iolevel then (sprint(prelabel, "[", elemnt, "]=", value, postlabel), newline()) (%o15) omniout_float_arr(iolevel, prelabel, elemnt, prelen, value, vallen, postlabel) := if glob_iolevel >= iolevel then (sprint(prelabel, "[", elemnt, "]=", value, postlabel), newline()) (%i16) dump_series(iolevel, dump_label, series_name, array_series, numb) := 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 (%o16) dump_series(iolevel, dump_label, series_name, array_series, numb) := 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 (%i17) dump_series_2(iolevel, dump_label, series_name, array_series2, numb, subnum) := 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 (%o17) dump_series_2(iolevel, dump_label, series_name, array_series2, numb, subnum) := 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 (%i18) 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)) (%o18) 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)) (%i19) logitem_time(fd, secs_in) := (secs : secs_in, printf(fd, ""), if secs >= 0.0 then (sec_in_millinium : sec_in_min min_in_hour hours_in_day days_in_year years_in_century secs centuries_in_millinium, milliniums : ----------------, sec_in_millinium millinium_int : floor(milliniums), centuries : (milliniums - millinium_int) centuries_in_millinium, cent_int : floor(centuries), years : (centuries - cent_int) years_in_century, years_int : floor(years), days : (years - years_int) days_in_year, days_int : floor(days), hours : (days - days_int) hours_in_day, hours_int : floor(hours), minutes : (hours - hours_int) min_in_hour, minutes_int : floor(minutes), seconds : (minutes - minutes_int) sec_in_min, sec_int : floor(seconds), if millinium_int > 0 then printf(fd, "~d Millinia ~d\ Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds", millinium_int, cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif cent_int > 0 then printf(fd, "~d Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds", cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif 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, "")) (%o19) logitem_time(fd, secs_in) := (secs : secs_in, printf(fd, ""), if secs >= 0.0 then (sec_in_millinium : sec_in_min min_in_hour hours_in_day days_in_year years_in_century secs centuries_in_millinium, milliniums : ----------------, sec_in_millinium millinium_int : floor(milliniums), centuries : (milliniums - millinium_int) centuries_in_millinium, cent_int : floor(centuries), years : (centuries - cent_int) years_in_century, years_int : floor(years), days : (years - years_int) days_in_year, days_int : floor(days), hours : (days - days_int) hours_in_day, hours_int : floor(hours), minutes : (hours - hours_int) min_in_hour, minutes_int : floor(minutes), seconds : (minutes - minutes_int) sec_in_min, sec_int : floor(seconds), if millinium_int > 0 then printf(fd, "~d Millinia ~d\ Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds", millinium_int, cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif cent_int > 0 then printf(fd, "~d Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds", cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif 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, "")) (%i20) omniout_timestr(secs_in) := (secs : convfloat(secs_in), if secs >= convfloat(0.0) then (sec_in_millinium : convfloat(sec_in_min) convfloat(min_in_hour) convfloat(hours_in_day) convfloat(days_in_year) convfloat(years_in_century) secs convfloat(centuries_in_millinium), milliniums : ---------------------------, convfloat(sec_in_millinium) millinium_int : floor(milliniums), centuries : (milliniums - millinium_int) convfloat(centuries_in_millinium), cent_int : floor(centuries), years : (centuries - cent_int) convfloat(years_in_century), years_int : floor(years), days : (years - years_int) convfloat(days_in_year), days_int : floor(days), hours : (days - days_int) convfloat(hours_in_day), hours_int : floor(hours), minutes : (hours - hours_int) convfloat(min_in_hour), minutes_int : floor(minutes), seconds : (minutes - minutes_int) convfloat(sec_in_min), sec_int : floor(seconds), if millinium_int > 0 then printf(true, "= ~d Millinia ~d Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", millinium_int, cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif cent_int > 0 then printf(true, "= ~d Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif 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~%")) (%o20) omniout_timestr(secs_in) := (secs : convfloat(secs_in), if secs >= convfloat(0.0) then (sec_in_millinium : convfloat(sec_in_min) convfloat(min_in_hour) convfloat(hours_in_day) convfloat(days_in_year) convfloat(years_in_century) secs convfloat(centuries_in_millinium), milliniums : ---------------------------, convfloat(sec_in_millinium) millinium_int : floor(milliniums), centuries : (milliniums - millinium_int) convfloat(centuries_in_millinium), cent_int : floor(centuries), years : (centuries - cent_int) convfloat(years_in_century), years_int : floor(years), days : (years - years_int) convfloat(days_in_year), days_int : floor(days), hours : (days - days_int) convfloat(hours_in_day), hours_int : floor(hours), minutes : (hours - hours_int) convfloat(min_in_hour), minutes_int : floor(minutes), seconds : (minutes - minutes_int) convfloat(sec_in_min), sec_int : floor(seconds), if millinium_int > 0 then printf(true, "= ~d Millinia ~d Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", millinium_int, cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif cent_int > 0 then printf(true, "= ~d Centuries ~d Years ~d Days ~d Hours ~d Minutes ~d Seconds~%", cent_int, years_int, days_int, hours_int, minutes_int, sec_int) elseif 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~%")) (%i21) mode_declare(ats, bfloat) modedeclare: bfloat is not a built-in type; assuming it is a Maxima extension type. (%o21) [ats] (%i22) ats(mmm_ats, array_a, array_b, jjj_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 : array_a array_b + ret_ats, iii_ats : 1 + iii_ats)), iii_ats lll_ats ret_ats) (%o22) ats(mmm_ats, array_a, array_b, jjj_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 : array_a array_b + ret_ats, iii_ats : 1 + iii_ats)), iii_ats lll_ats ret_ats) (%i23) mode_declare(att, bfloat) modedeclare: bfloat is not a built-in type; assuming it is a Maxima extension type. (%o23) [att] (%i24) att(mmm_att, array_aa, array_bb, jjj_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 : array_aa array_bb convfp(al_att) + ret_att, iii_att lll_att ret_att iii_att : 1 + iii_att), ret_att : ---------------), ret_att) convfp(mmm_att) (%o24) att(mmm_att, array_aa, array_bb, jjj_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 : array_aa array_bb convfp(al_att) + ret_att, iii_att lll_att ret_att iii_att : 1 + iii_att), ret_att : ---------------), ret_att) convfp(mmm_att) (%i25) 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 (%o25) 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 (%i26) logditto(file) := (printf(file, ""), printf(file, "ditto"), printf(file, "")) (%o26) logditto(file) := (printf(file, ""), printf(file, "ditto"), printf(file, "")) (%i27) logitem_integer(file, n) := (printf(file, ""), printf(file, "~d", n), printf(file, "")) (%o27) logitem_integer(file, n) := (printf(file, ""), printf(file, "~d", n), printf(file, "")) (%i28) logitem_str(file, str) := (printf(file, ""), printf(file, str), printf(file, "")) (%o28) logitem_str(file, str) := (printf(file, ""), printf(file, str), printf(file, "")) (%i29) log_revs(file, revs) := printf(file, revs) (%o29) log_revs(file, revs) := printf(file, revs) (%i30) logitem_float(file, x) := (printf(file, ""), printf(file, "~g", x), printf(file, "")) (%o30) logitem_float(file, x) := (printf(file, ""), printf(file, "~g", x), printf(file, "")) (%i31) 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, "")) (%o31) 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, "")) (%i32) logstart(file) := printf(file, "") (%o32) logstart(file) := printf(file, "") (%i33) logend(file) := printf(file, "~%") (%o33) logend(file) := printf(file, "~%") (%i34) chk_data() := (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()) (%o34) chk_data() := (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()) (%i35) mode_declare(comp_expect_sec, bfloat) modedeclare: bfloat is not a built-in type; assuming it is a Maxima extension type. (%o35) [comp_expect_sec] (%i36) comp_expect_sec(t_end2, t_start2, t2, clock_sec) := (ms2 : clock_sec, sub1 : t_end2 - t_start2, sub2 : t2 - t_start2, if sub1 = 0.0 then sec_left : 0.0 else (if abs(sub2) > 0.0 sub1 then (rrr : ----, sec_left : rrr ms2 - ms2) else sec_left : 0.0), sec_left) sub2 (%o36) comp_expect_sec(t_end2, t_start2, t2, clock_sec) := (ms2 : clock_sec, sub1 : t_end2 - t_start2, sub2 : t2 - t_start2, if sub1 = 0.0 then sec_left : 0.0 else (if abs(sub2) > 0.0 sub1 then (rrr : ----, sec_left : rrr ms2 - ms2) else sec_left : 0.0), sec_left) sub2 (%i37) mode_declare(comp_percent, bfloat) modedeclare: bfloat is not a built-in type; assuming it is a Maxima extension type. (%o37) [comp_percent] (%i38) comp_percent(t_end2, t_start2, t2) := (sub1 : t_end2 - t_start2, sub2 : t2 - t_start2, 100.0 sub2 if abs(sub2) > glob_small_float then rrr : ---------- else rrr : 0.0, rrr) sub1 (%o38) comp_percent(t_end2, t_start2, t2) := (sub1 : t_end2 - t_start2, sub2 : t2 - t_start2, 100.0 sub2 if abs(sub2) > glob_small_float then rrr : ---------- else rrr : 0.0, rrr) sub1 (%i39) mode_declare(factorial_1, bfloat) modedeclare: bfloat is not a built-in type; assuming it is a Maxima extension type. (%o39) [factorial_1] (%i40) factorial_1(nnn) := nnn! (%o40) factorial_1(nnn) := nnn! (%i41) mode_declare(factorial_3, bfloat) modedeclare: bfloat is not a built-in type; assuming it is a Maxima extension type. (%o41) [factorial_3] mmm2! (%i42) factorial_3(mmm2, nnn2) := ----- nnn2! mmm2! (%o42) factorial_3(mmm2, nnn2) := ----- nnn2! (%i43) convfp(mmm) := mmm (%o43) convfp(mmm) := mmm (%i44) convfloat(mmm) := mmm (%o44) convfloat(mmm) := mmm (%i45) elapsed_time_seconds() := (t : elapsed_real_time(), t) (%o45) elapsed_time_seconds() := (t : elapsed_real_time(), t) (%i46) arcsin(x) := asin(x) (%o46) arcsin(x) := asin(x) (%i47) arccos(x) := acos(x) (%o47) arccos(x) := acos(x) (%i48) arctan(x) := atan(x) (%o48) arctan(x) := atan(x) 1.0 (%i49) exact_soln_y(x) := ------- 1.0 - x 1.0 (%o49) exact_soln_y(x) := ------- 1.0 - x (%i50) mainprog() := (define_variable(DEBUGMASSIVE, 4, fixnum), define_variable(DEBUGL, 3, fixnum), define_variable(glob_max_terms, 30, fixnum), define_variable(ALWAYS, 1, fixnum), define_variable(INFO, 2, fixnum), define_variable(glob_iolevel, 5, fixnum), define_variable(glob_max_minutes, 0.0, float), define_variable(glob_log10abserr, 0.0, float), define_variable(glob_last_good_h, 0.1, float), define_variable(glob_html_log, true, boolean), define_variable(glob_percent_done, 0.0, float), define_variable(glob_orig_start_sec, 0.0, float), define_variable(glob_max_hours, 0.0, float), define_variable(glob_relerr, 1.0E-11, float), define_variable(glob_subiter_method, 3, fixnum), define_variable(MAX_UNCHANGED, 10, fixnum), define_variable(glob_small_float, 1.0E-51, float), define_variable(glob_max_trunc_err, 1.0E-11, float), define_variable(glob_max_rel_trunc_err, 1.0E-11, float), define_variable(glob_log10_relerr, 1.0E-11, float), define_variable(glob_h, 0.1, float), define_variable(glob_clock_sec, 0.0, float), define_variable(centuries_in_millinium, 10.0, float), define_variable(sec_in_min, 60.0, float), define_variable(glob_log10normmin, 0.1, float), define_variable(glob_log10_abserr, 1.0E-11, float), define_variable(glob_dump_analytic, false, boolean), define_variable(glob_disp_incr, 0.1, float), define_variable(glob_almost_1, 0.999, float), define_variable(min_in_hour, 60.0, float), define_variable(djd_debug2, true, boolean), define_variable(glob_display_flag, true, boolean), define_variable(glob_dump, false, boolean), define_variable(glob_start, 0, fixnum), define_variable(glob_smallish_float, 1.0E-101, float), define_variable(glob_unchanged_h_cnt, 0, fixnum), define_variable(glob_no_eqs, 0, fixnum), define_variable(glob_look_poles, false, boolean), define_variable(glob_max_opt_iter, 10, fixnum), define_variable(glob_optimal_clock_start_sec, 0.0, float), define_variable(glob_optimal_done, false, boolean), define_variable(glob_not_yet_finished, true, boolean), define_variable(years_in_century, 100.0, float), define_variable(glob_warned, false, boolean), define_variable(glob_large_float, 9.0E+100, float), define_variable(glob_hmin, 1.0E-11, float), define_variable(glob_not_yet_start_msg, true, boolean), define_variable(glob_initial_pass, true, boolean), define_variable(glob_clock_start_sec, 0.0, float), define_variable(djd_debug, true, boolean), define_variable(glob_current_iter, 0, fixnum), define_variable(glob_curr_iter_when_opt, 0, fixnum), define_variable(glob_hmin_init, 0.001, float), define_variable(glob_reached_optimal_h, false, boolean), define_variable(days_in_year, 365.0, float), define_variable(glob_iter, 0, fixnum), define_variable(glob_max_sec, 10000.0, float), define_variable(glob_optimal_start, 0.0, float), define_variable(glob_optimal_expect_sec, 0.1, float), define_variable(glob_log10relerr, 0.0, float), define_variable(glob_normmax, 0.0, float), define_variable(glob_warned2, false, boolean), define_variable(glob_max_iter, 1000, fixnum), define_variable(glob_abserr, 1.0E-11, float), define_variable(glob_hmax, 1.0, float), define_variable(hours_in_day, 24.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/nonlinear1postode.ode#################"), omniout_str(ALWAYS, "diff ( y , x , 1 ) = y * y;"), 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 : 0.0,"), omniout_str(ALWAYS, "x_end : 0.5 ,"), omniout_str(ALWAYS, "array_y_init[0 + 1] : exact_soln_y(x_start),"), omniout_str(ALWAYS, "glob_h : 0.01,"), omniout_str(ALWAYS, "glob_look_poles : true,"), omniout_str(ALWAYS, "glob_max_iter : 1000000,"), omniout_str(ALWAYS, "/* END SECOND INPUT BLOCK */"), omniout_str(ALWAYS, "/* BEGIN OVERRIDE BLOCK */"), omniout_str(ALWAYS, "glob_h : 0.001 ,"), omniout_str(ALWAYS, "glob_look_poles : true,"), omniout_str(ALWAYS, "glob_max_iter : 1000,"), omniout_str(ALWAYS, "glob_max_minutes : 15,"), omniout_str(ALWAYS, "/* END OVERRIDE BLOCK */"), omniout_str(ALWAYS, "!"), omniout_str(ALWAYS, "/* BEGIN USER DEF BLOCK */"), omniout_str(ALWAYS, "exact_soln_y (x) := ("), omniout_str(ALWAYS, "1.0/(1.0 - x) "), 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_m1, 1 + max_terms), array(array_last_rel_error, 1 + max_terms), array(array_1st_rel_error, 1 + max_terms), array(array_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_y_init, 1 + max_terms), array(array_type_pole, 1 + max_terms), array(array_norms, 1 + max_terms), array(array_poles, 1 + 1, 1 + 3), array(array_y_higher_work2, 1 + 2, 1 + max_terms), array(array_y_higher, 1 + 2, 1 + max_terms), array(array_y_set_initial, 1 + 2, 1 + max_terms), array(array_y_higher_work, 1 + 2, 1 + max_terms), array(array_complex_pole, 1 + 1, 1 + 3), array(array_real_pole, 1 + 1, 1 + 3), term : 1, while term <= max_terms do (array_m1 : 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_1st_rel_error : 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_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_y_init : 0.0, term : 1 + term), term term : 1, while term <= max_terms do (array_type_pole : 0.0, term term : 1 + term), term : 1, while term <= max_terms do (array_norms : 0.0, term : 1 + term), ord : 1, term 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 <= 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_higher : 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 <= 2 do (term : 1, while term <= max_terms do (array_y_higher_work : 0.0, term : 1 + term), ord, term ord : 1 + ord), 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 <= 1 do (term : 1, while term <= 3 do (array_real_pole : 0.0, term : 1 + term), ord, term ord : 1 + ord), 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_y, 1 + 1 + max_terms), term : 1, while term <= 1 + max_terms do (array_y : 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_tmp1, 1 + 1 + max_terms), term : 1, while term <= 1 + max_terms do (array_tmp1 : 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_const_0D0, 1 + 1 + max_terms), term : 1, while term <= 1 + max_terms do (array_const_0D0 : 0.0, term : 1 + term), term array_const_0D0 : 0.0, array(array_const_1, 1 + 1 + max_terms), term : 1, 1 while term <= 1 + max_terms do (array_const_1 : 0.0, term : 1 + term), term array_const_1 : 1, 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, x_start : 0.0, x_end : 0.5, 1 array_y_init : exact_soln_y(x_start), glob_h : 0.01, 1 + 0 glob_look_poles : true, glob_max_iter : 1000000, glob_h : 0.001, glob_look_poles : true, glob_max_iter : 1000, glob_max_minutes : 15, 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_log10_abserr glob_log10_relerr glob_abserr : 10.0 , glob_relerr : 10.0 , chk_data(), array_y_set_initial : true, array_y_set_initial : false, 1, 1 1, 2 array_y_set_initial : false, array_y_set_initial : false, 1, 3 1, 4 array_y_set_initial : false, array_y_set_initial : false, 1, 5 1, 6 array_y_set_initial : false, array_y_set_initial : false, 1, 7 1, 8 array_y_set_initial : false, array_y_set_initial : false, 1, 9 1, 10 array_y_set_initial : false, array_y_set_initial : false, 1, 11 1, 12 array_y_set_initial : false, array_y_set_initial : false, 1, 13 1, 14 array_y_set_initial : false, array_y_set_initial : false, 1, 15 1, 16 array_y_set_initial : false, array_y_set_initial : false, 1, 17 1, 18 array_y_set_initial : false, array_y_set_initial : false, 1, 19 1, 20 array_y_set_initial : false, array_y_set_initial : false, 1, 21 1, 22 array_y_set_initial : false, array_y_set_initial : false, 1, 23 1, 24 array_y_set_initial : false, array_y_set_initial : false, 1, 25 1, 26 array_y_set_initial : false, array_y_set_initial : false, 1, 27 1, 28 array_y_set_initial : false, array_y_set_initial : false, 1, 29 1, 30 if glob_html_log then html_log_file : openw("html/entry.html"), omniout_str(ALWAYS, "START of Soultion"), array_x : x_start, 1 array_x : glob_h, order_diff : 1, term_no : 1, 2 while term_no <= order_diff do (array_y : term_no term_no - 1 array_y_init glob_h 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, term_no - 1 array_y_init glob_h 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(), start_array_y(), if !array_y_higher ! > glob_small_float ! 1, 1! then (tmp : !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 (array_x <= x_end) and (convfloat(glob_clock_sec) - convfloat(glob_orig_start_sec) < 1 convfloat(glob_max_sec)) do (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(), if glob_look_poles then check_for_pole(), array_x : glob_h + array_x , 1 1 array_x : glob_h, order_diff : 1, 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 -------------------- calc_term - 1 glob_h -------------------------------------, 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 calc_term - 1 temp_sum glob_h ----------------------------, ord : 1, calc_term : 2, iii : glob_max_terms, convfp(calc_term - 1)! while iii >= calc_term do (array_y_higher_work : 1, iii array_y_higher 1, iii -------------------- calc_term - 1 glob_h -------------------------------------, 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 calc_term - 1 temp_sum glob_h ----------------------------, ord : 1, calc_term : 1, iii : glob_max_terms, convfp(calc_term - 1)! while iii >= calc_term do (array_y_higher_work : 1, iii array_y_higher 1, iii -------------------- calc_term - 1 glob_h -------------------------------------, 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 calc_term - 1 temp_sum glob_h ----------------------------, term_no : glob_max_terms, convfp(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 display_alot(current_iter)), 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 ) = y * y;"), 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, "2012-06-16T00:20:48-05:00"), logitem_str(html_log_file, "Maxima"), logitem_str(html_log_file, "nonlinear1"), logitem_str(html_log_file, "diff ( y , x , 1 ) = y * y;"), 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_integer(html_log_file, glob_max_terms), 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_optimal_expect_sec)), 0) else (logitem_str(html_log_file, "Done"), 0), log_revs(html_log_file, " 090 "), logitem_str(html_log_file, "nonlinear1 diffeq.max"), logitem_str(html_log_file, "nonlinear1 maxima results"), logitem_str(html_log_file, "Test of revised logic - mostly affecting systems of eqs"), logend(html_log_file)), if glob_html_log then close(html_log_file)) (%o50) mainprog() := (define_variable(DEBUGMASSIVE, 4, fixnum), define_variable(DEBUGL, 3, fixnum), define_variable(glob_max_terms, 30, fixnum), define_variable(ALWAYS, 1, fixnum), define_variable(INFO, 2, fixnum), define_variable(glob_iolevel, 5, fixnum), define_variable(glob_max_minutes, 0.0, float), define_variable(glob_log10abserr, 0.0, float), define_variable(glob_last_good_h, 0.1, float), define_variable(glob_html_log, true, boolean), define_variable(glob_percent_done, 0.0, float), define_variable(glob_orig_start_sec, 0.0, float), define_variable(glob_max_hours, 0.0, float), define_variable(glob_relerr, 1.0E-11, float), define_variable(glob_subiter_method, 3, fixnum), define_variable(MAX_UNCHANGED, 10, fixnum), define_variable(glob_small_float, 1.0E-51, float), define_variable(glob_max_trunc_err, 1.0E-11, float), define_variable(glob_max_rel_trunc_err, 1.0E-11, float), define_variable(glob_log10_relerr, 1.0E-11, float), define_variable(glob_h, 0.1, float), define_variable(glob_clock_sec, 0.0, float), define_variable(centuries_in_millinium, 10.0, float), define_variable(sec_in_min, 60.0, float), define_variable(glob_log10normmin, 0.1, float), define_variable(glob_log10_abserr, 1.0E-11, float), define_variable(glob_dump_analytic, false, boolean), define_variable(glob_disp_incr, 0.1, float), define_variable(glob_almost_1, 0.999, float), define_variable(min_in_hour, 60.0, float), define_variable(djd_debug2, true, boolean), define_variable(glob_display_flag, true, boolean), define_variable(glob_dump, false, boolean), define_variable(glob_start, 0, fixnum), define_variable(glob_smallish_float, 1.0E-101, float), define_variable(glob_unchanged_h_cnt, 0, fixnum), define_variable(glob_no_eqs, 0, fixnum), define_variable(glob_look_poles, false, boolean), define_variable(glob_max_opt_iter, 10, fixnum), define_variable(glob_optimal_clock_start_sec, 0.0, float), define_variable(glob_optimal_done, false, boolean), define_variable(glob_not_yet_finished, true, boolean), define_variable(years_in_century, 100.0, float), define_variable(glob_warned, false, boolean), define_variable(glob_large_float, 9.0E+100, float), define_variable(glob_hmin, 1.0E-11, float), define_variable(glob_not_yet_start_msg, true, boolean), define_variable(glob_initial_pass, true, boolean), define_variable(glob_clock_start_sec, 0.0, float), define_variable(djd_debug, true, boolean), define_variable(glob_current_iter, 0, fixnum), define_variable(glob_curr_iter_when_opt, 0, fixnum), define_variable(glob_hmin_init, 0.001, float), define_variable(glob_reached_optimal_h, false, boolean), define_variable(days_in_year, 365.0, float), define_variable(glob_iter, 0, fixnum), define_variable(glob_max_sec, 10000.0, float), define_variable(glob_optimal_start, 0.0, float), define_variable(glob_optimal_expect_sec, 0.1, float), define_variable(glob_log10relerr, 0.0, float), define_variable(glob_normmax, 0.0, float), define_variable(glob_warned2, false, boolean), define_variable(glob_max_iter, 1000, fixnum), define_variable(glob_abserr, 1.0E-11, float), define_variable(glob_hmax, 1.0, float), define_variable(hours_in_day, 24.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/nonlinear1postode.ode#################"), omniout_str(ALWAYS, "diff ( y , x , 1 ) = y * y;"), 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 : 0.0,"), omniout_str(ALWAYS, "x_end : 0.5 ,"), omniout_str(ALWAYS, "array_y_init[0 + 1] : exact_soln_y(x_start),"), omniout_str(ALWAYS, "glob_h : 0.01,"), omniout_str(ALWAYS, "glob_look_poles : true,"), omniout_str(ALWAYS, "glob_max_iter : 1000000,"), omniout_str(ALWAYS, "/* END SECOND INPUT BLOCK */"), omniout_str(ALWAYS, "/* BEGIN OVERRIDE BLOCK */"), omniout_str(ALWAYS, "glob_h : 0.001 ,"), omniout_str(ALWAYS, "glob_look_poles : true,"), omniout_str(ALWAYS, "glob_max_iter : 1000,"), omniout_str(ALWAYS, "glob_max_minutes : 15,"), omniout_str(ALWAYS, "/* END OVERRIDE BLOCK */"), omniout_str(ALWAYS, "!"), omniout_str(ALWAYS, "/* BEGIN USER DEF BLOCK */"), omniout_str(ALWAYS, "exact_soln_y (x) := ("), omniout_str(ALWAYS, "1.0/(1.0 - x) "), 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_m1, 1 + max_terms), array(array_last_rel_error, 1 + max_terms), array(array_1st_rel_error, 1 + max_terms), array(array_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_y_init, 1 + max_terms), array(array_type_pole, 1 + max_terms), array(array_norms, 1 + max_terms), array(array_poles, 1 + 1, 1 + 3), array(array_y_higher_work2, 1 + 2, 1 + max_terms), array(array_y_higher, 1 + 2, 1 + max_terms), array(array_y_set_initial, 1 + 2, 1 + max_terms), array(array_y_higher_work, 1 + 2, 1 + max_terms), array(array_complex_pole, 1 + 1, 1 + 3), array(array_real_pole, 1 + 1, 1 + 3), term : 1, while term <= max_terms do (array_m1 : 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_1st_rel_error : 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_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_y_init : 0.0, term : 1 + term), term term : 1, while term <= max_terms do (array_type_pole : 0.0, term term : 1 + term), term : 1, while term <= max_terms do (array_norms : 0.0, term : 1 + term), ord : 1, term 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 <= 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_higher : 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 <= 2 do (term : 1, while term <= max_terms do (array_y_higher_work : 0.0, term : 1 + term), ord, term ord : 1 + ord), 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 <= 1 do (term : 1, while term <= 3 do (array_real_pole : 0.0, term : 1 + term), ord, term ord : 1 + ord), 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_y, 1 + 1 + max_terms), term : 1, while term <= 1 + max_terms do (array_y : 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_tmp1, 1 + 1 + max_terms), term : 1, while term <= 1 + max_terms do (array_tmp1 : 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_const_0D0, 1 + 1 + max_terms), term : 1, while term <= 1 + max_terms do (array_const_0D0 : 0.0, term : 1 + term), term array_const_0D0 : 0.0, array(array_const_1, 1 + 1 + max_terms), term : 1, 1 while term <= 1 + max_terms do (array_const_1 : 0.0, term : 1 + term), term array_const_1 : 1, 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, x_start : 0.0, x_end : 0.5, 1 array_y_init : exact_soln_y(x_start), glob_h : 0.01, 1 + 0 glob_look_poles : true, glob_max_iter : 1000000, glob_h : 0.001, glob_look_poles : true, glob_max_iter : 1000, glob_max_minutes : 15, 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_log10_abserr glob_log10_relerr glob_abserr : 10.0 , glob_relerr : 10.0 , chk_data(), array_y_set_initial : true, array_y_set_initial : false, 1, 1 1, 2 array_y_set_initial : false, array_y_set_initial : false, 1, 3 1, 4 array_y_set_initial : false, array_y_set_initial : false, 1, 5 1, 6 array_y_set_initial : false, array_y_set_initial : false, 1, 7 1, 8 array_y_set_initial : false, array_y_set_initial : false, 1, 9 1, 10 array_y_set_initial : false, array_y_set_initial : false, 1, 11 1, 12 array_y_set_initial : false, array_y_set_initial : false, 1, 13 1, 14 array_y_set_initial : false, array_y_set_initial : false, 1, 15 1, 16 array_y_set_initial : false, array_y_set_initial : false, 1, 17 1, 18 array_y_set_initial : false, array_y_set_initial : false, 1, 19 1, 20 array_y_set_initial : false, array_y_set_initial : false, 1, 21 1, 22 array_y_set_initial : false, array_y_set_initial : false, 1, 23 1, 24 array_y_set_initial : false, array_y_set_initial : false, 1, 25 1, 26 array_y_set_initial : false, array_y_set_initial : false, 1, 27 1, 28 array_y_set_initial : false, array_y_set_initial : false, 1, 29 1, 30 if glob_html_log then html_log_file : openw("html/entry.html"), omniout_str(ALWAYS, "START of Soultion"), array_x : x_start, 1 array_x : glob_h, order_diff : 1, term_no : 1, 2 while term_no <= order_diff do (array_y : term_no term_no - 1 array_y_init glob_h 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, term_no - 1 array_y_init glob_h 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(), start_array_y(), if !array_y_higher ! > glob_small_float ! 1, 1! then (tmp : !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 (array_x <= x_end) and (convfloat(glob_clock_sec) - convfloat(glob_orig_start_sec) < 1 convfloat(glob_max_sec)) do (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(), if glob_look_poles then check_for_pole(), array_x : glob_h + array_x , 1 1 array_x : glob_h, order_diff : 1, 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 -------------------- calc_term - 1 glob_h -------------------------------------, 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 calc_term - 1 temp_sum glob_h ----------------------------, ord : 1, calc_term : 2, iii : glob_max_terms, convfp(calc_term - 1)! while iii >= calc_term do (array_y_higher_work : 1, iii array_y_higher 1, iii -------------------- calc_term - 1 glob_h -------------------------------------, 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 calc_term - 1 temp_sum glob_h ----------------------------, ord : 1, calc_term : 1, iii : glob_max_terms, convfp(calc_term - 1)! while iii >= calc_term do (array_y_higher_work : 1, iii array_y_higher 1, iii -------------------- calc_term - 1 glob_h -------------------------------------, 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 calc_term - 1 temp_sum glob_h ----------------------------, term_no : glob_max_terms, convfp(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 display_alot(current_iter)), 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 ) = y * y;"), 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, "2012-06-16T00:20:48-05:00"), logitem_str(html_log_file, "Maxima"), logitem_str(html_log_file, "nonlinear1"), logitem_str(html_log_file, "diff ( y , x , 1 ) = y * y;"), 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_integer(html_log_file, glob_max_terms), 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_optimal_expect_sec)), 0) else (logitem_str(html_log_file, "Done"), 0), log_revs(html_log_file, " 090 "), logitem_str(html_log_file, "nonlinear1 diffeq.max"), logitem_str(html_log_file, "nonlinear1 maxima results"), logitem_str(html_log_file, "Test of revised logic - mostly affecting systems of eqs"), logend(html_log_file)), if glob_html_log then close(html_log_file)) (%i51) mainprog() "##############ECHO OF PROBLEM#################" "##############temp/nonlinear1postode.ode#################" "diff ( y , x , 1 ) = y * y;" "!" "/* BEGIN FIRST INPUT BLOCK */" "Digits : 32," "max_terms : 30," "!" "/* END FIRST INPUT BLOCK */" "/* BEGIN SECOND INPUT BLOCK */" "x_start : 0.0," "x_end : 0.5 ," "array_y_init[0 + 1] : exact_soln_y(x_start)," "glob_h : 0.01," "glob_look_poles : true," "glob_max_iter : 1000000," "/* END SECOND INPUT BLOCK */" "/* BEGIN OVERRIDE BLOCK */" "glob_h : 0.001 ," "glob_look_poles : true," "glob_max_iter : 1000," "glob_max_minutes : 15," "/* END OVERRIDE BLOCK */" "!" "/* BEGIN USER DEF BLOCK */" "exact_soln_y (x) := (" "1.0/(1.0 - x) " ");" "/* END USER DEF BLOCK */" "#######END OF ECHO OF PROBLEM#################" "START of Soultion" x[1] = 0.0 " " y[1] (analytic) = 1. " " y[1] (numeric) = 1. " " absolute error = 0.0 " " relative error = 0.0 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.000E-3 " " y[1] (analytic) = 1.001001001001001 " " y[1] (numeric) = 1.001001001001001 " " absolute error = 0.0 " " relative error = 0.0 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.000E-3 " " y[1] (analytic) = 1.002004008016032 " " y[1] (numeric) = 1.0020040080160322 " " absolute error = 2.2204460492503130000000000000000E-16 " " relative error = 2.216005157151812700000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.000E-3 " " y[1] (analytic) = 1.0030090270812437 " " y[1] (numeric) = 1.003009027081244 " " absolute error = 2.2204460492503130000000000000000E-16 " " relative error = 2.213784711102562100000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.000E-3 " " y[1] (analytic) = 1.0040160642570282 " " y[1] (numeric) = 1.0040160642570284 " " absolute error = 2.2204460492503130000000000000000E-16 " " relative error = 2.211564265053311800000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.000E-3 " " y[1] (analytic) = 1.0050251256281406 " " y[1] (numeric) = 1.005025125628141 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.418687638008123600000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.000E-3 " " y[1] (analytic) = 1.0060362173038229 " " y[1] (numeric) = 1.0060362173038233 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.414246745909622400000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.000E-3 " " y[1] (analytic) = 1.0070493454179255 " " y[1] (numeric) = 1.007049345417926 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.409805853811122000000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.000E-3 " " y[1] (analytic) = 1.0080645161290323 " " y[1] (numeric) = 1.0080645161290327 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.40536496171262100000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.000000000000001000E-3 " " y[1] (analytic) = 1.0090817356205852 " " y[1] (numeric) = 1.0090817356205857 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.40092406961412100000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.000000000000000200E-2 " " y[1] (analytic) = 1.0101010101010102 " " y[1] (numeric) = 1.0101010101010106 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.3964831775156200000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.100000000000000300E-2 " " y[1] (analytic) = 1.0111223458038423 " " y[1] (numeric) = 1.0111223458038427 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.39204228541711870000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.200000000000000400E-2 " " y[1] (analytic) = 1.0121457489878543 " " y[1] (numeric) = 1.0121457489878547 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.387601393318618600000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.300000000000000600E-2 " " y[1] (analytic) = 1.0131712259371835 " " y[1] (numeric) = 1.013171225937184 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.383160501220117500000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.400000000000000700E-2 " " y[1] (analytic) = 1.0141987829614605 " " y[1] (numeric) = 1.014198782961461 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.37871960912161740000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.500000000000000800E-2 " " y[1] (analytic) = 1.015228426395939 " " y[1] (numeric) = 1.0152284263959397 " " absolute error = 6.6613381477509390000000000000000E-16 " " relative error = 6.56141807553467500000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.600000000000001000E-2 " " y[1] (analytic) = 1.016260162601626 " " y[1] (numeric) = 1.0162601626016265 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.36983782492461560000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.700000000000001000E-2 " " y[1] (analytic) = 1.017293997965412 " " y[1] (numeric) = 1.0172939979654125 " " absolute error = 4.4408920985006260000000000000000E-16 " " relative error = 4.365396932826115500000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.80000000000000100E-2 " " y[1] (analytic) = 1.0183299389002036 " " y[1] (numeric) = 1.0183299389002043 " " absolute error = 6.6613381477509390000000000000000E-16 " " relative error = 6.54143406109142200000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 1.90000000000000100E-2 " " y[1] (analytic) = 1.019367991845056 " " y[1] (numeric) = 1.0193679918450567 " " absolute error = 6.6613381477509390000000000000000E-16 " " relative error = 6.53477272294367100000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.00000000000000120E-2 " " y[1] (analytic) = 1.0204081632653061 " " y[1] (numeric) = 1.0204081632653068 " " absolute error = 6.6613381477509390000000000000000E-16 " " relative error = 6.5281113847959200000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.10000000000000130E-2 " " y[1] (analytic) = 1.0214504596527068 " " y[1] (numeric) = 1.0214504596527076 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.69526672886422600000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.200000000000001400E-2 " " y[1] (analytic) = 1.0224948875255624 " " y[1] (numeric) = 1.0224948875255633 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.68638494466722500000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.300000000000001500E-2 " " y[1] (analytic) = 1.0235414534288638 " " y[1] (numeric) = 1.0235414534288647 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.67750316047022400000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.400000000000001600E-2 " " y[1] (analytic) = 1.0245901639344261 " " y[1] (numeric) = 1.024590163934427 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.66862137627322300000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.500000000000001700E-2 " " y[1] (analytic) = 1.0256410256410258 " " y[1] (numeric) = 1.0256410256410264 " " absolute error = 6.6613381477509390000000000000000E-16 " " relative error = 6.49480469405716500000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.600000000000002000E-2 " " y[1] (analytic) = 1.026694045174538 " " y[1] (numeric) = 1.0266940451745388 " " absolute error = 6.6613381477509390000000000000000E-16 " " relative error = 6.48814335590941500000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.700000000000002000E-2 " " y[1] (analytic) = 1.027749229188078 " " y[1] (numeric) = 1.027749229188079 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.64197602368221900000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.800000000000002000E-2 " " y[1] (analytic) = 1.02880658436214 " " y[1] (numeric) = 1.0288065843621408 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.63309423948521700000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 2.90000000000000200E-2 " " y[1] (analytic) = 1.0298661174047374 " " y[1] (numeric) = 1.0298661174047383 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.62421245528821600000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.00000000000000200E-2 " " y[1] (analytic) = 1.0309278350515465 " " y[1] (numeric) = 1.0309278350515474 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.61533067109121400000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.10000000000000200E-2 " " y[1] (analytic) = 1.0319917440660475 " " y[1] (numeric) = 1.0319917440660484 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.60644888689421200000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.20000000000000230E-2 " " y[1] (analytic) = 1.0330578512396695 " " y[1] (numeric) = 1.0330578512396704 " " absolute error = 8.8817841970012520000000000000000E-16 " " relative error = 8.59756710269721100000000000000E-14 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.30000000000000240E-2 " " y[1] (analytic) = 1.0341261633919339 " " y[1] (numeric) = 1.034126163391935 " " absolute error = 1.1102230246251565000000000000000E-15 " " relative error = 1.07358566481252640000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.40000000000000250E-2 " " y[1] (analytic) = 1.0351966873706004 " " y[1] (numeric) = 1.0351966873706016 " " absolute error = 1.1102230246251565000000000000000E-15 " " relative error = 1.07247544178790120000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.500000000000002600E-2 " " y[1] (analytic) = 1.0362694300518136 " " y[1] (numeric) = 1.0362694300518147 " " absolute error = 1.1102230246251565000000000000000E-15 " " relative error = 1.07136521876327590000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.600000000000002600E-2 " " y[1] (analytic) = 1.037344398340249 " " y[1] (numeric) = 1.0373443983402502 " " absolute error = 1.1102230246251565000000000000000E-15 " " relative error = 1.07025499573865080000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.700000000000003000E-2 " " y[1] (analytic) = 1.0384215991692627 " " y[1] (numeric) = 1.038421599169264 " " absolute error = 1.3322676295501878000000000000000E-15 " " relative error = 1.2829737272568312000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.80000000000000300E-2 " " y[1] (analytic) = 1.0395010395010396 " " y[1] (numeric) = 1.039501039501041 " " absolute error = 1.3322676295501878000000000000000E-15 " " relative error = 1.28164145962728040000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 3.90000000000000300E-2 " " y[1] (analytic) = 1.040582726326743 " " y[1] (numeric) = 1.0405827263267444 " " absolute error = 1.3322676295501878000000000000000E-15 " " relative error = 1.28030919199773020000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.00000000000000300E-2 " " y[1] (analytic) = 1.0416666666666667 " " y[1] (numeric) = 1.041666666666668 " " absolute error = 1.3322676295501878000000000000000E-15 " " relative error = 1.27897692436818030000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.10000000000000300E-2 " " y[1] (analytic) = 1.0427528675703859 " " y[1] (numeric) = 1.0427528675703872 " " absolute error = 1.3322676295501878000000000000000E-15 " " relative error = 1.277644656738630100000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.20000000000000300E-2 " " y[1] (analytic) = 1.0438413361169103 " " y[1] (numeric) = 1.0438413361169117 " " absolute error = 1.3322676295501878000000000000000E-15 " " relative error = 1.27631238910908000000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.30000000000000300E-2 " " y[1] (analytic) = 1.044932079414838 " " y[1] (numeric) = 1.0449320794148396 " " absolute error = 1.5543122344752192000000000000000E-15 " " relative error = 1.48747680839278470000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.40000000000000340E-2 " " y[1] (analytic) = 1.0460251046025104 " " y[1] (numeric) = 1.046025104602512 " " absolute error = 1.5543122344752192000000000000000E-15 " " relative error = 1.48592249615830950000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.50000000000000340E-2 " " y[1] (analytic) = 1.0471204188481675 " " y[1] (numeric) = 1.047120418848169 " " absolute error = 1.5543122344752192000000000000000E-15 " " relative error = 1.48436818392383430000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.600000000000003500E-2 " " y[1] (analytic) = 1.0482180293501049 " " y[1] (numeric) = 1.0482180293501064 " " absolute error = 1.5543122344752192000000000000000E-15 " " relative error = 1.4828138716893590000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.700000000000003600E-2 " " y[1] (analytic) = 1.0493179433368311 " " y[1] (numeric) = 1.0493179433368327 " " absolute error = 1.5543122344752192000000000000000E-15 " " relative error = 1.48125955945488360000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.800000000000003700E-2 " " y[1] (analytic) = 1.050420168067227 " " y[1] (numeric) = 1.0504201680672285 " " absolute error = 1.5543122344752192000000000000000E-15 " " relative error = 1.47970524722040860000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 4.90000000000000400E-2 " " y[1] (analytic) = 1.0515247108307046 " " y[1] (numeric) = 1.0515247108307062 " " absolute error = 1.5543122344752192000000000000000E-15 " " relative error = 1.47815093498593310000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.00000000000000300E-2 " " y[1] (analytic) = 1.0526315789473684 " " y[1] (numeric) = 1.0526315789473701 " " absolute error = 1.7763568394002505000000000000000E-15 " " relative error = 1.6875389974302380000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.10000000000000300E-2 " " y[1] (analytic) = 1.053740779768177 " " y[1] (numeric) = 1.0537407797681788 " " absolute error = 1.7763568394002505000000000000000E-15 " " relative error = 1.68576264059083770000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.20000000000000400E-2 " " y[1] (analytic) = 1.0548523206751055 " " y[1] (numeric) = 1.0548523206751073 " " absolute error = 1.7763568394002505000000000000000E-15 " " relative error = 1.68398628375143740000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.30000000000000400E-2 " " y[1] (analytic) = 1.0559662090813096 " " y[1] (numeric) = 1.0559662090813113 " " absolute error = 1.7763568394002505000000000000000E-15 " " relative error = 1.6822099269120370000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.40000000000000400E-2 " " y[1] (analytic) = 1.0570824524312896 " " y[1] (numeric) = 1.0570824524312916 " " absolute error = 1.9984014443252818000000000000000E-15 " " relative error = 1.89048776633171660000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.50000000000000400E-2 " " y[1] (analytic) = 1.0582010582010584 " " y[1] (numeric) = 1.0582010582010601 " " absolute error = 1.7763568394002505000000000000000E-15 " " relative error = 1.67865721323323640000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.60000000000000400E-2 " " y[1] (analytic) = 1.0593220338983051 " " y[1] (numeric) = 1.0593220338983071 " " absolute error = 1.9984014443252818000000000000000E-15 " " relative error = 1.88649096344306570000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.700000000000004000E-2 " " y[1] (analytic) = 1.0604453870625663 " " y[1] (numeric) = 1.0604453870625683 " " absolute error = 1.9984014443252818000000000000000E-15 " " relative error = 1.88449256199874040000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.800000000000004000E-2 " " y[1] (analytic) = 1.0615711252653928 " " y[1] (numeric) = 1.0615711252653948 " " absolute error = 1.9984014443252818000000000000000E-15 " " relative error = 1.88249416055441540000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 5.900000000000004000E-2 " " y[1] (analytic) = 1.0626992561105209 " " y[1] (numeric) = 1.0626992561105228 " " absolute error = 1.9984014443252818000000000000000E-15 " " relative error = 1.88049575911009000000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.000000000000004000E-2 " " y[1] (analytic) = 1.0638297872340425 " " y[1] (numeric) = 1.0638297872340448 " " absolute error = 2.220446049250313000000000000000E-15 " " relative error = 2.08721928629529430000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.10000000000000400E-2 " " y[1] (analytic) = 1.0649627263045793 " " y[1] (numeric) = 1.0649627263045816 " " absolute error = 2.220446049250313000000000000000E-15 " " relative error = 2.0849988402460440000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.20000000000000400E-2 " " y[1] (analytic) = 1.0660980810234542 " " y[1] (numeric) = 1.0660980810234564 " " absolute error = 2.220446049250313000000000000000E-15 " " relative error = 2.08277839419679370000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.30000000000000400E-2 " " y[1] (analytic) = 1.0672358591248667 " " y[1] (numeric) = 1.067235859124869 " " absolute error = 2.220446049250313000000000000000E-15 " " relative error = 2.0805579481475430000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.40000000000000500E-2 " " y[1] (analytic) = 1.0683760683760684 " " y[1] (numeric) = 1.0683760683760708 " " absolute error = 2.4424906541753444000000000000000E-15 " " relative error = 2.2861712523081220000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.50000000000000500E-2 " " y[1] (analytic) = 1.0695187165775402 " " y[1] (numeric) = 1.0695187165775426 " " absolute error = 2.4424906541753444000000000000000E-15 " " relative error = 2.2837287616539470000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.60000000000000500E-2 " " y[1] (analytic) = 1.0706638115631693 " " y[1] (numeric) = 1.0706638115631717 " " absolute error = 2.4424906541753444000000000000000E-15 " " relative error = 2.28128627099977140000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.70000000000000500E-2 " " y[1] (analytic) = 1.0718113612004287 " " y[1] (numeric) = 1.0718113612004314 " " absolute error = 2.6645352591003757000000000000000E-15 " " relative error = 2.48601139674065020000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.80000000000000500E-2 " " y[1] (analytic) = 1.072961373390558 " " y[1] (numeric) = 1.0729613733905607 " " absolute error = 2.6645352591003757000000000000000E-15 " " relative error = 2.483346861481550000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 6.90000000000000500E-2 " " y[1] (analytic) = 1.0741138560687433 " " y[1] (numeric) = 1.0741138560687462 " " absolute error = 2.886579864025407000000000000000E-15 " " relative error = 2.6874058534076540000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.00000000000000500E-2 " " y[1] (analytic) = 1.0752688172043012 " " y[1] (numeric) = 1.0752688172043041 " " absolute error = 2.886579864025407000000000000000E-15 " " relative error = 2.6845192735436280000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.10000000000000500E-2 " " y[1] (analytic) = 1.0764262648008611 " " y[1] (numeric) = 1.0764262648008642 " " absolute error = 3.1086244689504383000000000000000E-15 " " relative error = 2.88791213165495800000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.20000000000000500E-2 " " y[1] (analytic) = 1.0775862068965518 " " y[1] (numeric) = 1.077586206896555 " " absolute error = 3.1086244689504383000000000000000E-15 " " relative error = 2.88480350718600600000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.30000000000000500E-2 " " y[1] (analytic) = 1.0787486515641855 " " y[1] (numeric) = 1.0787486515641889 " " absolute error = 3.3306690738754696000000000000000E-15 " " relative error = 3.08753023148256030000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.40000000000000500E-2 " " y[1] (analytic) = 1.0799136069114472 " " y[1] (numeric) = 1.0799136069114506 " " absolute error = 3.3306690738754696000000000000000E-15 " " relative error = 3.08419956240868430000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.50000000000000600E-2 " " y[1] (analytic) = 1.0810810810810811 " " y[1] (numeric) = 1.0810810810810847 " " absolute error = 3.552713678800501000000000000000E-15 " " relative error = 3.28626015289046340000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.60000000000000600E-2 " " y[1] (analytic) = 1.0822510822510822 " " y[1] (numeric) = 1.082251082251086 " " absolute error = 3.774758283725532000000000000000E-15 " " relative error = 3.48787665416239200000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.70000000000000600E-2 " " y[1] (analytic) = 1.0834236186348862 " " y[1] (numeric) = 1.08342361863489 " " absolute error = 3.774758283725532000000000000000E-15 " " relative error = 3.4841018958786660000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.80000000000000600E-2 " " y[1] (analytic) = 1.0845986984815619 " " y[1] (numeric) = 1.0845986984815656 " " absolute error = 3.774758283725532000000000000000E-15 " " relative error = 3.48032713759494000000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 7.90000000000000600E-2 " " y[1] (analytic) = 1.0857763300760044 " " y[1] (numeric) = 1.0857763300760082 " " absolute error = 3.774758283725532000000000000000E-15 " " relative error = 3.47655237931121500000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.00000000000000600E-2 " " y[1] (analytic) = 1.0869565217391306 " " y[1] (numeric) = 1.0869565217391344 " " absolute error = 3.774758283725532000000000000000E-15 " " relative error = 3.4727776210274890000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.10000000000000600E-2 " " y[1] (analytic) = 1.0881392818280742 " " y[1] (numeric) = 1.088139281828078 " " absolute error = 3.774758283725532000000000000000E-15 " " relative error = 3.46900286274376360000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.20000000000000600E-2 " " y[1] (analytic) = 1.0893246187363834 " " y[1] (numeric) = 1.0893246187363874 " " absolute error = 3.9968028886505635000000000000000E-15 " " relative error = 3.66906505178121730000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.30000000000000600E-2 " " y[1] (analytic) = 1.0905125408942205 " " y[1] (numeric) = 1.0905125408942242 " " absolute error = 3.774758283725532000000000000000E-15 " " relative error = 3.4614533461763125000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.40000000000000600E-2 " " y[1] (analytic) = 1.091703056768559 " " y[1] (numeric) = 1.091703056768563 " " absolute error = 3.9968028886505635000000000000000E-15 " " relative error = 3.6610714460039160000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.50000000000000600E-2 " " y[1] (analytic) = 1.092896174863388 " " y[1] (numeric) = 1.092896174863392 " " absolute error = 3.9968028886505635000000000000000E-15 " " relative error = 3.6570746431152650000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.60000000000000700E-2 " " y[1] (analytic) = 1.0940919037199126 " " y[1] (numeric) = 1.0940919037199166 " " absolute error = 3.9968028886505635000000000000000E-15 " " relative error = 3.65307784022661450000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.70000000000000700E-2 " " y[1] (analytic) = 1.095290251916758 " " y[1] (numeric) = 1.0952902519167622 " " absolute error = 4.218847493575595000000000000000E-15 " " relative error = 3.85180776163451750000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.80000000000000700E-2 " " y[1] (analytic) = 1.0964912280701755 " " y[1] (numeric) = 1.0964912280701797 " " absolute error = 4.218847493575595000000000000000E-15 " " relative error = 3.84758891414094250000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 8.90000000000000700E-2 " " y[1] (analytic) = 1.0976948408342482 " " y[1] (numeric) = 1.0976948408342524 " " absolute error = 4.218847493575595000000000000000E-15 " " relative error = 3.84337006664736640000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.00000000000000700E-2 " " y[1] (analytic) = 1.098901098901099 " " y[1] (numeric) = 1.0989010989011032 " " absolute error = 4.218847493575595000000000000000E-15 " " relative error = 3.8391512191537910000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.10000000000000700E-2 " " y[1] (analytic) = 1.1001100110011002 " " y[1] (numeric) = 1.1001100110011044 " " absolute error = 4.218847493575595000000000000000E-15 " " relative error = 3.83493237166021570000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.20000000000000700E-2 " " y[1] (analytic) = 1.1013215859030838 " " y[1] (numeric) = 1.101321585903088 " " absolute error = 4.218847493575595000000000000000E-15 " " relative error = 3.83071352416663960000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.30000000000000700E-2 " " y[1] (analytic) = 1.1025358324145536 " " y[1] (numeric) = 1.102535832414558 " " absolute error = 4.440892098500626000000000000000E-15 " " relative error = 4.02788913334006740000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.40000000000000700E-2 " " y[1] (analytic) = 1.1037527593818985 " " y[1] (numeric) = 1.103752759381903 " " absolute error = 4.440892098500626000000000000000E-15 " " relative error = 4.0234482412415670000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.50000000000000700E-2 " " y[1] (analytic) = 1.1049723756906078 " " y[1] (numeric) = 1.1049723756906122 " " absolute error = 4.440892098500626000000000000000E-15 " " relative error = 4.01900734914306670000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.60000000000000700E-2 " " y[1] (analytic) = 1.106194690265487 " " y[1] (numeric) = 1.1061946902654913 " " absolute error = 4.440892098500626000000000000000E-15 " " relative error = 4.01456645704456550000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.70000000000000800E-2 " " y[1] (analytic) = 1.1074197120708749 " " y[1] (numeric) = 1.1074197120708795 " " absolute error = 4.6629367034256575000000000000000E-15 " " relative error = 4.21063184319336870000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.80000000000000800E-2 " " y[1] (analytic) = 1.108647450110865 " " y[1] (numeric) = 1.1086474501108694 " " absolute error = 4.440892098500626000000000000000E-15 " " relative error = 4.0056846728475637000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 9.90000000000000800E-2 " " y[1] (analytic) = 1.109877913429523 " " y[1] (numeric) = 1.1098779134295276 " " absolute error = 4.6629367034256575000000000000000E-15 " " relative error = 4.20130596978651700000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10000000000000007 " " y[1] (analytic) = 1.1111111111111112 " " y[1] (numeric) = 1.111111111111116 " " absolute error = 4.884981308350689000000000000000E-15 " " relative error = 4.396483177515620000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10100000000000008 " " y[1] (analytic) = 1.1123470522803116 " " y[1] (numeric) = 1.1123470522803165 " " absolute error = 4.884981308350689000000000000000E-15 " " relative error = 4.39159819620726870000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10200000000000008 " " y[1] (analytic) = 1.11358574610245 " " y[1] (numeric) = 1.113585746102455 " " absolute error = 4.884981308350689000000000000000E-15 " " relative error = 4.3867132148989180000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10300000000000008 " " y[1] (analytic) = 1.1148272017837237 " " y[1] (numeric) = 1.1148272017837286 " " absolute error = 4.884981308350689000000000000000E-15 " " relative error = 4.3818282335905673000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10400000000000008 " " y[1] (analytic) = 1.1160714285714286 " " y[1] (numeric) = 1.1160714285714337 " " absolute error = 5.10702591327572000000000000000E-15 " " relative error = 4.5758952182950446000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10500000000000008 " " y[1] (analytic) = 1.11731843575419 " " y[1] (numeric) = 1.1173184357541952 " " absolute error = 5.10702591327572000000000000000E-15 " " relative error = 4.5707881923817684000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10600000000000008 " " y[1] (analytic) = 1.1185682326621924 " " y[1] (numeric) = 1.1185682326621977 " " absolute error = 5.329070518200751000000000000000E-15 " " relative error = 4.7641890432714720000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10700000000000008 " " y[1] (analytic) = 1.1198208286674134 " " y[1] (numeric) = 1.1198208286674187 " " absolute error = 5.329070518200751000000000000000E-15 " " relative error = 4.7588599727532704000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10800000000000008 " " y[1] (analytic) = 1.1210762331838566 " " y[1] (numeric) = 1.1210762331838622 " " absolute error = 5.551115123125783000000000000000E-15 " " relative error = 4.9515946898281970000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.10900000000000008 " " y[1] (analytic) = 1.1223344556677892 " " y[1] (numeric) = 1.1223344556677948 " " absolute error = 5.551115123125783000000000000000E-15 " " relative error = 4.9460435747050710000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11000000000000008 " " y[1] (analytic) = 1.1235955056179776 " " y[1] (numeric) = 1.1235955056179834 " " absolute error = 5.773159728050814000000000000000E-15 " " relative error = 5.1381121579652240000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11100000000000008 " " y[1] (analytic) = 1.124859392575928 " " y[1] (numeric) = 1.124859392575934 " " absolute error = 5.995204332975845000000000000000E-15 " " relative error = 5.3297366520155250000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11200000000000009 " " y[1] (analytic) = 1.1261261261261262 " " y[1] (numeric) = 1.1261261261261322 " " absolute error = 5.995204332975845000000000000000E-15 " " relative error = 5.323741447682551000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11300000000000009 " " y[1] (analytic) = 1.1273957158962797 " " y[1] (numeric) = 1.1273957158962857 " " absolute error = 5.995204332975845000000000000000E-15 " " relative error = 5.3177462433495740000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11400000000000009 " " y[1] (analytic) = 1.1286681715575622 " " y[1] (numeric) = 1.1286681715575684 " " absolute error = 6.217248937900877000000000000000E-15 " " relative error = 5.5084825589801760000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11500000000000009 " " y[1] (analytic) = 1.1299435028248588 " " y[1] (numeric) = 1.1299435028248652 " " absolute error = 6.439293542825908000000000000000E-15 " " relative error = 5.6987747854009290000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11600000000000009 " " y[1] (analytic) = 1.1312217194570138 " " y[1] (numeric) = 1.1312217194570202 " " absolute error = 6.439293542825908000000000000000E-15 " " relative error = 5.6923354918581020000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11700000000000009 " " y[1] (analytic) = 1.1325028312570782 " " y[1] (numeric) = 1.1325028312570848 " " absolute error = 6.661338147750939000000000000000E-15 " " relative error = 5.8819615844640790000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11800000000000009 " " y[1] (analytic) = 1.1337868480725626 " " y[1] (numeric) = 1.1337868480725692 " " absolute error = 6.661338147750939000000000000000E-15 " " relative error = 5.8753002463163270000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.11900000000000009 " " y[1] (analytic) = 1.135073779795687 " " y[1] (numeric) = 1.1350737797956938 " " absolute error = 6.8833827526759700000000000000000E-15 " " relative error = 6.0642602051075290000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.12000000000000009 " " y[1] (analytic) = 1.1363636363636365 " " y[1] (numeric) = 1.1363636363636436 " " absolute error = 7.105427357601002000000000000000E-15 " " relative error = 6.2527760746888800000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1210000000000001 " " y[1] (analytic) = 1.137656427758817 " " y[1] (numeric) = 1.1376564277588241 " " absolute error = 7.105427357601002000000000000000E-15 " " relative error = 6.245670647331280000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1220000000000001 " " y[1] (analytic) = 1.1389521640091118 " " y[1] (numeric) = 1.138952164009119 " " absolute error = 7.105427357601002000000000000000E-15 " " relative error = 6.2385652199736790000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1230000000000001 " " y[1] (analytic) = 1.1402508551881416 " " y[1] (numeric) = 1.1402508551881487 " " absolute error = 7.105427357601002000000000000000E-15 " " relative error = 6.2314597926160780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1240000000000001 " " y[1] (analytic) = 1.1415525114155252 " " y[1] (numeric) = 1.1415525114155325 " " absolute error = 7.327471962526033000000000000000E-15 " " relative error = 6.4188654391728040000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.12500000000000008 " " y[1] (analytic) = 1.142857142857143 " " y[1] (numeric) = 1.1428571428571503 " " absolute error = 7.327471962526033000000000000000E-15 " " relative error = 6.4115379672102780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.12600000000000008 " " y[1] (analytic) = 1.1441647597254005 " " y[1] (numeric) = 1.144164759725408 " " absolute error = 7.549516567451064000000000000000E-15 " " relative error = 6.598277479952230000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.12700000000000009 " " y[1] (analytic) = 1.1454753722794961 " " y[1] (numeric) = 1.1454753722795037 " " absolute error = 7.549516567451064000000000000000E-15 " " relative error = 6.5907279633847790000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.12800000000000009 " " y[1] (analytic) = 1.1467889908256883 " " y[1] (numeric) = 1.1467889908256959 " " absolute error = 7.549516567451064000000000000000E-15 " " relative error = 6.5831784468173270000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1290000000000001 " " y[1] (analytic) = 1.1481056257175661 " " y[1] (numeric) = 1.148105625717574 " " absolute error = 7.771561172376096000000000000000E-15 " " relative error = 6.7690297811395780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1300000000000001 " " y[1] (analytic) = 1.149425287356322 " " y[1] (numeric) = 1.1494252873563298 " " absolute error = 7.771561172376096000000000000000E-15 " " relative error = 6.7612582199672020000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1310000000000001 " " y[1] (analytic) = 1.1507479861910244 " " y[1] (numeric) = 1.1507479861910321 " " absolute error = 7.771561172376096000000000000000E-15 " " relative error = 6.7534866587948260000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1320000000000001 " " y[1] (analytic) = 1.1520737327188941 " " y[1] (numeric) = 1.1520737327189021 " " absolute error = 7.993605777301127000000000000000E-15 " " relative error = 6.9384498146973770000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1330000000000001 " " y[1] (analytic) = 1.1534025374855825 " " y[1] (numeric) = 1.1534025374855907 " " absolute error = 8.215650382226158000000000000000E-15 " " relative error = 7.1229688813900780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1340000000000001 " " y[1] (analytic) = 1.1547344110854505 " " y[1] (numeric) = 1.1547344110854587 " " absolute error = 8.215650382226158000000000000000E-15 " " relative error = 7.1147532310078520000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1350000000000001 " " y[1] (analytic) = 1.1560693641618498 " " y[1] (numeric) = 1.1560693641618582 " " absolute error = 8.43769498715119000000000000000E-15 " " relative error = 7.2986061638857780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1360000000000001 " " y[1] (analytic) = 1.1574074074074077 " " y[1] (numeric) = 1.157407407407416 " " absolute error = 8.43769498715119000000000000000E-15 " " relative error = 7.2901684688986270000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1370000000000001 " " y[1] (analytic) = 1.1587485515643108 " " y[1] (numeric) = 1.1587485515643194 " " absolute error = 8.659739592076221000000000000000E-15 " " relative error = 7.4733552679617780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1380000000000001 " " y[1] (analytic) = 1.1600928074245942 " " y[1] (numeric) = 1.160092807424603 " " absolute error = 8.881784197001252000000000000000E-15 " " relative error = 7.6560979778150780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1390000000000001 " " y[1] (analytic) = 1.16144018583043 " " y[1] (numeric) = 1.161440185830439 " " absolute error = 9.103828801926284000000000000000E-15 " " relative error = 7.8383965984585290000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1400000000000001 " " y[1] (analytic) = 1.1627906976744187 " " y[1] (numeric) = 1.162790697674428 " " absolute error = 9.325873406851315000000000000000E-15 " " relative error = 8.0202511298921310000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1410000000000001 " " y[1] (analytic) = 1.1641443538998837 " " y[1] (numeric) = 1.164144353899893 " " absolute error = 9.325873406851315000000000000000E-15 " " relative error = 8.0109252564852780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1420000000000001 " " y[1] (analytic) = 1.1655011655011658 " " y[1] (numeric) = 1.165501165501175 " " absolute error = 9.325873406851315000000000000000E-15 " " relative error = 8.0015993830784260000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1430000000000001 " " y[1] (analytic) = 1.166861143523921 " " y[1] (numeric) = 1.1668611435239304 " " absolute error = 9.547918011776346000000000000000E-15 " " relative error = 8.1825657360923280000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1440000000000001 " " y[1] (analytic) = 1.1682242990654208 " " y[1] (numeric) = 1.1682242990654306 " " absolute error = 9.769962616701378000000000000000E-15 " " relative error = 8.3630879998963770000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1450000000000001 " " y[1] (analytic) = 1.169590643274854 " " y[1] (numeric) = 1.169590643274864 " " absolute error = 9.992007221626409000000000000000E-15 " " relative error = 8.5431661744905780000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1460000000000001 " " y[1] (analytic) = 1.1709601873536302 " " y[1] (numeric) = 1.1709601873536402 " " absolute error = 9.992007221626409000000000000000E-15 " " relative error = 8.5331741672689520000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1470000000000001 " " y[1] (analytic) = 1.172332942555686 " " y[1] (numeric) = 1.172332942555696 " " absolute error = 9.992007221626409000000000000000E-15 " " relative error = 8.5231821600473250000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1480000000000001 " " y[1] (analytic) = 1.1737089201877937 " " y[1] (numeric) = 1.173708920187804 " " absolute error = 1.02140518265514400000000000000E-14 " " relative error = 8.7023721562218250000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1490000000000001 " " y[1] (analytic) = 1.175088131609871 " " y[1] (numeric) = 1.1750881316098813 " " absolute error = 1.043609643147647100000000000000E-14 " " relative error = 8.8811180631864760000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1500000000000001 " " y[1] (analytic) = 1.1764705882352944 " " y[1] (numeric) = 1.1764705882353048 " " absolute error = 1.043609643147647100000000000000E-14 " " relative error = 8.8706819667549990000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1510000000000001 " " y[1] (analytic) = 1.1778563015312133 " " y[1] (numeric) = 1.177856301531224 " " absolute error = 1.065814103640150300000000000000E-14 " " relative error = 9.0487617399048750000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1520000000000001 " " y[1] (analytic) = 1.1792452830188682 " " y[1] (numeric) = 1.1792452830188789 " " absolute error = 1.065814103640150300000000000000E-14 " " relative error = 9.0381035988684720000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1530000000000001 " " y[1] (analytic) = 1.180637544273908 " " y[1] (numeric) = 1.180637544273919 " " absolute error = 1.088018564132653400000000000000E-14 " " relative error = 9.2155172382035730000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1540000000000001 " " y[1] (analytic) = 1.182033096926714 " " y[1] (numeric) = 1.1820330969267252 " " absolute error = 1.110223024625156500000000000000E-14 " " relative error = 9.3924867883288220000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1550000000000001 " " y[1] (analytic) = 1.1834319526627222 " " y[1] (numeric) = 1.1834319526627333 " " absolute error = 1.110223024625156500000000000000E-14 " " relative error = 9.381384558082570000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1560000000000001 " " y[1] (analytic) = 1.184834123222749 " " y[1] (numeric) = 1.1848341232227604 " " absolute error = 1.132427485117659700000000000000E-14 " " relative error = 9.5576879743930440000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1570000000000001 " " y[1] (analytic) = 1.1862396204033216 " " y[1] (numeric) = 1.1862396204033332 " " absolute error = 1.154631945610162800000000000000E-14 " " relative error = 9.7335473014936710000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1580000000000001 " " y[1] (analytic) = 1.1876484560570073 " " y[1] (numeric) = 1.187648456057019 " " absolute error = 1.17683640610266600000000000000E-14 " " relative error = 9.9089625393844460000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.1590000000000001 " " y[1] (analytic) = 1.189060642092747 " " y[1] (numeric) = 1.1890606420927587 " " absolute error = 1.17683640610266600000000000000E-14 " " relative error = 9.8971941753234190000000000000E-13 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16000000000000011 " " y[1] (analytic) = 1.1904761904761907 " " y[1] (numeric) = 1.1904761904762027 " " absolute error = 1.19904086659516900000000000000E-14 " " relative error = 1.0071943279399419000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16100000000000012 " " y[1] (analytic) = 1.191895113230036 " " y[1] (numeric) = 1.1918951132300482 " " absolute error = 1.221245327087672200000000000000E-14 " " relative error = 1.0246248294265568000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16200000000000012 " " y[1] (analytic) = 1.1933174224343677 " " y[1] (numeric) = 1.19331742243438 " " absolute error = 1.243449787580175300000000000000E-14 " " relative error = 1.0420109219921868000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16300000000000012 " " y[1] (analytic) = 1.1947431302270013 " " y[1] (numeric) = 1.194743130227014 " " absolute error = 1.265654248072678500000000000000E-14 " " relative error = 1.0593526056368317000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16400000000000012 " " y[1] (analytic) = 1.196172248803828 " " y[1] (numeric) = 1.1961722488038407 " " absolute error = 1.265654248072678500000000000000E-14 " " relative error = 1.058086951388758900000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16500000000000012 " " y[1] (analytic) = 1.1976047904191618 " " y[1] (numeric) = 1.1976047904191747 " " absolute error = 1.287858708565181600000000000000E-14 " " relative error = 1.0753620216519265000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16600000000000012 " " y[1] (analytic) = 1.1990407673860914 " " y[1] (numeric) = 1.1990407673861043 " " absolute error = 1.287858708565181600000000000000E-14 " " relative error = 1.0740741629433612000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16700000000000012 " " y[1] (analytic) = 1.2004801920768309 " " y[1] (numeric) = 1.200480192076844 " " absolute error = 1.310063169057684700000000000000E-14 " " relative error = 1.0912826198250514000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16800000000000012 " " y[1] (analytic) = 1.201923076923077 " " y[1] (numeric) = 1.2019230769230904 " " absolute error = 1.332267629550187800000000000000E-14 " " relative error = 1.108446667785756200000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.16900000000000012 " " y[1] (analytic) = 1.2033694344163661 " " y[1] (numeric) = 1.2033694344163794 " " absolute error = 1.332267629550187800000000000000E-14 " " relative error = 1.107114400156206000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17000000000000012 " " y[1] (analytic) = 1.204819277108434 " " y[1] (numeric) = 1.2048192771084476 " " absolute error = 1.35447209004269100000000000000E-14 " " relative error = 1.1242118347354331000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17100000000000012 " " y[1] (analytic) = 1.2062726176115803 " " y[1] (numeric) = 1.2062726176115943 " " absolute error = 1.398881011027697200000000000000E-14 " " relative error = 1.159672358141961000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17200000000000013 " " y[1] (analytic) = 1.207729468599034 " " y[1] (numeric) = 1.207729468599048 " " absolute error = 1.398881011027697200000000000000E-14 " " relative error = 1.158273477130933000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17300000000000013 " " y[1] (analytic) = 1.2091898428053207 " " y[1] (numeric) = 1.209189842805335 " " absolute error = 1.421085471520200400000000000000E-14 " " relative error = 1.1752376849472053000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17400000000000013 " " y[1] (analytic) = 1.2106537530266346 " " y[1] (numeric) = 1.210653753026649 " " absolute error = 1.443289932012703500000000000000E-14 " " relative error = 1.1921574838424930000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17500000000000013 " " y[1] (analytic) = 1.2121212121212124 " " y[1] (numeric) = 1.212121212121227 " " absolute error = 1.465494392505206600000000000000E-14 " " relative error = 1.2090328738167953000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17600000000000013 " " y[1] (analytic) = 1.2135922330097089 " " y[1] (numeric) = 1.2135922330097237 " " absolute error = 1.487698852997709800000000000000E-14 " " relative error = 1.2258638548701127000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17700000000000013 " " y[1] (analytic) = 1.2150668286755775 " " y[1] (numeric) = 1.2150668286755923 " " absolute error = 1.487698852997709800000000000000E-14 " " relative error = 1.2243761560171149000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17800000000000013 " " y[1] (analytic) = 1.2165450121654504 " " y[1] (numeric) = 1.2165450121654655 " " absolute error = 1.50990331349021300000000000000E-14 " " relative error = 1.2411405236889547000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.17900000000000013 " " y[1] (analytic) = 1.2180267965895253 " " y[1] (numeric) = 1.2180267965895406 " " absolute error = 1.53210777398271600000000000000E-14 " " relative error = 1.2578604824398096000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18000000000000013 " " y[1] (analytic) = 1.2195121951219514 " " y[1] (numeric) = 1.219512195121967 " " absolute error = 1.55431223447521920000000000000E-14 " " relative error = 1.2745360322696794000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18100000000000013 " " y[1] (analytic) = 1.2210012210012213 " " y[1] (numeric) = 1.2210012210012369 " " absolute error = 1.55431223447521920000000000000E-14 " " relative error = 1.272981720035204200000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18200000000000013 " " y[1] (analytic) = 1.2224938875305627 " " y[1] (numeric) = 1.2224938875305784 " " absolute error = 1.576516694967722300000000000000E-14 " " relative error = 1.2895906564835966000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18300000000000013 " " y[1] (analytic) = 1.2239902080783356 " " y[1] (numeric) = 1.2239902080783516 " " absolute error = 1.598721155460225400000000000000E-14 " " relative error = 1.3061551840110040000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18400000000000014 " " y[1] (analytic) = 1.2254901960784317 " " y[1] (numeric) = 1.2254901960784477 " " absolute error = 1.598721155460225400000000000000E-14 " " relative error = 1.3045564628555437000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18500000000000014 " " y[1] (analytic) = 1.226993865030675 " " y[1] (numeric) = 1.2269938650306913 " " absolute error = 1.620925615952728500000000000000E-14 " " relative error = 1.3210543770014735000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18600000000000014 " " y[1] (analytic) = 1.2285012285012287 " " y[1] (numeric) = 1.228501228501245 " " absolute error = 1.643130076445231700000000000000E-14 " " relative error = 1.3375078822264186000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18700000000000014 " " y[1] (analytic) = 1.2300123001230014 " " y[1] (numeric) = 1.230012300123018 " " absolute error = 1.665334536937734800000000000000E-14 " " relative error = 1.3539169785303780000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18800000000000014 " " y[1] (analytic) = 1.2315270935960594 " " y[1] (numeric) = 1.2315270935960763 " " absolute error = 1.68753899743023800000000000000E-14 " " relative error = 1.370281665913353000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.18900000000000014 " " y[1] (analytic) = 1.2330456226880397 " " y[1] (numeric) = 1.2330456226880568 " " absolute error = 1.70974345792274100000000000000E-14 " " relative error = 1.3866019443753427000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19000000000000014 " " y[1] (analytic) = 1.234567901234568 " " y[1] (numeric) = 1.2345679012345854 " " absolute error = 1.731947918415244200000000000000E-14 " " relative error = 1.4028778139163475000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19100000000000014 " " y[1] (analytic) = 1.236093943139679 " " y[1] (numeric) = 1.2360939431396962 " " absolute error = 1.731947918415244200000000000000E-14 " " relative error = 1.4011458659979323000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19200000000000014 " " y[1] (analytic) = 1.2376237623762378 " " y[1] (numeric) = 1.2376237623762554 " " absolute error = 1.754152378907747300000000000000E-14 " " relative error = 1.4173551221574596000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19300000000000014 " " y[1] (analytic) = 1.2391573729863696 " " y[1] (numeric) = 1.2391573729863872 " " absolute error = 1.754152378907747300000000000000E-14 " " relative error = 1.4156009697785518000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19400000000000014 " " y[1] (analytic) = 1.240694789081886 " " y[1] (numeric) = 1.240694789081904 " " absolute error = 1.798561299892753600000000000000E-14 " " relative error = 1.4496404077135594000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19500000000000015 " " y[1] (analytic) = 1.2422360248447208 " " y[1] (numeric) = 1.2422360248447388 " " absolute error = 1.798561299892753600000000000000E-14 " " relative error = 1.4478418464136664000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19600000000000015 " " y[1] (analytic) = 1.2437810945273635 " " y[1] (numeric) = 1.2437810945273817 " " absolute error = 1.820765760385256700000000000000E-14 " " relative error = 1.463895671349746000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19700000000000015 " " y[1] (analytic) = 1.2453300124533004 " " y[1] (numeric) = 1.2453300124533189 " " absolute error = 1.8429702208777599000000000000E-14 " " relative error = 1.4799050873648410000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19800000000000015 " " y[1] (analytic) = 1.2468827930174566 " " y[1] (numeric) = 1.2468827930174753 " " absolute error = 1.86517468137026300000000000000E-14 " " relative error = 1.4958700944589506000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.19900000000000015 " " y[1] (analytic) = 1.248439450686642 " " y[1] (numeric) = 1.2484394506866607 " " absolute error = 1.86517468137026300000000000000E-14 " " relative error = 1.4940049197775800000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20000000000000015 " " y[1] (analytic) = 1.2500000000000002 " " y[1] (numeric) = 1.2500000000000193 " " absolute error = 1.909583602355269200000000000000E-14 " " relative error = 1.5276668818842150000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20100000000000015 " " y[1] (analytic) = 1.251564455569462 " " y[1] (numeric) = 1.2515644555694814 " " absolute error = 1.931788062847772400000000000000E-14 " " relative error = 1.5434986622153699000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20200000000000015 " " y[1] (analytic) = 1.2531328320802009 " " y[1] (numeric) = 1.2531328320802202 " " absolute error = 1.931788062847772400000000000000E-14 " " relative error = 1.5415668741525218000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20300000000000015 " " y[1] (analytic) = 1.254705144291092 " " y[1] (numeric) = 1.2547051442911115 " " absolute error = 1.953992523340275500000000000000E-14 " " relative error = 1.5573320411021993000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20400000000000015 " " y[1] (analytic) = 1.2562814070351762 " " y[1] (numeric) = 1.256281407035196 " " absolute error = 1.976196983832778600000000000000E-14 " " relative error = 1.5730527991308912000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20500000000000015 " " y[1] (analytic) = 1.2578616352201262 " " y[1] (numeric) = 1.2578616352201462 " " absolute error = 1.998401444325281800000000000000E-14 " " relative error = 1.5887291482385987000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20600000000000016 " " y[1] (analytic) = 1.2594458438287157 " " y[1] (numeric) = 1.259445843828736 " " absolute error = 2.02060590481778500000000000000E-14 " " relative error = 1.6043610884253207000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20700000000000016 " " y[1] (analytic) = 1.261034047919294 " " y[1] (numeric) = 1.2610340479193145 " " absolute error = 2.04281036531028800000000000000E-14 " " relative error = 1.619948619691058000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20800000000000016 " " y[1] (analytic) = 1.262626262626263 " " y[1] (numeric) = 1.2626262626262836 " " absolute error = 2.065014825802791200000000000000E-14 " " relative error = 1.63549174203581000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.20900000000000016 " " y[1] (analytic) = 1.2642225031605565 " " y[1] (numeric) = 1.2642225031605776 " " absolute error = 2.109423746787797400000000000000E-14 " " relative error = 1.6685541837091475000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21000000000000016 " " y[1] (analytic) = 1.2658227848101269 " " y[1] (numeric) = 1.2658227848101482 " " absolute error = 2.131628207280300600000000000000E-14 " " relative error = 1.683986283751437200000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21100000000000016 " " y[1] (analytic) = 1.2674271229404312 " " y[1] (numeric) = 1.2674271229404528 " " absolute error = 2.153832667772803700000000000000E-14 " " relative error = 1.6993739748727418000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21200000000000016 " " y[1] (analytic) = 1.269035532994924 " " y[1] (numeric) = 1.2690355329949459 " " absolute error = 2.176037128265306800000000000000E-14 " " relative error = 1.7147172570730615000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21300000000000016 " " y[1] (analytic) = 1.270648030495553 " " y[1] (numeric) = 1.270648030495575 " " absolute error = 2.1982415887578100000000000000E-14 " " relative error = 1.730016130352396200000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21400000000000016 " " y[1] (analytic) = 1.2722646310432573 " " y[1] (numeric) = 1.2722646310432795 " " absolute error = 2.22044604925031300000000000000E-14 " " relative error = 1.7452705947107458000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21500000000000016 " " y[1] (analytic) = 1.2738853503184717 " " y[1] (numeric) = 1.2738853503184941 " " absolute error = 2.242650509742816200000000000000E-14 " " relative error = 1.7604806501481102000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21600000000000016 " " y[1] (analytic) = 1.275510204081633 " " y[1] (numeric) = 1.2755102040816557 " " absolute error = 2.264854970235319300000000000000E-14 " " relative error = 1.7756462966644898000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21700000000000016 " " y[1] (analytic) = 1.2771392081736912 " " y[1] (numeric) = 1.2771392081737143 " " absolute error = 2.309263891220325600000000000000E-14 " " relative error = 1.8081536268255144000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21800000000000017 " " y[1] (analytic) = 1.2787723785166243 " " y[1] (numeric) = 1.2787723785166476 " " absolute error = 2.331468351712828700000000000000E-14 " " relative error = 1.8232082510394318000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.21900000000000017 " " y[1] (analytic) = 1.2804097311139568 " " y[1] (numeric) = 1.2804097311139804 " " absolute error = 2.35367281220533200000000000000E-14 " " relative error = 1.8382184663323634000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22000000000000017 " " y[1] (analytic) = 1.2820512820512824 " " y[1] (numeric) = 1.2820512820513061 " " absolute error = 2.37587727269783500000000000000E-14 " " relative error = 1.853184272704311000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22100000000000017 " " y[1] (analytic) = 1.283697047496791 " " y[1] (numeric) = 1.283697047496815 " " absolute error = 2.39808173319033800000000000000E-14 " " relative error = 1.868105670155273000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22200000000000017 " " y[1] (analytic) = 1.2853470437017998 " " y[1] (numeric) = 1.285347043701824 " " absolute error = 2.420286193682841300000000000000E-14 " " relative error = 1.8829826586852502000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22300000000000017 " " y[1] (analytic) = 1.2870012870012872 " " y[1] (numeric) = 1.2870012870013117 " " absolute error = 2.442490654175344400000000000000E-14 " " relative error = 1.897815238294242300000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22400000000000017 " " y[1] (analytic) = 1.2886597938144333 " " y[1] (numeric) = 1.288659793814458 " " absolute error = 2.464695114667847500000000000000E-14 " " relative error = 1.9126034089822494000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22500000000000017 " " y[1] (analytic) = 1.2903225806451617 " " y[1] (numeric) = 1.2903225806451866 " " absolute error = 2.486899575160350700000000000000E-14 " " relative error = 1.927347170749271000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22600000000000017 " " y[1] (analytic) = 1.2919896640826878 " " y[1] (numeric) = 1.2919896640827129 " " absolute error = 2.509104035652854000000000000000E-14 " " relative error = 1.942046523595308300000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22700000000000017 " " y[1] (analytic) = 1.2936610608020702 " " y[1] (numeric) = 1.2936610608020958 " " absolute error = 2.5535129566378600000000000000E-14 " " relative error = 1.9738655154810655000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22800000000000017 " " y[1] (analytic) = 1.2953367875647672 " " y[1] (numeric) = 1.295336787564793 " " absolute error = 2.575717417130363000000000000000E-14 " " relative error = 1.9884538460246398000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.22900000000000018 " " y[1] (analytic) = 1.2970168612191961 " " y[1] (numeric) = 1.2970168612192223 " " absolute error = 2.620126338115369400000000000000E-14 " " relative error = 2.0201174066869493000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23000000000000018 " " y[1] (analytic) = 1.2987012987012991 " " y[1] (numeric) = 1.2987012987013256 " " absolute error = 2.642330798607872600000000000000E-14 " " relative error = 2.034594714928061000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23100000000000018 " " y[1] (analytic) = 1.3003901170351109 " " y[1] (numeric) = 1.3003901170351377 " " absolute error = 2.686739719592879000000000000000E-14 " " relative error = 2.0661028443669233000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23200000000000018 " " y[1] (analytic) = 1.3020833333333337 " " y[1] (numeric) = 1.3020833333333608 " " absolute error = 2.70894418008538200000000000000E-14 " " relative error = 2.0804691303055725000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23300000000000018 " " y[1] (analytic) = 1.3037809647979144 " " y[1] (numeric) = 1.3037809647979417 " " absolute error = 2.73114864057788500000000000000E-14 " " relative error = 2.094791007323237000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23400000000000018 " " y[1] (analytic) = 1.305483028720627 " " y[1] (numeric) = 1.3054830287206547 " " absolute error = 2.775557561562891400000000000000E-14 " " relative error = 2.1260770921571745000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23500000000000018 " " y[1] (analytic) = 1.3071895424836606 " " y[1] (numeric) = 1.3071895424836886 " " absolute error = 2.797762022055394500000000000000E-14 " " relative error = 2.140287946872376000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23600000000000018 " " y[1] (analytic) = 1.3089005235602098 " " y[1] (numeric) = 1.3089005235602382 " " absolute error = 2.84217094304040100000000000000E-14 " " relative error = 2.1714186004828656000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23700000000000018 " " y[1] (analytic) = 1.3106159895150724 " " y[1] (numeric) = 1.3106159895151013 " " absolute error = 2.88657986402540700000000000000E-14 " " relative error = 2.202460436251385300000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23800000000000018 " " y[1] (analytic) = 1.3123359580052496 " " y[1] (numeric) = 1.312335958005279 " " absolute error = 2.93098878501041300000000000000E-14 " " relative error = 2.2334134541779344000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.23900000000000018 " " y[1] (analytic) = 1.3140604467805523 " " y[1] (numeric) = 1.3140604467805819 " " absolute error = 2.953193245502916400000000000000E-14 " " relative error = 2.247380059827718800000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.24000000000000019 " " y[1] (analytic) = 1.3157894736842108 " " y[1] (numeric) = 1.3157894736842408 " " absolute error = 2.997602166487922700000000000000E-14 " " relative error = 2.278177646530821000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2410000000000002 " " y[1] (analytic) = 1.317523056653492 " " y[1] (numeric) = 1.317523056653522 " " absolute error = 3.01980662698042600000000000000E-14 " " relative error = 2.2920332298781423000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2420000000000002 " " y[1] (analytic) = 1.319261213720317 " " y[1] (numeric) = 1.3192612137203477 " " absolute error = 3.06421554796543200000000000000E-14 " " relative error = 2.3226753853577967000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2430000000000002 " " y[1] (analytic) = 1.3210039630118895 " " y[1] (numeric) = 1.3210039630119204 " " absolute error = 3.08642000845793500000000000000E-14 " " relative error = 2.336419946402656000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2440000000000002 " " y[1] (analytic) = 1.3227513227513232 " " y[1] (numeric) = 1.3227513227513543 " " absolute error = 3.10862446895043830000000000000E-14 " " relative error = 2.3501200985265305000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2450000000000002 " " y[1] (analytic) = 1.3245033112582785 " " y[1] (numeric) = 1.32450331125831 " " absolute error = 3.153033389935444600000000000000E-14 " " relative error = 2.3805402094012604000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2460000000000002 " " y[1] (analytic) = 1.3262599469496026 " " y[1] (numeric) = 1.3262599469496343 " " absolute error = 3.17523785042794770000000000000E-14 " " relative error = 2.394129339222672000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2470000000000002 " " y[1] (analytic) = 1.3280212483399738 " " y[1] (numeric) = 1.328021248340006 " " absolute error = 3.21964677141295400000000000000E-14 " " relative error = 2.4243940188739535000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2480000000000002 " " y[1] (analytic) = 1.3297872340425536 " " y[1] (numeric) = 1.329787234042586 " " absolute error = 3.24185123190545700000000000000E-14 " " relative error = 2.437872126392903200000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2490000000000002 " " y[1] (analytic) = 1.3315579227696408 " " y[1] (numeric) = 1.3315579227696737 " " absolute error = 3.286260152890463400000000000000E-14 " " relative error = 2.4679813748207374000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.25000000000000017 " " y[1] (analytic) = 1.3333333333333337 " " y[1] (numeric) = 1.3333333333333668 " " absolute error = 3.308464613382966500000000000000E-14 " " relative error = 2.481348460037224000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.25100000000000017 " " y[1] (analytic) = 1.3351134846461952 " " y[1] (numeric) = 1.3351134846462287 " " absolute error = 3.35287353436797300000000000000E-14 " " relative error = 2.511302277241611000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.25200000000000017 " " y[1] (analytic) = 1.3368983957219256 " " y[1] (numeric) = 1.3368983957219593 " " absolute error = 3.37507799486047600000000000000E-14 " " relative error = 2.5245583401556354000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.25300000000000017 " " y[1] (analytic) = 1.3386880856760377 " " y[1] (numeric) = 1.3386880856760721 " " absolute error = 3.44169137633798500000000000000E-14 " " relative error = 2.5709434581244744000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.25400000000000017 " " y[1] (analytic) = 1.340482573726542 " " y[1] (numeric) = 1.3404825737265766 " " absolute error = 3.463895836830488400000000000000E-14 " " relative error = 2.5840662942755440000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.25500000000000017 " " y[1] (analytic) = 1.3422818791946312 " " y[1] (numeric) = 1.3422818791946662 " " absolute error = 3.508304757815494700000000000000E-14 " " relative error = 2.613687044572543000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.25600000000000017 " " y[1] (analytic) = 1.3440860215053767 " " y[1] (numeric) = 1.344086021505412 " " absolute error = 3.53050921830799800000000000000E-14 " " relative error = 2.6266988584211500000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2570000000000002 " " y[1] (analytic) = 1.3458950201884254 " " y[1] (numeric) = 1.3458950201884614 " " absolute error = 3.59712259978550700000000000000E-14 " " relative error = 2.6726620916406313000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2580000000000002 " " y[1] (analytic) = 1.3477088948787066 " " y[1] (numeric) = 1.3477088948787428 " " absolute error = 3.619327060278010300000000000000E-14 " " relative error = 2.685540678726283000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2590000000000002 " " y[1] (analytic) = 1.3495276653171393 " " y[1] (numeric) = 1.349527665317176 " " absolute error = 3.663735981263016600000000000000E-14 " " relative error = 2.714828362115895000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2600000000000002 " " y[1] (analytic) = 1.3513513513513518 " " y[1] (numeric) = 1.3513513513513886 " " absolute error = 3.6859404417555197000000000000E-14 " " relative error = 2.7275959268990840000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2610000000000002 " " y[1] (analytic) = 1.3531799729364007 " " y[1] (numeric) = 1.3531799729364382 " " absolute error = 3.75255382323302900000000000000E-14 " " relative error = 2.773137275369208000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2620000000000002 " " y[1] (analytic) = 1.3550135501355018 " " y[1] (numeric) = 1.3550135501355394 " " absolute error = 3.75255382323302900000000000000E-14 " " relative error = 2.7693847215459744000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2630000000000002 " " y[1] (analytic) = 1.35685210312076 " " y[1] (numeric) = 1.3568521031207983 " " absolute error = 3.819167204710538500000000000000E-14 " " relative error = 2.8147262298716663000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2640000000000002 " " y[1] (analytic) = 1.3586956521739135 " " y[1] (numeric) = 1.358695652173952 " " absolute error = 3.841371665203041600000000000000E-14 " " relative error = 2.8272495455894375000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2650000000000002 " " y[1] (analytic) = 1.360544217687075 " " y[1] (numeric) = 1.360544217687114 " " absolute error = 3.90798504668055100000000000000E-14 " " relative error = 2.8723690093102044000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2660000000000002 " " y[1] (analytic) = 1.362397820163488 " " y[1] (numeric) = 1.3623978201635274 " " absolute error = 3.93018950717305400000000000000E-14 " " relative error = 2.884759098265021000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2670000000000002 " " y[1] (analytic) = 1.3642564802182813 " " y[1] (numeric) = 1.364256480218321 " " absolute error = 3.974598428158060400000000000000E-14 " " relative error = 2.9133806478398580000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2680000000000002 " " y[1] (analytic) = 1.3661202185792354 " " y[1] (numeric) = 1.3661202185792753 " " absolute error = 3.996802888650563500000000000000E-14 " " relative error = 2.9256597144922120000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2690000000000002 " " y[1] (analytic) = 1.3679890560875516 " " y[1] (numeric) = 1.367989056087592 " " absolute error = 4.0412118096355700000000000000E-14 " " relative error = 2.954125832843601000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2700000000000002 " " y[1] (analytic) = 1.3698630136986305 " " y[1] (numeric) = 1.3698630136986714 " " absolute error = 4.08562073062057600000000000000E-14 " " relative error = 2.9825031333530200000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2710000000000002 " " y[1] (analytic) = 1.3717421124828535 " " y[1] (numeric) = 1.3717421124828948 " " absolute error = 4.130029651605582300000000000000E-14 " " relative error = 3.010791616020469000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2720000000000002 " " y[1] (analytic) = 1.373626373626374 " " y[1] (numeric) = 1.3736263736264156 " " absolute error = 4.152234112098085500000000000000E-14 " " relative error = 3.022826433607405000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2730000000000002 " " y[1] (analytic) = 1.3755158184319123 " " y[1] (numeric) = 1.3755158184319545 " " absolute error = 4.21884749357559500000000000000E-14 " " relative error = 3.0671021278294570000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2740000000000002 " " y[1] (analytic) = 1.3774104683195596 " " y[1] (numeric) = 1.3774104683196022 " " absolute error = 4.26325641456060100000000000000E-14 " " relative error = 3.0951241569709960000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2750000000000002 " " y[1] (analytic) = 1.3793103448275865 " " y[1] (numeric) = 1.3793103448276296 " " absolute error = 4.307665335545607400000000000000E-14 " " relative error = 3.1230573682705650000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2760000000000002 " " y[1] (analytic) = 1.3812154696132601 " " y[1] (numeric) = 1.3812154696133034 " " absolute error = 4.329869796038110500000000000000E-14 " " relative error = 3.1348257323315910000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2770000000000002 " " y[1] (analytic) = 1.3831258644536655 " " y[1] (numeric) = 1.3831258644537094 " " absolute error = 4.3964831775156200000000000000E-14 " " relative error = 3.1786573373437926000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2780000000000002 " " y[1] (analytic) = 1.385041551246538 " " y[1] (numeric) = 1.3850415512465821 " " absolute error = 4.41868763800812300000000000000E-14 " " relative error = 3.1902924746418637000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2790000000000002 " " y[1] (analytic) = 1.386962552011096 " " y[1] (numeric) = 1.3869625520111408 " " absolute error = 4.485301019485632400000000000000E-14 " " relative error = 3.2339020350491404000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2800000000000002 " " y[1] (analytic) = 1.3888888888888893 " " y[1] (numeric) = 1.3888888888889346 " " absolute error = 4.52970994047063870000000000000E-14 " " relative error = 3.261391157138859300000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2810000000000002 " " y[1] (analytic) = 1.3908205841446457 " " y[1] (numeric) = 1.3908205841446915 " " absolute error = 4.57411886145564500000000000000E-14 " " relative error = 3.288791461386607600000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2820000000000002 " " y[1] (analytic) = 1.3927576601671314 " " y[1] (numeric) = 1.3927576601671776 " " absolute error = 4.61852778244065100000000000000E-14 " " relative error = 3.3161029477923865000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2830000000000002 " " y[1] (analytic) = 1.3947001394700143 " " y[1] (numeric) = 1.394700139470061 " " absolute error = 4.685141163918160600000000000000E-14 " " relative error = 3.3592462145293206000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2840000000000002 " " y[1] (analytic) = 1.3966480446927378 " " y[1] (numeric) = 1.3966480446927851 " " absolute error = 4.72955008490316700000000000000E-14 " " relative error = 3.3863578607906664000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2850000000000002 " " y[1] (analytic) = 1.398601398601399 " " y[1] (numeric) = 1.3986013986014467 " " absolute error = 4.77395900588817300000000000000E-14 " " relative error = 3.413380689210043000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2860000000000002 " " y[1] (analytic) = 1.4005602240896364 " " y[1] (numeric) = 1.4005602240896844 " " absolute error = 4.79616346638067600000000000000E-14 " " relative error = 3.4244607149958020000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2870000000000002 " " y[1] (analytic) = 1.4025245441795233 " " y[1] (numeric) = 1.4025245441795722 " " absolute error = 4.88498130835068900000000000000E-14 " " relative error = 3.4829916728540405000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2880000000000002 " " y[1] (analytic) = 1.4044943820224725 " " y[1] (numeric) = 1.4044943820225215 " " absolute error = 4.90718576884319200000000000000E-14 " " relative error = 3.493916267416351000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2890000000000002 " " y[1] (analytic) = 1.4064697609001409 " " y[1] (numeric) = 1.4064697609001908 " " absolute error = 4.996003610813204400000000000000E-14 " " relative error = 3.552158567288188000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2900000000000002 " " y[1] (analytic) = 1.4084507042253527 " " y[1] (numeric) = 1.4084507042254029 " " absolute error = 5.018208071305708000000000000000E-14 " " relative error = 3.5629277306270507000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2910000000000002 " " y[1] (analytic) = 1.4104372355430186 " " y[1] (numeric) = 1.4104372355430697 " " absolute error = 5.1070259132757200000000000000E-14 " " relative error = 3.6208813725124850000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2920000000000002 " " y[1] (analytic) = 1.412429378531074 " " y[1] (numeric) = 1.4124293785311253 " " absolute error = 5.12923037376822300000000000000E-14 " " relative error = 3.6314951046279010000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2930000000000002 " " y[1] (analytic) = 1.4144271570014146 " " y[1] (numeric) = 1.4144271570014668 " " absolute error = 5.21804821573823600000000000000E-14 " " relative error = 3.689160088526932000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2940000000000002 " " y[1] (analytic) = 1.4164305949008504 " " y[1] (numeric) = 1.4164305949009028 " " absolute error = 5.24025267623073900000000000000E-14 " " relative error = 3.6996183894189005000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2950000000000002 " " y[1] (analytic) = 1.418439716312057 " " y[1] (numeric) = 1.4184397163121103 " " absolute error = 5.329070518200751000000000000000E-14 " " relative error = 3.756994715331529000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2960000000000002 " " y[1] (analytic) = 1.420454545454546 " " y[1] (numeric) = 1.4204545454545996 " " absolute error = 5.351274978693255000000000000000E-14 " " relative error = 3.76729758500005000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2970000000000002 " " y[1] (analytic) = 1.4224751066856334 " " y[1] (numeric) = 1.4224751066856878 " " absolute error = 5.44009282066326700000000000000E-14 " " relative error = 3.824385252926276000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2980000000000002 " " y[1] (analytic) = 1.4245014245014251 " " y[1] (numeric) = 1.4245014245014798 " " absolute error = 5.4622972811557700000000000000E-14 " " relative error = 3.834532691371349000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.2990000000000002 " " y[1] (analytic) = 1.4265335235378034 " " y[1] (numeric) = 1.426533523537859 " " absolute error = 5.55111512312578300000000000000E-14 " " relative error = 3.891331701311173000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3000000000000002 " " y[1] (analytic) = 1.428571428571429 " " y[1] (numeric) = 1.428571428571485 " " absolute error = 5.59552404411078900000000000000E-14 " " relative error = 3.916866830877550600000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3010000000000002 " " y[1] (analytic) = 1.4306151645207443 " " y[1] (numeric) = 1.430615164520801 " " absolute error = 5.662137425588298000000000000000E-14 " " relative error = 3.95783406048622000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3020000000000002 " " y[1] (analytic) = 1.4326647564469919 " " y[1] (numeric) = 1.432664756447049 " " absolute error = 5.70654634657330500000000000000E-14 " " relative error = 3.9831693499081655000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3030000000000002 " " y[1] (analytic) = 1.434720229555237 " " y[1] (numeric) = 1.4347202295552948 " " absolute error = 5.77315972805081400000000000000E-14 " " relative error = 4.0238923304514170000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3040000000000002 " " y[1] (analytic) = 1.4367816091954029 " " y[1] (numeric) = 1.436781609195461 " " absolute error = 5.8175686490358200000000000000E-14 " " relative error = 4.0490277797289290000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3050000000000002 " " y[1] (analytic) = 1.4388489208633097 " " y[1] (numeric) = 1.4388489208633688 " " absolute error = 5.90638649100583300000000000000E-14 " " relative error = 4.104938611249053000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3060000000000002 " " y[1] (analytic) = 1.4409221902017297 " " y[1] (numeric) = 1.4409221902017892 " " absolute error = 5.95079541199083900000000000000E-14 " " relative error = 4.129852015921640600000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3070000000000002 " " y[1] (analytic) = 1.4430014430014433 " " y[1] (numeric) = 1.4430014430015037 " " absolute error = 6.03961325396085200000000000000E-14 " " relative error = 4.185451984994869000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3080000000000002 " " y[1] (analytic) = 1.4450867052023126 " " y[1] (numeric) = 1.4450867052023735 " " absolute error = 6.08402217494585800000000000000E-14 " " relative error = 4.210143345062532000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3090000000000002 " " y[1] (analytic) = 1.4471780028943564 " " y[1] (numeric) = 1.447178002894418 " " absolute error = 6.15063555642336700000000000000E-14 " " relative error = 4.250089169488545000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3100000000000002 " " y[1] (analytic) = 1.4492753623188412 " " y[1] (numeric) = 1.4492753623189032 " " absolute error = 6.19504447740837300000000000000E-14 " " relative error = 4.274580689411776000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3110000000000002 " " y[1] (analytic) = 1.4513788098693763 " " y[1] (numeric) = 1.4513788098694391 " " absolute error = 6.28386231937838600000000000000E-14 " " relative error = 4.3295811380517063000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3120000000000002 " " y[1] (analytic) = 1.4534883720930238 " " y[1] (numeric) = 1.4534883720930871 " " absolute error = 6.32827124036339200000000000000E-14 " " relative error = 4.353850613370012000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3130000000000002 " " y[1] (analytic) = 1.4556040756914124 " " y[1] (numeric) = 1.4556040756914765 " " absolute error = 6.41708908233340500000000000000E-14 " " relative error = 4.408540199563048000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3140000000000002 " " y[1] (analytic) = 1.4577259475218665 " " y[1] (numeric) = 1.457725947521931 " " absolute error = 6.46149800331841100000000000000E-14 " " relative error = 4.432587630276428300000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3150000000000002 " " y[1] (analytic) = 1.4598540145985406 " " y[1] (numeric) = 1.459854014598606 " " absolute error = 6.55031584528842400000000000000E-14 " " relative error = 4.4869663540225685000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3160000000000002 " " y[1] (analytic) = 1.461988304093568 " " y[1] (numeric) = 1.461988304093634 " " absolute error = 6.5947247662734300000000000000E-14 " " relative error = 4.510791740131024400000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3170000000000002 " " y[1] (analytic) = 1.464128843338214 " " y[1] (numeric) = 1.464128843338281 " " absolute error = 6.70574706873594600000000000000E-14 " " relative error = 4.58002524794664970000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3180000000000002 " " y[1] (analytic) = 1.4662756598240476 " " y[1] (numeric) = 1.466275659824115 " " absolute error = 6.72795152922844900000000000000E-14 " " relative error = 4.5884629429338003000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.31900000000000023 " " y[1] (analytic) = 1.468428781204112 " " y[1] (numeric) = 1.4684287812041803 " " absolute error = 6.83897383169096400000000000000E-14 " " relative error = 4.657341179381545600000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32000000000000023 " " y[1] (analytic) = 1.4705882352941182 " " y[1] (numeric) = 1.470588235294187 " " absolute error = 6.8833827526759700000000000000E-14 " " relative error = 4.680700271819658300000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32100000000000023 " " y[1] (analytic) = 1.472754050073638 " " y[1] (numeric) = 1.4727540500737077 " " absolute error = 6.97220059464598300000000000000E-14 " " relative error = 4.734124203764621400000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32200000000000023 " " y[1] (analytic) = 1.4749262536873162 " " y[1] (numeric) = 1.4749262536873864 " " absolute error = 7.01660951563098900000000000000E-14 " " relative error = 4.757261251597809000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32300000000000023 " " y[1] (analytic) = 1.477104874446086 " " y[1] (numeric) = 1.4771048744461572 " " absolute error = 7.12763181809350500000000000000E-14 " " relative error = 4.825406740849302300000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32400000000000023 " " y[1] (analytic) = 1.479289940828403 " " y[1] (numeric) = 1.4792899408284748 " " absolute error = 7.17204073907851100000000000000E-14 " " relative error = 4.848299539617071000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32500000000000023 " " y[1] (analytic) = 1.4814814814814818 " " y[1] (numeric) = 1.4814814814815547 " " absolute error = 7.28306304154102700000000000000E-14 " " relative error = 4.916067553040192600000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32600000000000023 " " y[1] (analytic) = 1.4836795252225525 " " y[1] (numeric) = 1.4836795252226258 " " absolute error = 7.32747196252603300000000000000E-14 " " relative error = 4.938716102742544700000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32700000000000023 " " y[1] (analytic) = 1.4858841010401194 " " y[1] (numeric) = 1.4858841010401935 " " absolute error = 7.41628980449604600000000000000E-14 " " relative error = 4.991163038425837000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32800000000000024 " " y[1] (analytic) = 1.4880952380952388 " " y[1] (numeric) = 1.4880952380953136 " " absolute error = 7.48290318597355500000000000000E-14 " " relative error = 5.028510940974227000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.32900000000000024 " " y[1] (analytic) = 1.4903129657228023 " " y[1] (numeric) = 1.4903129657228782 " " absolute error = 7.59392548843607100000000000000E-14 " " relative error = 5.095524002740602000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33000000000000024 " " y[1] (analytic) = 1.4925373134328366 " " y[1] (numeric) = 1.492537313432913 " " absolute error = 7.63833440942107700000000000000E-14 " " relative error = 5.117684054312119000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33100000000000024 " " y[1] (analytic) = 1.494768310911809 " " y[1] (numeric) = 1.4947683109118866 " " absolute error = 7.74935671188359300000000000000E-14 " " relative error = 5.184319640250122000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33200000000000024 " " y[1] (analytic) = 1.4970059880239528 " " y[1] (numeric) = 1.4970059880240307 " " absolute error = 7.79376563286859900000000000000E-14 " " relative error = 5.2062354427562220000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33300000000000024 " " y[1] (analytic) = 1.499250374812594 " " y[1] (numeric) = 1.499250374812673 " " absolute error = 7.90478793533111500000000000000E-14 " " relative error = 5.272493552865852000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33400000000000024 " " y[1] (analytic) = 1.501501501501502 " " y[1] (numeric) = 1.5015015015015818 " " absolute error = 7.97140131680862400000000000000E-14 " " relative error = 5.3089532769945410000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33500000000000024 " " y[1] (analytic) = 1.503759398496241 " " y[1] (numeric) = 1.5037593984963218 " " absolute error = 8.0824236192711400000000000000E-14 " " relative error = 5.374811706815307000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33600000000000024 " " y[1] (analytic) = 1.506024096385543 " " y[1] (numeric) = 1.5060240963856242 " " absolute error = 8.12683254025614600000000000000E-14 " " relative error = 5.396216806730078000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33700000000000024 " " y[1] (analytic) = 1.5082956259426852 " " y[1] (numeric) = 1.5082956259427676 " " absolute error = 8.23785484271866200000000000000E-14 " " relative error = 5.461697760722470000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33800000000000024 " " y[1] (analytic) = 1.510574018126889 " " y[1] (numeric) = 1.510574018126972 " " absolute error = 8.30446822419617100000000000000E-14 " " relative error = 5.4975579644178630000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.33900000000000025 " " y[1] (analytic) = 1.5128593040847205 " " y[1] (numeric) = 1.5128593040848048 " " absolute error = 8.4376949871511900000000000000E-14 " " relative error = 5.577316386506935000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34000000000000025 " " y[1] (analytic) = 1.5151515151515158 " " y[1] (numeric) = 1.5151515151516008 " " absolute error = 8.50430836862869900000000000000E-14 " " relative error = 5.612843523294939000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34100000000000025 " " y[1] (analytic) = 1.5174506828528078 " " y[1] (numeric) = 1.517450682852894 " " absolute error = 8.61533067109121500000000000000E-14 " " relative error = 5.677502912249108000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34200000000000025 " " y[1] (analytic) = 1.5197568389057758 " " y[1] (numeric) = 1.5197568389058627 " " absolute error = 8.68194405256872400000000000000E-14 " " relative error = 5.712719186590218000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34300000000000025 " " y[1] (analytic) = 1.5220700152207005 " " y[1] (numeric) = 1.5220700152207887 " " absolute error = 8.81517081552374300000000000000E-14 " " relative error = 5.791567225799098000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34400000000000025 " " y[1] (analytic) = 1.5243902439024397 " " y[1] (numeric) = 1.5243902439025285 " " absolute error = 8.88178419700125200000000000000E-14 " " relative error = 5.826450433232819000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34500000000000025 " " y[1] (analytic) = 1.526717557251909 " " y[1] (numeric) = 1.5267175572519989 " " absolute error = 8.99280649946376800000000000000E-14 " " relative error = 5.890288257148766000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34600000000000025 " " y[1] (analytic) = 1.5290519877675848 " " y[1] (numeric) = 1.5290519877676756 " " absolute error = 9.0816243414337800000000000000E-14 " " relative error = 5.93938231929769000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34700000000000025 " " y[1] (analytic) = 1.5313935681470143 " " y[1] (numeric) = 1.5313935681471063 " " absolute error = 9.19264664389629600000000000000E-14 " " relative error = 6.002798258464279000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34800000000000025 " " y[1] (analytic) = 1.5337423312883443 " " y[1] (numeric) = 1.533742331288437 " " absolute error = 9.28146448586630900000000000000E-14 " " relative error = 6.0515148447848310000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.34900000000000025 " " y[1] (analytic) = 1.536098310291859 " " y[1] (numeric) = 1.5360983102919532 " " absolute error = 9.41469124882132700000000000000E-14 " " relative error = 6.128964002982683000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35000000000000026 " " y[1] (analytic) = 1.5384615384615392 " " y[1] (numeric) = 1.538461538461634 " " absolute error = 9.48130463029883700000000000000E-14 " " relative error = 6.162848009694241000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35100000000000026 " " y[1] (analytic) = 1.540832049306626 " " y[1] (numeric) = 1.5408320493067222 " " absolute error = 9.61453139325385600000000000000E-14 " " relative error = 6.23983087422175000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35200000000000026 " " y[1] (analytic) = 1.5432098765432105 " " y[1] (numeric) = 1.5432098765433075 " " absolute error = 9.70334923522386800000000000000E-14 " " relative error = 6.287770304425064000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35300000000000026 " " y[1] (analytic) = 1.5455950540958274 " " y[1] (numeric) = 1.5455950540959256 " " absolute error = 9.81437153768638400000000000000E-14 " " relative error = 6.349898384883088000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35400000000000026 " " y[1] (analytic) = 1.547987616099072 " " y[1] (numeric) = 1.547987616099171 " " absolute error = 9.90318937965639600000000000000E-14 " " relative error = 6.3974603392580290000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35500000000000026 " " y[1] (analytic) = 1.5503875968992253 " " y[1] (numeric) = 1.5503875968993257 " " absolute error = 1.00364161426114150000000000000E-13 " " relative error = 6.47348841198436000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35600000000000026 " " y[1] (analytic) = 1.5527950310559013 " " y[1] (numeric) = 1.5527950310560026 " " absolute error = 1.01252339845814280000000000000E-13 " " relative error = 6.5206506860704360000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35700000000000026 " " y[1] (analytic) = 1.5552099533437018 " " y[1] (numeric) = 1.5552099533438044 " " absolute error = 1.02584607475364460000000000000E-13 " " relative error = 6.596190260665933000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35800000000000026 " " y[1] (analytic) = 1.557632398753895 " " y[1] (numeric) = 1.5576323987539982 " " absolute error = 1.03250741290139560000000000000E-13 " " relative error = 6.6286975908269560000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.35900000000000026 " " y[1] (analytic) = 1.5600624024961003 " " y[1] (numeric) = 1.560062402496205 " " absolute error = 1.04805053524614780000000000000E-13 " " relative error = 6.718003930927805000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36000000000000026 " " y[1] (analytic) = 1.5625000000000009 " " y[1] (numeric) = 1.5625000000001064 " " absolute error = 1.05471187339389870000000000000E-13 " " relative error = 6.750155989720947000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36100000000000027 " " y[1] (analytic) = 1.5649452269170585 " " y[1] (numeric) = 1.5649452269171653 " " absolute error = 1.06803454968940060000000000000E-13 " " relative error = 6.824740772515266000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36200000000000027 " " y[1] (analytic) = 1.567398119122258 " " y[1] (numeric) = 1.5673981191223656 " " absolute error = 1.07691633388640180000000000000E-13 " " relative error = 6.87072621019523900000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36300000000000027 " " y[1] (analytic) = 1.569858712715856 " " y[1] (numeric) = 1.5698587127159653 " " absolute error = 1.0924594562311540000000000000E-13 " " relative error = 6.958966736192449000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36400000000000027 " " y[1] (analytic) = 1.572327044025158 " " y[1] (numeric) = 1.5723270440252681 " " absolute error = 1.10134124042815530000000000000E-13 " " relative error = 7.004530289123064000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36500000000000027 " " y[1] (analytic) = 1.5748031496062997 " " y[1] (numeric) = 1.5748031496064112 " " absolute error = 1.11466391672365720000000000000E-13 " " relative error = 7.0781158711952210000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36600000000000027 " " y[1] (analytic) = 1.5772870662460576 " " y[1] (numeric) = 1.57728706624617 " " absolute error = 1.12354570092065840000000000000E-13 " " relative error = 7.123279743836971000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36700000000000027 " " y[1] (analytic) = 1.5797788309636656 " " y[1] (numeric) = 1.5797788309637795 " " absolute error = 1.13908882326541060000000000000E-13 " " relative error = 7.210432251270047000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36800000000000027 " " y[1] (analytic) = 1.582278481012659 " " y[1] (numeric) = 1.582278481012774 " " absolute error = 1.14797060746241190000000000000E-13 " " relative error = 7.255174239162440000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.36900000000000027 " " y[1] (analytic) = 1.5847860538827263 " " y[1] (numeric) = 1.5847860538828427 " " absolute error = 1.1635137298071640000000000000E-13 " " relative error = 7.341771635083202000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3700000000000003 " " y[1] (analytic) = 1.587301587301588 " " y[1] (numeric) = 1.5873015873017053 " " absolute error = 1.17239551400416530000000000000E-13 " " relative error = 7.386091738226238000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3710000000000003 " " y[1] (analytic) = 1.5898251192368844 " " y[1] (numeric) = 1.5898251192370032 " " absolute error = 1.18793863634891750000000000000E-13 " " relative error = 7.47213402263469000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3720000000000003 " " y[1] (analytic) = 1.59235668789809 " " y[1] (numeric) = 1.5923566878982098 " " absolute error = 1.1990408665951690000000000000E-13 " " relative error = 7.529976642217658000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3730000000000003 " " y[1] (analytic) = 1.5948963317384375 " " y[1] (numeric) = 1.594896331738559 " " absolute error = 1.21458398893992130000000000000E-13 " " relative error = 7.615441610653303000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3740000000000003 " " y[1] (analytic) = 1.5974440894568698 " " y[1] (numeric) = 1.5974440894569923 " " absolute error = 1.22568621918617280000000000000E-13 " " relative error = 7.672795732105439000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3750000000000003 " " y[1] (analytic) = 1.6000000000000005 " " y[1] (numeric) = 1.6000000000001247 " " absolute error = 1.2412293415309250000000000000E-13 " " relative error = 7.757683384568278000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3760000000000003 " " y[1] (analytic) = 1.6025641025641035 " " y[1] (numeric) = 1.6025641025642285 " " absolute error = 1.25011112572792630000000000000E-13 " " relative error = 7.800693424542254000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3770000000000003 " " y[1] (analytic) = 1.6051364365971112 " " y[1] (numeric) = 1.605136436597238 " " absolute error = 1.26787469412192880000000000000E-13 " " relative error = 7.898859344379614000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3780000000000003 " " y[1] (analytic) = 1.607717041800644 " " y[1] (numeric) = 1.6077170418007718 " " absolute error = 1.27897692436818030000000000000E-13 " " relative error = 7.955236469570078000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3790000000000003 " " y[1] (analytic) = 1.6103059581320456 " " y[1] (numeric) = 1.6103059581321753 " " absolute error = 1.29674049276218280000000000000E-13 " " relative error = 8.052758460053152000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3800000000000003 " " y[1] (analytic) = 1.6129032258064524 " " y[1] (numeric) = 1.6129032258065832 " " absolute error = 1.30784272300843440000000000000E-13 " " relative error = 8.10862488265229000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3810000000000003 " " y[1] (analytic) = 1.6155088852988697 " " y[1] (numeric) = 1.615508885299002 " " absolute error = 1.32338584535318660000000000000E-13 " " relative error = 8.191758382736222000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3820000000000003 " " y[1] (analytic) = 1.6181229773462793 " " y[1] (numeric) = 1.6181229773464127 " " absolute error = 1.33448807559943820000000000000E-13 " " relative error = 8.247136307204522000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3830000000000003 " " y[1] (analytic) = 1.6207455429497575 " " y[1] (numeric) = 1.6207455429498927 " " absolute error = 1.35225164399344070000000000000E-13 " " relative error = 8.343392643439527000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3840000000000003 " " y[1] (analytic) = 1.6233766233766243 " " y[1] (numeric) = 1.6233766233767606 " " absolute error = 1.36335387423969220000000000000E-13 " " relative error = 8.3982598653165000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3850000000000003 " " y[1] (analytic) = 1.6260162601626023 " " y[1] (numeric) = 1.6260162601627404 " " absolute error = 1.38111744263369470000000000000E-13 " " relative error = 8.49387227219721900000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3860000000000003 " " y[1] (analytic) = 1.6286644951140075 " " y[1] (numeric) = 1.6286644951141467 " " absolute error = 1.39221967287994630000000000000E-13 " " relative error = 8.548228791482866000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3870000000000003 " " y[1] (analytic) = 1.6313213703099516 " " y[1] (numeric) = 1.6313213703100928 " " absolute error = 1.4122036873231990000000000000E-13 " " relative error = 8.656808603291208000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3880000000000003 " " y[1] (analytic) = 1.633986928104576 " " y[1] (numeric) = 1.6339869281047186 " " absolute error = 1.4255263636187010000000000000E-13 " " relative error = 8.724221345346446000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3890000000000003 " " y[1] (analytic) = 1.6366612111292969 " " y[1] (numeric) = 1.6366612111294412 " " absolute error = 1.44328993201270350000000000000E-13 " " relative error = 8.818501484597615000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3900000000000003 " " y[1] (analytic) = 1.639344262295083 " " y[1] (numeric) = 1.6393442622952286 " " absolute error = 1.45661260830820540000000000000E-13 " " relative error = 8.885336910680048000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3910000000000003 " " y[1] (analytic) = 1.6420361247947461 " " y[1] (numeric) = 1.6420361247948938 " " absolute error = 1.47659662275145820000000000000E-13 " " relative error = 8.992473432556377000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3920000000000003 " " y[1] (analytic) = 1.6447368421052642 " " y[1] (numeric) = 1.6447368421054132 " " absolute error = 1.489919299046960000000000000E-13 " " relative error = 9.058709338205512000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3930000000000003 " " y[1] (analytic) = 1.647446457990116 " " y[1] (numeric) = 1.647446457990267 " " absolute error = 1.5099033134902130000000000000E-13 " " relative error = 9.165113112885589000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3940000000000003 " " y[1] (analytic) = 1.650165016501651 " " y[1] (numeric) = 1.6501650165018036 " " absolute error = 1.5254464358349650000000000000E-13 " " relative error = 9.244205401159883000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3950000000000003 " " y[1] (analytic) = 1.6528925619834718 " " y[1] (numeric) = 1.652892561983626 " " absolute error = 1.54321000422896760000000000000E-13 " " relative error = 9.33642052558525100000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3960000000000003 " " y[1] (analytic) = 1.6556291390728486 " " y[1] (numeric) = 1.6556291390730045 " " absolute error = 1.55875312657371980000000000000E-13 " " relative error = 9.414868884505262000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3970000000000003 " " y[1] (analytic) = 1.6583747927031516 " " y[1] (numeric) = 1.6583747927033095 " " absolute error = 1.57873714101697260000000000000E-13 " " relative error = 9.51978496033234000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3980000000000003 " " y[1] (analytic) = 1.6611295681063132 " " y[1] (numeric) = 1.6611295681064726 " " absolute error = 1.59428026336172480000000000000E-13 " " relative error = 9.597567185437578000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.3990000000000003 " " y[1] (analytic) = 1.6638935108153086 " " y[1] (numeric) = 1.66389351081547 " " absolute error = 1.61426427780497760000000000000E-13 " " relative error = 9.701728309607912000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4000000000000003 " " y[1] (analytic) = 1.6666666666666676 " " y[1] (numeric) = 1.6666666666668306 " " absolute error = 1.62980740014972980000000000000E-13 " " relative error = 9.778844400898374000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4010000000000003 " " y[1] (analytic) = 1.6694490818030057 " " y[1] (numeric) = 1.6694490818031709 " " absolute error = 1.6520118606422330000000000000E-13 " " relative error = 9.895551045246972000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4020000000000003 " " y[1] (analytic) = 1.6722408026755864 " " y[1] (numeric) = 1.672240802675753 " " absolute error = 1.66533453693773480000000000000E-13 " " relative error = 9.958700530887648000000000000E-12 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4030000000000003 " " y[1] (analytic) = 1.675041876046902 " " y[1] (numeric) = 1.6750418760470707 " " absolute error = 1.6875389974302380000000000000E-13 " " relative error = 1.007460781465851600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4040000000000003 " " y[1] (analytic) = 1.6778523489932895 " " y[1] (numeric) = 1.67785234899346 " " absolute error = 1.70530256582424040000000000000E-13 " " relative error = 1.016360329231246800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4050000000000003 " " y[1] (analytic) = 1.6806722689075637 " " y[1] (numeric) = 1.6806722689077365 " " absolute error = 1.72750702631674360000000000000E-13 " " relative error = 1.02786668065846200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4060000000000003 " " y[1] (analytic) = 1.6835016835016845 " " y[1] (numeric) = 1.6835016835018588 " " absolute error = 1.74305014866149580000000000000E-13 " " relative error = 1.035371788304927800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4070000000000003 " " y[1] (analytic) = 1.6863406408094441 " " y[1] (numeric) = 1.6863406408096209 " " absolute error = 1.76747505520324920000000000000E-13 " " relative error = 1.048112707735526400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4080000000000003 " " y[1] (analytic) = 1.6891891891891901 " " y[1] (numeric) = 1.6891891891893684 " " absolute error = 1.78301817754800140000000000000E-13 " " relative error = 1.055546761108416300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4090000000000003 " " y[1] (analytic) = 1.6920473773265658 " " y[1] (numeric) = 1.6920473773267464 " " absolute error = 1.80522263804050450000000000000E-13 " " relative error = 1.066886579081937800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4100000000000003 " " y[1] (analytic) = 1.6949152542372892 " " y[1] (numeric) = 1.6949152542374715 " " absolute error = 1.8229862064345070000000000000E-13 " " relative error = 1.075561861796358500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4110000000000003 " " y[1] (analytic) = 1.6977928692699498 " " y[1] (numeric) = 1.6977928692701345 " " absolute error = 1.84741111297626050000000000000E-13 " " relative error = 1.08812514554301700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4120000000000003 " " y[1] (analytic) = 1.7006802721088445 " " y[1] (numeric) = 1.700680272109031 " " absolute error = 1.8651746813702630000000000000E-13 " " relative error = 1.09672271264571400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4130000000000003 " " y[1] (analytic) = 1.703577512776832 " " y[1] (numeric) = 1.703577512777021 " " absolute error = 1.88959958791201640000000000000E-13 " " relative error = 1.109194958104353200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4140000000000003 " " y[1] (analytic) = 1.7064846416382262 " " y[1] (numeric) = 1.706484641638417 " " absolute error = 1.9073631563060190000000000000E-13 " " relative error = 1.117714809595326500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4150000000000003 " " y[1] (analytic) = 1.7094017094017102 " " y[1] (numeric) = 1.7094017094019032 " " absolute error = 1.9295676167985220000000000000E-13 " " relative error = 1.12879705582713500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4160000000000003 " " y[1] (analytic) = 1.7123287671232887 " " y[1] (numeric) = 1.7123287671234837 " " absolute error = 1.9495516312417750000000000000E-13 " " relative error = 1.138538152645195800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4170000000000003 " " y[1] (analytic) = 1.7152658662092632 " " y[1] (numeric) = 1.7152658662094606 " " absolute error = 1.97397653778352830000000000000E-13 " " relative error = 1.150828321527796500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4180000000000003 " " y[1] (analytic) = 1.718213058419245 " " y[1] (numeric) = 1.7182130584194444 " " absolute error = 1.99396055222678110000000000000E-13 " " relative error = 1.16048504139598600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4190000000000003 " " y[1] (analytic) = 1.7211703958691917 " " y[1] (numeric) = 1.7211703958693938 " " absolute error = 2.0206059048177850000000000000E-13 " " relative error = 1.173972030699132600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4200000000000003 " " y[1] (analytic) = 1.7241379310344838 " " y[1] (numeric) = 1.7241379310346878 " " absolute error = 2.04058991926103770000000000000E-13 " " relative error = 1.183542153171401200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4210000000000003 " " y[1] (analytic) = 1.7271157167530233 " " y[1] (numeric) = 1.7271157167532298 " " absolute error = 2.06501482580279120000000000000E-13 " " relative error = 1.195643584139815400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4220000000000003 " " y[1] (analytic) = 1.730103806228375 " " y[1] (numeric) = 1.7301038062285834 " " absolute error = 2.0849988402460440000000000000E-13 " " relative error = 1.205129329662212600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4230000000000003 " " y[1] (analytic) = 1.7331022530329296 " " y[1] (numeric) = 1.733102253033141 " " absolute error = 2.1138646388862980000000000000E-13 " " relative error = 1.219699896637393600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4240000000000003 " " y[1] (analytic) = 1.7361111111111123 " " y[1] (numeric) = 1.7361111111113257 " " absolute error = 2.13384865332955100000000000000E-13 " " relative error = 1.229096824317820500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4250000000000003 " " y[1] (analytic) = 1.7391304347826095 " " y[1] (numeric) = 1.7391304347828258 " " absolute error = 2.1627144519698050000000000000E-13 " " relative error = 1.243560809882637100000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4260000000000003 " " y[1] (analytic) = 1.7421602787456458 " " y[1] (numeric) = 1.742160278745864 " " absolute error = 2.18269846641305780000000000000E-13 " " relative error = 1.252868919721094300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4270000000000003 " " y[1] (analytic) = 1.74520069808028 " " y[1] (numeric) = 1.7452006980805013 " " absolute error = 2.21156426505331180000000000000E-13 " " relative error = 1.267226323875547000000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4280000000000003 " " y[1] (analytic) = 1.7482517482517494 " " y[1] (numeric) = 1.7482517482519728 " " absolute error = 2.2337687255458150000000000000E-13 " " relative error = 1.277715711012205300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4290000000000003 " " y[1] (analytic) = 1.751313485113836 " " y[1] (numeric) = 1.7513134851140626 " " absolute error = 2.26485497023531930000000000000E-13 " " relative error = 1.293232188004366800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4300000000000003 " " y[1] (analytic) = 1.754385964912282 " " y[1] (numeric) = 1.7543859649125106 " " absolute error = 2.28705943072782250000000000000E-13 " " relative error = 1.303623875514858000000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4310000000000003 " " y[1] (analytic) = 1.7574692442882258 " " y[1] (numeric) = 1.7574692442884576 " " absolute error = 2.3181456754173269000000000000E-13 " " relative error = 1.319024889312458400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43200000000000033 " " y[1] (analytic) = 1.7605633802816913 " " y[1] (numeric) = 1.7605633802819256 " " absolute error = 2.34257058195908030000000000000E-13 " " relative error = 1.330580090552756800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43300000000000033 " " y[1] (analytic) = 1.7636684303350978 " " y[1] (numeric) = 1.7636684303353352 " " absolute error = 2.37365682664858470000000000000E-13 " " relative error = 1.34586342070974700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43400000000000033 " " y[1] (analytic) = 1.766784452296821 " " y[1] (numeric) = 1.7667844522970606 " " absolute error = 2.3958612871410878000000000000E-13 " " relative error = 1.356057488521855000000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43500000000000033 " " y[1] (analytic) = 1.7699115044247795 " " y[1] (numeric) = 1.7699115044250224 " " absolute error = 2.42916797787984250000000000000E-13 " " relative error = 1.372479907502110500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43600000000000033 " " y[1] (analytic) = 1.7730496453900721 " " y[1] (numeric) = 1.7730496453903175 " " absolute error = 2.4535928844215960000000000000E-13 " " relative error = 1.383826386813779300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43700000000000033 " " y[1] (analytic) = 1.7761989342806404 " " y[1] (numeric) = 1.7761989342808888 " " absolute error = 2.48467912911110030000000000000E-13 " " relative error = 1.398874349689548700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43800000000000033 " " y[1] (analytic) = 1.7793594306049834 " " y[1] (numeric) = 1.7793594306052345 " " absolute error = 2.5113244817021040000000000000E-13 " " relative error = 1.411364358716581700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.43900000000000033 " " y[1] (analytic) = 1.782531194295901 " " y[1] (numeric) = 1.7825311942961555 " " absolute error = 2.5446311724408590000000000000E-13 " " relative error = 1.427538087739321200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44000000000000034 " " y[1] (analytic) = 1.785714285714287 " " y[1] (numeric) = 1.785714285714544 " " absolute error = 2.57127652503186250000000000000E-13 " " relative error = 1.43991485401784220000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44100000000000034 " " y[1] (analytic) = 1.7889087656529525 " " y[1] (numeric) = 1.788908765653213 " " absolute error = 2.6045832157706170000000000000E-13 " " relative error = 1.455962017615774200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44200000000000034 " " y[1] (analytic) = 1.7921146953405032 " " y[1] (numeric) = 1.792114695340766 " " absolute error = 2.62900812231237070000000000000E-13 " " relative error = 1.466986532250301700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44300000000000034 " " y[1] (analytic) = 1.7953321364452433 " " y[1] (numeric) = 1.7953321364455097 " " absolute error = 2.66453525910037570000000000000E-13 " " relative error = 1.484146139318908700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44400000000000034 " " y[1] (analytic) = 1.798561151079138 " " y[1] (numeric) = 1.798561151079407 " " absolute error = 2.69118061169137950000000000000E-13 " " relative error = 1.49629642010040600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44500000000000034 " " y[1] (analytic) = 1.8018018018018027 " " y[1] (numeric) = 1.8018018018020754 " " absolute error = 2.72670774847938450000000000000E-13 " " relative error = 1.513322800406057800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44600000000000034 " " y[1] (analytic) = 1.80505415162455 " " y[1] (numeric) = 1.8050541516248255 " " absolute error = 2.75557354711963850000000000000E-13 " " relative error = 1.526587745104278600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44700000000000034 " " y[1] (analytic) = 1.8083182640144675 " " y[1] (numeric) = 1.8083182640147466 " " absolute error = 2.79110068390764350000000000000E-13 " " relative error = 1.54347867820092600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44800000000000034 " " y[1] (analytic) = 1.811594202898552 " " y[1] (numeric) = 1.811594202898834 " " absolute error = 2.81996648254789760000000000000E-13 " " relative error = 1.556621498366438000000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.44900000000000034 " " y[1] (analytic) = 1.8148820326678776 " " y[1] (numeric) = 1.8148820326681634 " " absolute error = 2.8577140653851530000000000000E-13 " " relative error = 1.574600450027218200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45000000000000034 " " y[1] (analytic) = 1.8181818181818195 " " y[1] (numeric) = 1.8181818181821083 " " absolute error = 2.88880031007465730000000000000E-13 " " relative error = 1.588840170541060700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45100000000000035 " " y[1] (analytic) = 1.8214936247723144 " " y[1] (numeric) = 1.8214936247726068 " " absolute error = 2.92432744686266230000000000000E-13 " " relative error = 1.605455768327600800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45200000000000035 " " y[1] (analytic) = 1.8248175182481765 " " y[1] (numeric) = 1.824817518248472 " " absolute error = 2.95541369155216670000000000000E-13 " " relative error = 1.61956670297058620000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45300000000000035 " " y[1] (analytic) = 1.8281535648994525 " " y[1] (numeric) = 1.828153564899752 " " absolute error = 2.99538172043867230000000000000E-13 " " relative error = 1.63847380107995300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45400000000000035 " " y[1] (analytic) = 1.8315018315018328 " " y[1] (numeric) = 1.8315018315021354 " " absolute error = 3.02646796512817700000000000000E-13 " " relative error = 1.652451508959983400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45500000000000035 " " y[1] (analytic) = 1.834862385321102 " " y[1] (numeric) = 1.8348623853214083 " " absolute error = 3.0642155479654320000000000000E-13 " " relative error = 1.669997473641159600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45600000000000035 " " y[1] (analytic) = 1.8382352941176485 " " y[1] (numeric) = 1.838235294117958 " " absolute error = 3.09530179265493640000000000000E-13 " " relative error = 1.68384417520428400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45700000000000035 " " y[1] (analytic) = 1.841620626151014 " " y[1] (numeric) = 1.8416206261513275 " " absolute error = 3.1352698215414420000000000000E-13 " " relative error = 1.702451513097002200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45800000000000035 " " y[1] (analytic) = 1.8450184501845033 " " y[1] (numeric) = 1.8450184501848201 " " absolute error = 3.1685765122801970000000000000E-13 " " relative error = 1.717368469655865300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.45900000000000035 " " y[1] (analytic) = 1.8484288354898346 " " y[1] (numeric) = 1.8484288354901557 " " absolute error = 3.21076498721595270000000000000E-13 " " relative error = 1.737023858083829600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46000000000000035 " " y[1] (analytic) = 1.8518518518518532 " " y[1] (numeric) = 1.8518518518521778 " " absolute error = 3.24629212400395800000000000000E-13 " " relative error = 1.75299774696213580000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46100000000000035 " " y[1] (analytic) = 1.855287569573285 " " y[1] (numeric) = 1.8552875695736137 " " absolute error = 3.28848059893971370000000000000E-13 " " relative error = 1.772491042828504600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46200000000000035 " " y[1] (analytic) = 1.8587360594795552 " " y[1] (numeric) = 1.8587360594798876 " " absolute error = 3.32400773572771870000000000000E-13 " " relative error = 1.788316161821511500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46300000000000036 " " y[1] (analytic) = 1.862197392923651 " " y[1] (numeric) = 1.8621973929239877 " " absolute error = 3.36619621066347460000000000000E-13 " " relative error = 1.80764736512628480000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46400000000000036 " " y[1] (analytic) = 1.8656716417910462 " " y[1] (numeric) = 1.8656716417913866 " " absolute error = 3.403943793500730000000000000E-13 " " relative error = 1.8245138733163901000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46500000000000036 " " y[1] (analytic) = 1.869158878504674 " " y[1] (numeric) = 1.8691588785050188 " " absolute error = 3.4483527144857360000000000000E-13 " " relative error = 1.84486870224986800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46600000000000036 " " y[1] (analytic) = 1.872659176029964 " " y[1] (numeric) = 1.8726591760303126 " " absolute error = 3.48610029732299150000000000000E-13 " " relative error = 1.86157755877047580000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46700000000000036 " " y[1] (analytic) = 1.8761726078799261 " " y[1] (numeric) = 1.8761726078802792 " " absolute error = 3.5305092183079980000000000000E-13 " " relative error = 1.881761413358161400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46800000000000036 " " y[1] (analytic) = 1.8796992481203023 " " y[1] (numeric) = 1.8796992481206594 " " absolute error = 3.57047724719450340000000000000E-13 " " relative error = 1.899493895507474400000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.46900000000000036 " " y[1] (analytic) = 1.8832391713747656 " " y[1] (numeric) = 1.8832391713751275 " " absolute error = 3.61932706027801030000000000000E-13 " " relative error = 1.921862669007622600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47000000000000036 " " y[1] (analytic) = 1.8867924528301903 " " y[1] (numeric) = 1.886792452830556 " " absolute error = 3.65707464311526560000000000000E-13 " " relative error = 1.938249560851089000000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47100000000000036 " " y[1] (analytic) = 1.8903591682419671 " " y[1] (numeric) = 1.8903591682423377 " " absolute error = 3.70592445619877250000000000000E-13 " " relative error = 1.960434037329149600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47200000000000036 " " y[1] (analytic) = 1.8939393939393954 " " y[1] (numeric) = 1.8939393939397702 " " absolute error = 3.74811293113452850000000000000E-13 " " relative error = 1.979003627639029600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47300000000000036 " " y[1] (analytic) = 1.8975332068311206 " " y[1] (numeric) = 1.8975332068315003 " " absolute error = 3.79696274421803540000000000000E-13 " " relative error = 2.000999366202903800000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47400000000000037 " " y[1] (analytic) = 1.901140684410648 " " y[1] (numeric) = 1.9011406844110317 " " absolute error = 3.8369307731045410000000000000E-13 " " relative error = 2.01822558665298700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47500000000000037 " " y[1] (analytic) = 1.904761904761906 " " y[1] (numeric) = 1.9047619047622948 " " absolute error = 3.8880010322372980000000000000E-13 " " relative error = 2.041200541924580200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47600000000000037 " " y[1] (analytic) = 1.9083969465648871 " " y[1] (numeric) = 1.9083969465652804 " " absolute error = 3.93240995322230450000000000000E-13 " " relative error = 2.060582815488486000000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47700000000000037 " " y[1] (analytic) = 1.9120458891013397 " " y[1] (numeric) = 1.9120458891017382 " " absolute error = 3.9857006584043120000000000000E-13 " " relative error = 2.08452144434545380000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47800000000000037 " " y[1] (analytic) = 1.915708812260538 " " y[1] (numeric) = 1.9157088122609411 " " absolute error = 4.03233002543856860000000000000E-13 " " relative error = 2.10487627327893100000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.47900000000000037 " " y[1] (analytic) = 1.9193857965451067 " " y[1] (numeric) = 1.9193857965455152 " " absolute error = 4.0856207306205760000000000000E-13 " " relative error = 2.12860840065331870000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.48000000000000037 " " y[1] (analytic) = 1.9230769230769247 " " y[1] (numeric) = 1.9230769230773377 " " absolute error = 4.13002965160558230000000000000E-13 " " relative error = 2.14761541883490100000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.48100000000000037 " " y[1] (analytic) = 1.926782273603084 " " y[1] (numeric) = 1.9267822736035025 " " absolute error = 4.185540802836840000000000000E-13 " " relative error = 2.172295676672318700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4820000000000004 " " y[1] (analytic) = 1.9305019305019322 " " y[1] (numeric) = 1.9305019305023554 " " absolute error = 4.23217016987109700000000000000E-13 " " relative error = 2.192264147993226200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4830000000000004 " " y[1] (analytic) = 1.9342359767891695 " " y[1] (numeric) = 1.9342359767895985 " " absolute error = 4.2899017671516050000000000000E-13 " " relative error = 2.217879213617378300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4840000000000004 " " y[1] (analytic) = 1.9379844961240327 " " y[1] (numeric) = 1.9379844961244665 " " absolute error = 4.3387515802351120000000000000E-13 " " relative error = 2.238795815401315700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4850000000000004 " " y[1] (analytic) = 1.9417475728155351 " " y[1] (numeric) = 1.941747572815975 " " absolute error = 4.398703623564870000000000000E-13 " " relative error = 2.26533236613590680000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4860000000000004 " " y[1] (analytic) = 1.9455252918287953 " " y[1] (numeric) = 1.9455252918292403 " " absolute error = 4.44977388269762740000000000000E-13 " " relative error = 2.287183775706578500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4870000000000004 " " y[1] (analytic) = 1.9493177387914242 " " y[1] (numeric) = 1.949317738791875 " " absolute error = 4.50750547997813560000000000000E-13 " " relative error = 2.312350311228782200000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4880000000000004 " " y[1] (analytic) = 1.9531250000000016 " " y[1] (numeric) = 1.9531250000004576 " " absolute error = 4.5607961851601430000000000000E-13 " " relative error = 2.335127646801991300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4890000000000004 " " y[1] (analytic) = 1.9569471624266157 " " y[1] (numeric) = 1.9569471624270778 " " absolute error = 4.6207482284899015000000000000E-13 " " relative error = 2.361202344758338300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4900000000000004 " " y[1] (analytic) = 1.960784313725492 " " y[1] (numeric) = 1.9607843137259593 " " absolute error = 4.6740389336719090000000000000E-13 " " relative error = 2.383759856172671700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4910000000000004 " " y[1] (analytic) = 1.964636542239687 " " y[1] (numeric) = 1.9646365422401606 " " absolute error = 4.7362114230509180000000000000E-13 " " relative error = 2.410731614332915500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4920000000000004 " " y[1] (analytic) = 1.9685039370078756 " " y[1] (numeric) = 1.9685039370083548 " " absolute error = 4.7917225742821756000000000000E-13 " " relative error = 2.43419506773534330000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4930000000000004 " " y[1] (analytic) = 1.9723865877712043 " " y[1] (numeric) = 1.97238658777169 " " absolute error = 4.8561155097104347000000000000E-13 " " relative error = 2.46205056342318900000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4940000000000004 " " y[1] (analytic) = 1.9762845849802388 " " y[1] (numeric) = 1.97628458498073 " " absolute error = 4.9116266609416925000000000000E-13 " " relative error = 2.485283090436494500000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4950000000000004 " " y[1] (analytic) = 1.9801980198019815 " " y[1] (numeric) = 1.9801980198024793 " " absolute error = 4.9782400424192020000000000000E-13 " " relative error = 2.514011221421695300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4960000000000004 " " y[1] (analytic) = 1.984126984126986 " " y[1] (numeric) = 1.9841269841274896 " " absolute error = 5.035971639699710000000000000E-13 " " relative error = 2.538129706408651700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4970000000000004 " " y[1] (analytic) = 1.9880715705765422 " " y[1] (numeric) = 1.9880715705770526 " " absolute error = 5.104805467226470000000000000E-13 " " relative error = 2.567717150014912600000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4980000000000004 " " y[1] (analytic) = 1.992031872509962 " " y[1] (numeric) = 1.9920318725104784 " " absolute error = 5.1647575105562280000000000000E-13 " " relative error = 2.592708270299224300000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.4990000000000004 " " y[1] (analytic) = 1.9960079840319374 " " y[1] (numeric) = 1.996007984032461 " " absolute error = 5.2358117841322380000000000000E-13 " " relative error = 2.623141703850249700000000000E-11 "%" h = 1.000E-3 " " " " "TOP MAIN SOLVE Loop" "NO POLE" x[1] = 0.5000000000000003 " " y[1] (analytic) = 2.0000000000000013 " " y[1] (numeric) = 2.0000000000005316 " " absolute error = 5.3024251656097480000000000000E-13 " " relative error = 2.651212582804872000000000000E-11 "%" h = 1.000E-3 " " "Finished!" "diff ( y , x , 1 ) = y * y;" Iterations = 500 "Total Elapsed Time "= 3 Minutes 2 Seconds "Elapsed Time(since restart) "= 3 Minutes 2 Seconds "Time to Timeout "= 11 Minutes 57 Seconds Percent Done = 100.20000000000007 "%" (%o51) true (%o51) diffeq.max