##############ECHO OF PROBLEM################# ##############temp/mtest6_rev_sm_hpostode.ode################# diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1)+x1; diff (x1,t,1) = 4.0 * x2 - 2.0 * diff (x2,t ,1) - 2.0 * x1; ! /* BEGIN FIRST INPUT BLOCK */ Digits=64; max_terms=30; ! /* END FIRST INPUT BLOCK */ /* BEGIN SECOND INPUT BLOCK */ /* # problem from Boyce DePrima - */ /* # _Elementary Differential Equations and Boundary Value Problems_ */ /* # page 269 */ /* # */ t_start=1.5; /* # did poorly with t_start := 0.5; */ t_end=5.0; array_x1_init[0 + 1] = exact_soln_x1(t_start); /* # I think following line should be omitted */ /* # diff(x1,1,exact_soln_x1p(t_start)); */ array_x2_init[0 + 1] = exact_soln_x2(t_start); array_x2_init[1 + 1] = exact_soln_x2p(t_start); glob_max_h=0.0000001; glob_look_poles=true; glob_max_iter=100; /* END SECOND INPUT BLOCK */ /* BEGIN OVERRIDE BLOCK */ glob_desired_digits_correct=10; glob_display_interval=0.001; glob_look_poles=true; glob_max_iter=10000000; glob_max_minutes=3; glob_subiter_method=3; /* END OVERRIDE BLOCK */ ! /* BEGIN USER DEF BLOCK */ double exact_soln_x1 (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(2.0 * c1 + 6.0 * c3 * exp(-t)); } double exact_soln_x1p (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( - 6.0 * c3 * exp(-t)); } double exact_soln_x2 (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return(c1 + c2 * exp(2.0 * t) + c3 * exp(-t)); } double exact_soln_x2p (double t) { double c1,c2,c3; c1 = 1.0; c2 = 0.0002; c3 = 0.0003; return( 2.0 * c2 * exp(2.0 * t) - c3 * exp(-t)); } /* END USER DEF BLOCK */ #######END OF ECHO OF PROBLEM################# START of Optimize min_size = 0 min_size = 1 opt_iter = 1 glob_desired_digits_correct = 10 desired_abs_gbl_error = 1e-10 range = 3.5 estimated_steps = 35000000 step_error = 2.857142857142857e-18 est_needed_step_err = 2.857142857142857e-18 hn_div_ho = 0.5 hn_div_ho_2 = 0.25 hn_div_ho_3 = 0.125 value3 = 7.130136767447456e-204 max_value3 = 7.130136767447456e-204 value3 = 7.130136767447456e-204 best_h = 1e-07 START of Soultion TOP MAIN SOLVE Loop t[1] = 1.5 x2[1] (analytic) = 1.004084046432682 x2[1] (numeric) = 1.004084046432682 absolute error = 0 relative error = 0 % Correct digits = 16 h = 1e-07 x1[1] (analytic) = 2.000401634288267 x1[1] (numeric) = 2.000401634288267 absolute error = 0 relative error = 0 % Correct digits = 16 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.277e+04 Order of pole = 5.966e+07 TOP MAIN SOLVE Loop t[1] = 1.501000000000584 x2[1] (analytic) = 1.00409202178144 x2[1] (numeric) = 1.004092021982522 absolute error = 2.01081595818664e-10 relative error = 2.002621188662659e-08 % Correct digits = 9 h = 1e-07 x1[1] (analytic) = 2.000401232854729 x1[1] (numeric) = 2.0004012324531 absolute error = 4.016285082286686e-10 relative error = 2.007739755566503e-08 % Correct digits = 9 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.279e+04 Order of pole = 5.972e+07 TOP MAIN SOLVE Loop t[1] = 1.502000000001168 x2[1] (analytic) = 1.004100013297674 x2[1] (numeric) = 1.004100014103079 absolute error = 8.054052980099868e-10 relative error = 8.02116609245794e-08 % Correct digits = 9 h = 1e-07 x1[1] (analytic) = 2.000400831822423 x1[1] (numeric) = 2.000400830215887 absolute error = 1.606536681464377e-09 relative error = 8.031073852337759e-08 % Correct digits = 9 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.28e+04 Order of pole = 5.978e+07 TOP MAIN SOLVE Loop t[1] = 1.503000000001752 x2[1] (analytic) = 1.00410802101355 x2[1] (numeric) = 1.004108022828134 absolute error = 1.814583150405724e-09 relative error = 1.8071593020182e-07 % Correct digits = 8 h = 1e-07 x1[1] (analytic) = 2.00040043119095 x1[1] (numeric) = 2.000400427576242 absolute error = 3.61470808840636e-09 relative error = 1.806992256172592e-07 % Correct digits = 8 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.282e+04 Order of pole = 5.984e+07 TOP MAIN SOLVE Loop t[1] = 1.504000000002335 x2[1] (analytic) = 1.004116044961301 x2[1] (numeric) = 1.004116048191537 absolute error = 3.230236078621829e-09 relative error = 3.216994783452866e-07 % Correct digits = 8 h = 1e-07 x1[1] (analytic) = 2.000400030959907 x1[1] (numeric) = 2.000400024533758 absolute error = 6.426149834481976e-09 relative error = 3.212432381036476e-07 % Correct digits = 8 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.284e+04 Order of pole = 5.99e+07 TOP MAIN SOLVE Loop t[1] = 1.505000000002919 x2[1] (analytic) = 1.004124085173221 x2[1] (numeric) = 1.004124090227195 absolute error = 5.053973239910192e-09 relative error = 5.033215829135631e-07 % Correct digits = 8 h = 1e-07 x1[1] (analytic) = 2.000399631128897 x1[1] (numeric) = 2.000399621088023 absolute error = 1.004087302192147e-08 relative error = 5.019433550012729e-07 % Correct digits = 8 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.285e+04 Order of pole = 5.996e+07 TOP MAIN SOLVE Loop t[1] = 1.506000000003503 x2[1] (analytic) = 1.004132141681672 x2[1] (numeric) = 1.004132148969098 absolute error = 7.287425773938594e-09 relative error = 7.257437015942907e-07 % Correct digits = 8 h = 1e-07 x1[1] (analytic) = 2.000399231697516 x1[1] (numeric) = 2.000399217238638 absolute error = 1.445887809481405e-08 relative error = 7.227996224805793e-07 % Correct digits = 8 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.287e+04 Order of pole = 6.002e+07 TOP MAIN SOLVE Loop t[1] = 1.507000000004087 x2[1] (analytic) = 1.004140214519079 x2[1] (numeric) = 1.004140224451297 absolute error = 9.932217492902851e-09 relative error = 9.891265531736287e-07 % Correct digits = 8 h = 1e-07 x1[1] (analytic) = 2.000398832665368 x1[1] (numeric) = 2.000398812985209 absolute error = 1.968015928e-08 relative error = 9.838117758636058e-07 % Correct digits = 8 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.289e+04 Order of pole = 6.008e+07 TOP MAIN SOLVE Loop t[1] = 1.508000000004671 x2[1] (analytic) = 1.004148303717934 x2[1] (numeric) = 1.00414831670791 absolute error = 1.298997598375706e-08 relative error = 1.293631223163024e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000398434032052 x1[1] (numeric) = 2.000398408327325 absolute error = 2.570472767970955e-08 relative error = 1.284980394025728e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.29e+04 Order of pole = 6.014e+07 TOP MAIN SOLVE Loop t[1] = 1.509000000005255 x2[1] (analytic) = 1.004156409310791 x2[1] (numeric) = 1.004156425773129 absolute error = 1.646233815932874e-08 relative error = 1.639419716558674e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.00039803579717 x1[1] (numeric) = 2.000398003264577 absolute error = 3.253259306390532e-08 relative error = 1.62630598919484e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.292e+04 Order of pole = 6.02e+07 TOP MAIN SOLVE Loop t[1] = 1.510000000005839 x2[1] (analytic) = 1.004164531330274 x2[1] (numeric) = 1.004164551681214 absolute error = 2.035094026631157e-08 relative error = 2.026653962707836e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000397637960325 x1[1] (numeric) = 2.000397597796562 absolute error = 4.016376253801468e-08 relative error = 2.007788940351232e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.294e+04 Order of pole = 6.026e+07 TOP MAIN SOLVE Loop t[1] = 1.511000000006423 x2[1] (analytic) = 1.004172669809068 x2[1] (numeric) = 1.004172694466493 absolute error = 2.465742499069279e-08 relative error = 2.45549652286206e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000397240521117 x1[1] (numeric) = 2.000397191922881 absolute error = 4.859823565794841e-08 relative error = 2.42942924902697e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.295e+04 Order of pole = 6.032e+07 TOP MAIN SOLVE Loop t[1] = 1.512000000007006 x2[1] (analytic) = 1.004180824779927 x2[1] (numeric) = 1.004180854163366 absolute error = 2.938343857117331e-08 relative error = 2.926110302655191e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000396843479149 x1[1] (numeric) = 2.000396785643119 absolute error = 5.78360301872749e-08 relative error = 2.891227826908823e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.297e+04 Order of pole = 6.038e+07 TOP MAIN SOLVE Loop t[1] = 1.51300000000759 x2[1] (analytic) = 1.004188996275669 x2[1] (numeric) = 1.004189030806293 absolute error = 3.453062413782959e-08 relative error = 3.438657888693921e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000396446834026 x1[1] (numeric) = 2.00039637895688 absolute error = 6.787714523781574e-08 relative error = 3.393184653234267e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.298e+04 Order of pole = 6.044e+07 TOP MAIN SOLVE Loop t[1] = 1.514000000008174 x2[1] (analytic) = 1.004197184329178 x2[1] (numeric) = 1.004197224429819 absolute error = 4.010064125203883e-08 relative error = 3.993303494355722e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000396050585349 x1[1] (numeric) = 2.000395971863745 absolute error = 7.872160345812063e-08 relative error = 3.935300883796756e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 Complex estimate of poles used for equation 2 Radius of convergence = 3221 Order of pole = 1.422e+07 TOP MAIN SOLVE Loop t[1] = 1.515000000008758 x2[1] (analytic) = 1.004205388973404 x2[1] (numeric) = 1.004205435068545 absolute error = 4.609514125952785e-08 relative error = 4.590210505308158e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000395654732722 x1[1] (numeric) = 2.000395564363321 absolute error = 9.036940173956509e-08 relative error = 4.51757638673933e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.302e+04 Order of pole = 6.056e+07 TOP MAIN SOLVE Loop t[1] = 1.516000000009342 x2[1] (analytic) = 1.004213610241365 x2[1] (numeric) = 1.004213662757148 absolute error = 5.251578372167387e-08 relative error = 5.229543115737259e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000395259275751 x1[1] (numeric) = 2.000395156455189 absolute error = 1.028205613984312e-07 relative error = 5.140012251161688e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.303e+04 Order of pole = 6.063e+07 TOP MAIN SOLVE Loop t[1] = 1.517000000009926 x2[1] (analytic) = 1.004221848166142 x2[1] (numeric) = 1.004221907530377 absolute error = 5.936423508323685e-08 relative error = 5.911466195606551e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000394864214039 x1[1] (numeric) = 2.000394748138945 absolute error = 1.160750944251276e-07 relative error = 5.802609099915575e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.305e+04 Order of pole = 6.069e+07 TOP MAIN SOLVE Loop t[1] = 1.51800000001051 x2[1] (analytic) = 1.004230102780886 x2[1] (numeric) = 1.004230169423047 absolute error = 6.664216156693215e-08 relative error = 6.636144583038145e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000394469547191 x1[1] (numeric) = 2.000394339414178 absolute error = 1.301330128100631e-07 relative error = 6.505367555806129e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.307e+04 Order of pole = 6.075e+07 TOP MAIN SOLVE Loop t[1] = 1.519000000011093 x2[1] (analytic) = 1.004238374118811 x2[1] (numeric) = 1.004238448470042 absolute error = 7.435123072774275e-08 relative error = 7.403743239047572e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000394075274813 x1[1] (numeric) = 2.000393930280482 absolute error = 1.44994330320003e-07 relative error = 7.248288330392288e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.308e+04 Order of pole = 6.081e+07 TOP MAIN SOLVE Loop t[1] = 1.520000000011677 x2[1] (analytic) = 1.004246662213201 x2[1] (numeric) = 1.004246744706326 absolute error = 8.249312499764017e-08 relative error = 8.214428596241322e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000393681396509 x1[1] (numeric) = 2.000393520737448 absolute error = 1.606590611658021e-07 relative error = 8.031372157386701e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.31e+04 Order of pole = 6.087e+07 TOP MAIN SOLVE Loop t[1] = 1.521000000012261 x2[1] (analytic) = 1.004254967097405 x2[1] (numeric) = 1.004255058166921 absolute error = 9.106951592841028e-08 relative error = 9.068365993909718e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000393287911888 x1[1] (numeric) = 2.00039311078467 absolute error = 1.771272177819583e-07 relative error = 8.854619681655336e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.312e+04 Order of pole = 6.093e+07 TOP MAIN SOLVE Loop t[1] = 1.522000000012845 x2[1] (analytic) = 1.004263288804838 x2[1] (numeric) = 1.004263388886935 absolute error = 1.000820961660764e-07 relative error = 9.965722861898391e-06 % Correct digits = 7 h = 1e-07 x1[1] (analytic) = 2.000392894820554 x1[1] (numeric) = 2.00039270042173 absolute error = 1.943988241492889e-07 relative error = 9.718032125220455e-06 % Correct digits = 7 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.313e+04 Order of pole = 6.099e+07 TOP MAIN SOLVE Loop t[1] = 1.523000000013429 x2[1] (analytic) = 1.004271627368986 x2[1] (numeric) = 1.00427173690152 absolute error = 1.095325348199339e-07 relative error = 1.090666427636612e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000392502122115 x1[1] (numeric) = 2.000392289648221 absolute error = 2.124738944786486e-07 relative error = 1.062161022165629e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.315e+04 Order of pole = 6.105e+07 TOP MAIN SOLVE Loop t[1] = 1.524000000014013 x2[1] (analytic) = 1.004279982823397 x2[1] (numeric) = 1.004280102245927 absolute error = 1.194225305312102e-07 relative error = 1.189135824408946e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000392109816179 x1[1] (numeric) = 2.000391878463738 absolute error = 2.313524407604461e-07 relative error = 1.156535459349046e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.317e+04 Order of pole = 6.111e+07 TOP MAIN SOLVE Loop t[1] = 1.525000000014597 x2[1] (analytic) = 1.00428835520169 x2[1] (numeric) = 1.004288484955469 absolute error = 1.297537792765979e-07 relative error = 1.291997249639917e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000391717902352 x1[1] (numeric) = 2.000391466867867 absolute error = 2.51034485199142e-07 relative error = 1.25492663738071e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.318e+04 Order of pole = 6.117e+07 TOP MAIN SOLVE Loop t[1] = 1.526000000015181 x2[1] (analytic) = 1.004296744537551 x2[1] (numeric) = 1.004296885065528 absolute error = 1.405279774768786e-07 relative error = 1.39926747986809e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000391326380243 x1[1] (numeric) = 2.000391054860185 absolute error = 2.715200579928023e-07 relative error = 1.357334709524684e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.32e+04 Order of pole = 6.123e+07 TOP MAIN SOLVE Loop t[1] = 1.527000000015764 x2[1] (analytic) = 1.004305150864732 x2[1] (numeric) = 1.004305302611556 absolute error = 1.517468244394138e-07 relative error = 1.510963319353246e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.00039093524946 x1[1] (numeric) = 2.000390642440292 absolute error = 2.928091684673007e-07 relative error = 1.46375972469994e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.322e+04 Order of pole = 6.13e+07 TOP MAIN SOLVE Loop t[1] = 1.528000000016348 x2[1] (analytic) = 1.004313574217055 x2[1] (numeric) = 1.004313737629083 absolute error = 1.634120276872153e-07 relative error = 1.627101653132672e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000390544509613 x1[1] (numeric) = 2.000390229607775 absolute error = 3.149018383830082e-07 relative error = 1.57420179398121e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.323e+04 Order of pole = 6.136e+07 TOP MAIN SOLVE Loop t[1] = 1.529000000016932 x2[1] (analytic) = 1.004322014628408 x2[1] (numeric) = 1.004322190153705 absolute error = 1.755252967416965e-07 relative error = 1.747699385108466e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000390154160311 x1[1] (numeric) = 2.000389816362218 absolute error = 3.3779809305301e-07 relative error = 1.688661046198786e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.325e+04 Order of pole = 6.142e+07 TOP MAIN SOLVE Loop t[1] = 1.530000000017516 x2[1] (analytic) = 1.004330472132748 x2[1] (numeric) = 1.00433066022109 absolute error = 1.880883420124491e-07 relative error = 1.872773426987968e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000389764201162 x1[1] (numeric) = 2.000389402703209 absolute error = 3.614979529054096e-07 relative error = 1.807137585758297e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.327e+04 Order of pole = 6.148e+07 TOP MAIN SOLVE Loop t[1] = 1.5310000000181 x2[1] (analytic) = 1.004338946764101 x2[1] (numeric) = 1.004339147866982 absolute error = 2.011028814585814e-07 relative error = 2.002340764604605e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000389374631778 x1[1] (numeric) = 2.000388988630336 absolute error = 3.860014419210245e-07 relative error = 1.929631534820954e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.328e+04 Order of pole = 6.154e+07 TOP MAIN SOLVE Loop t[1] = 1.532000000018684 x2[1] (analytic) = 1.004347438556559 x2[1] (numeric) = 1.004347653127193 absolute error = 2.145706345935139e-07 relative error = 2.136418398217785e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000388985451768 x1[1] (numeric) = 2.000388574143177 absolute error = 4.113085916301884e-07 relative error = 2.056143053283701e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.33e+04 Order of pole = 6.16e+07 TOP MAIN SOLVE Loop t[1] = 1.533000000019268 x2[1] (analytic) = 1.004355947544284 x2[1] (numeric) = 1.004356176037611 absolute error = 2.284933269258715e-07 relative error = 2.275023386724126e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000388596660744 x1[1] (numeric) = 2.000388159241327 absolute error = 4.374194171319346e-07 relative error = 2.18667221889847e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 Complex estimate of poles used for equation 2 Radius of convergence = 1.119e+04 Order of pole = 2.939e+06 TOP MAIN SOLVE Loop t[1] = 1.534000000019851 x2[1] (analytic) = 1.004364473761507 x2[1] (numeric) = 1.004364716634195 absolute error = 2.428726884051713e-07 relative error = 2.418172832174897e-05 % Correct digits = 6 h = 1e-07 x1[1] (analytic) = 2.000388208258317 x1[1] (numeric) = 2.000387743924367 absolute error = 4.643339495125076e-07 relative error = 2.321219189333207e-05 % Correct digits = 6 h = 1e-07 Complex estimate of poles used for equation 1 NO POLE for equation 2 Radius of convergence = 3.333e+04 Order of pole = 6.173e+07 Finished! Maximum Time Reached before Solution Completed! diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1)+x1; diff (x1,t,1) = 4.0 * x2 - 2.0 * diff (x2,t ,1) - 2.0 * x1; Iterations = 342126 Total Elapsed Time = 3 Minutes 0 Seconds Elapsed Time(since restart) = 3 Minutes 0 Seconds Expected Time Remaining = 5 Hours 3 Minutes 54 Seconds Optimized Time Remaining = 5 Hours 3 Minutes 54 Seconds Expected Total Time = 5 Hours 6 Minutes 54 Seconds Time to Timeout = 0 Seconds Percent Done = 0.9775 %