Demonstration Problem for DLSOIBT Galerkin method solution of system of 3 PDEs: u(i) = -(u(1)+u(2)+u(3)) u(i) + eta(i) u(i) (i=1,2,3) t x xx x interval is -1 to 1, zero boundary conditions x discretized using piecewise linear basis functions Fixed parameters are as follows: Diffusion coefficients are eta = 0.10E+00 0.20E-01 0.10E-01 t0 = 0.00000E+00 tlast = 0.40000E+00 Uniform mesh, number of intervals = 100 Block size mb = 3 Number of blocks nb = 99 ODE system size neq = 297 Initial profiles: Values of PDE component i = 1 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.1000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.2000E+00 0.1000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 Values of PDE component i = 2 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.1500E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.1500E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 Values of PDE component i = 3 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.2500E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.2500E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 ****************************************************************************************** Run with rtol = 0.1E-02 atol = 0.1E-02 mf = 11 At time t = 0.10000E+00 current h = 0.22309E-01 current order = 2 current nst = 17 At time t = 0.20000E+00 current h = 0.41111E-01 current order = 2 current nst = 20 At time t = 0.30000E+00 current h = 0.51073E-01 current order = 2 current nst = 22 At time t = 0.40000E+00 current h = 0.58735E-01 current order = 2 current nst = 24 Values of PDE component i = 1 0.1714E-02 0.3444E-02 0.5203E-02 0.7006E-02 0.8866E-02 0.1079E-01 0.1280E-01 0.1490E-01 0.1710E-01 0.1940E-01 0.2181E-01 0.2433E-01 0.2695E-01 0.2968E-01 0.3252E-01 0.3546E-01 0.3848E-01 0.4161E-01 0.4476E-01 0.4807E-01 0.5130E-01 0.5475E-01 0.5807E-01 0.6152E-01 0.6497E-01 0.6846E-01 0.7191E-01 0.7531E-01 0.7899E-01 0.8240E-01 0.8592E-01 0.8937E-01 0.9290E-01 0.9639E-01 0.9989E-01 0.1034E+00 0.1069E+00 0.1104E+00 0.1139E+00 0.1174E+00 0.1209E+00 0.1244E+00 0.1279E+00 0.1313E+00 0.1348E+00 0.1383E+00 0.1417E+00 0.1451E+00 0.1484E+00 0.1517E+00 0.1549E+00 0.1580E+00 0.1610E+00 0.1639E+00 0.1666E+00 0.1692E+00 0.1717E+00 0.1739E+00 0.1760E+00 0.1778E+00 0.1794E+00 0.1808E+00 0.1819E+00 0.1828E+00 0.1834E+00 0.1836E+00 0.1836E+00 0.1832E+00 0.1825E+00 0.1814E+00 0.1799E+00 0.1779E+00 0.1755E+00 0.1725E+00 0.1691E+00 0.1650E+00 0.1603E+00 0.1548E+00 0.1488E+00 0.1424E+00 0.1348E+00 0.1266E+00 0.1176E+00 0.1082E+00 0.9850E-01 0.8871E-01 0.7919E-01 0.7011E-01 0.6161E-01 0.5373E-01 0.4646E-01 0.3977E-01 0.3361E-01 0.2792E-01 0.2263E-01 0.1770E-01 0.1303E-01 0.8572E-02 0.4251E-02 Values of PDE component i = 2 0.7744E-05 0.1814E-04 0.3450E-04 0.6154E-04 0.1063E-03 0.1791E-03 0.2948E-03 0.4744E-03 0.7460E-03 0.1146E-02 0.1721E-02 0.2524E-02 0.3617E-02 0.5068E-02 0.6946E-02 0.9315E-02 0.1223E-01 0.1575E-01 0.1989E-01 0.2467E-01 0.3009E-01 0.3613E-01 0.4275E-01 0.4991E-01 0.5755E-01 0.6565E-01 0.7409E-01 0.8288E-01 0.9193E-01 0.1012E+00 0.1106E+00 0.1202E+00 0.1299E+00 0.1396E+00 0.1494E+00 0.1592E+00 0.1690E+00 0.1787E+00 0.1885E+00 0.1981E+00 0.2076E+00 0.2170E+00 0.2260E+00 0.2348E+00 0.2431E+00 0.2510E+00 0.2582E+00 0.2649E+00 0.2709E+00 0.2762E+00 0.2808E+00 0.2847E+00 0.2880E+00 0.2907E+00 0.2929E+00 0.2946E+00 0.2960E+00 0.2970E+00 0.2978E+00 0.2984E+00 0.2989E+00 0.2992E+00 0.2994E+00 0.2996E+00 0.2997E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2998E+00 0.2997E+00 0.2995E+00 0.2990E+00 0.2983E+00 0.2967E+00 0.2937E+00 0.2890E+00 0.2817E+00 0.2700E+00 0.2505E+00 0.2224E+00 0.1867E+00 0.1467E+00 0.1076E+00 0.7410E-01 0.4856E-01 0.3063E-01 0.1876E-01 0.1121E-01 0.6560E-02 0.3762E-02 0.2117E-02 0.1169E-02 0.6331E-03 0.3338E-03 0.1667E-03 0.6863E-04 Values of PDE component i = 3 0.9946E-07 0.2111E-06 0.3530E-06 0.5634E-06 0.9413E-06 0.1748E-05 0.3651E-05 0.8255E-05 0.1918E-04 0.4406E-04 0.9777E-04 0.2074E-03 0.4186E-03 0.8031E-03 0.1465E-02 0.2544E-02 0.4215E-02 0.6676E-02 0.1014E-01 0.1480E-01 0.2085E-01 0.2842E-01 0.3759E-01 0.4839E-01 0.6081E-01 0.7478E-01 0.9019E-01 0.1069E+00 0.1248E+00 0.1437E+00 0.1634E+00 0.1839E+00 0.2049E+00 0.2264E+00 0.2481E+00 0.2701E+00 0.2921E+00 0.3140E+00 0.3356E+00 0.3565E+00 0.3766E+00 0.3955E+00 0.4129E+00 0.4286E+00 0.4425E+00 0.4545E+00 0.4646E+00 0.4730E+00 0.4798E+00 0.4851E+00 0.4892E+00 0.4923E+00 0.4946E+00 0.4962E+00 0.4974E+00 0.4983E+00 0.4988E+00 0.4992E+00 0.4995E+00 0.4997E+00 0.4998E+00 0.4999E+00 0.4999E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5001E+00 0.5002E+00 0.4993E+00 0.4986E+00 0.4997E+00 0.4870E+00 0.4577E+00 0.4177E+00 0.3451E+00 0.2443E+00 0.1485E+00 0.7868E-01 0.3717E-01 0.1608E-01 0.6517E-02 0.2512E-02 0.9293E-03 0.3318E-03 0.1145E-03 0.3809E-04 0.1208E-04 0.3554E-05 0.9159E-06 0.1828E-06 Final statistics for mf = 11: 24 steps, 38 res, 9 jacobians, rwork size = 7447, iwork size = 317 Final output is correct to within 0.75E+01 times local error tolerance. ****************************************************************************************** Run with rtol = 0.1E-02 atol = 0.1E-02 mf = 12 At time t = 0.10000E+00 current h = 0.22305E-01 current order = 2 current nst = 17 At time t = 0.20000E+00 current h = 0.41173E-01 current order = 2 current nst = 20 At time t = 0.30000E+00 current h = 0.51349E-01 current order = 2 current nst = 22 At time t = 0.40000E+00 current h = 0.59258E-01 current order = 2 current nst = 24 Values of PDE component i = 1 0.1712E-02 0.3440E-02 0.5197E-02 0.6997E-02 0.8855E-02 0.1078E-01 0.1279E-01 0.1488E-01 0.1708E-01 0.1937E-01 0.2178E-01 0.2429E-01 0.2691E-01 0.2963E-01 0.3247E-01 0.3540E-01 0.3841E-01 0.4153E-01 0.4468E-01 0.4797E-01 0.5119E-01 0.5464E-01 0.5793E-01 0.6141E-01 0.6483E-01 0.6828E-01 0.7182E-01 0.7520E-01 0.7881E-01 0.8223E-01 0.8579E-01 0.8925E-01 0.9278E-01 0.9628E-01 0.9978E-01 0.1033E+00 0.1068E+00 0.1103E+00 0.1138E+00 0.1173E+00 0.1208E+00 0.1243E+00 0.1278E+00 0.1313E+00 0.1348E+00 0.1383E+00 0.1417E+00 0.1451E+00 0.1485E+00 0.1517E+00 0.1549E+00 0.1580E+00 0.1611E+00 0.1639E+00 0.1667E+00 0.1693E+00 0.1717E+00 0.1739E+00 0.1760E+00 0.1778E+00 0.1794E+00 0.1808E+00 0.1819E+00 0.1828E+00 0.1834E+00 0.1837E+00 0.1836E+00 0.1833E+00 0.1826E+00 0.1814E+00 0.1800E+00 0.1779E+00 0.1756E+00 0.1726E+00 0.1692E+00 0.1652E+00 0.1605E+00 0.1553E+00 0.1493E+00 0.1427E+00 0.1352E+00 0.1271E+00 0.1182E+00 0.1088E+00 0.9905E-01 0.8926E-01 0.7972E-01 0.7062E-01 0.6210E-01 0.5418E-01 0.4687E-01 0.4013E-01 0.3392E-01 0.2818E-01 0.2285E-01 0.1786E-01 0.1315E-01 0.8653E-02 0.4292E-02 Values of PDE component i = 2 0.7637E-05 0.1793E-04 0.3419E-04 0.6115E-04 0.1058E-03 0.1786E-03 0.2943E-03 0.4739E-03 0.7456E-03 0.1146E-02 0.1720E-02 0.2523E-02 0.3617E-02 0.5068E-02 0.6944E-02 0.9312E-02 0.1223E-01 0.1574E-01 0.1988E-01 0.2465E-01 0.3006E-01 0.3608E-01 0.4269E-01 0.4983E-01 0.5745E-01 0.6552E-01 0.7394E-01 0.8271E-01 0.9173E-01 0.1010E+00 0.1104E+00 0.1200E+00 0.1297E+00 0.1395E+00 0.1493E+00 0.1591E+00 0.1690E+00 0.1788E+00 0.1885E+00 0.1982E+00 0.2077E+00 0.2171E+00 0.2261E+00 0.2349E+00 0.2432E+00 0.2510E+00 0.2583E+00 0.2650E+00 0.2710E+00 0.2763E+00 0.2808E+00 0.2848E+00 0.2880E+00 0.2907E+00 0.2929E+00 0.2946E+00 0.2960E+00 0.2970E+00 0.2978E+00 0.2984E+00 0.2989E+00 0.2992E+00 0.2994E+00 0.2996E+00 0.2997E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2998E+00 0.2997E+00 0.2995E+00 0.2990E+00 0.2982E+00 0.2967E+00 0.2939E+00 0.2892E+00 0.2818E+00 0.2698E+00 0.2505E+00 0.2233E+00 0.1883E+00 0.1484E+00 0.1090E+00 0.7518E-01 0.4929E-01 0.3109E-01 0.1904E-01 0.1137E-01 0.6649E-02 0.3811E-02 0.2144E-02 0.1184E-02 0.6409E-03 0.3380E-03 0.1688E-03 0.6954E-04 Values of PDE component i = 3 0.4971E-08 0.1654E-07 0.4938E-07 0.1430E-06 0.4036E-06 0.1108E-05 0.2946E-05 0.7554E-05 0.1859E-04 0.4372E-04 0.9783E-04 0.2080E-03 0.4199E-03 0.8052E-03 0.1468E-02 0.2548E-02 0.4220E-02 0.6682E-02 0.1015E-01 0.1482E-01 0.2087E-01 0.2845E-01 0.3764E-01 0.4847E-01 0.6093E-01 0.7494E-01 0.9039E-01 0.1072E+00 0.1251E+00 0.1440E+00 0.1638E+00 0.1843E+00 0.2053E+00 0.2268E+00 0.2485E+00 0.2704E+00 0.2924E+00 0.3143E+00 0.3358E+00 0.3567E+00 0.3768E+00 0.3956E+00 0.4130E+00 0.4287E+00 0.4426E+00 0.4546E+00 0.4647E+00 0.4731E+00 0.4798E+00 0.4851E+00 0.4892E+00 0.4923E+00 0.4946E+00 0.4963E+00 0.4974E+00 0.4983E+00 0.4988E+00 0.4992E+00 0.4995E+00 0.4997E+00 0.4998E+00 0.4999E+00 0.4999E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5001E+00 0.5002E+00 0.4993E+00 0.4985E+00 0.4997E+00 0.4872E+00 0.4581E+00 0.4191E+00 0.3472E+00 0.2462E+00 0.1500E+00 0.7966E-01 0.3769E-01 0.1632E-01 0.6615E-02 0.2551E-02 0.9450E-03 0.3386E-03 0.1179E-03 0.4004E-04 0.1329E-04 0.4317E-05 0.1366E-05 0.3917E-06 Final statistics for mf = 12: 24 steps, 128 res, 9 jacobians, rwork size = 7447, iwork size = 317 Final output is correct to within 0.71E+01 times local error tolerance. ****************************************************************************************** Run with rtol = 0.1E-02 atol = 0.1E-02 mf = 21 At time t = 0.10000E+00 current h = 0.17227E-01 current order = 2 current nst = 20 At time t = 0.20000E+00 current h = 0.33122E-01 current order = 2 current nst = 24 At time t = 0.30000E+00 current h = 0.33122E-01 current order = 2 current nst = 27 At time t = 0.40000E+00 current h = 0.33122E-01 current order = 2 current nst = 30 Values of PDE component i = 1 0.1714E-02 0.3443E-02 0.5202E-02 0.7005E-02 0.8866E-02 0.1080E-01 0.1281E-01 0.1491E-01 0.1711E-01 0.1942E-01 0.2183E-01 0.2435E-01 0.2698E-01 0.2971E-01 0.3255E-01 0.3549E-01 0.3851E-01 0.4162E-01 0.4480E-01 0.4804E-01 0.5134E-01 0.5468E-01 0.5807E-01 0.6148E-01 0.6492E-01 0.6838E-01 0.7185E-01 0.7533E-01 0.7882E-01 0.8230E-01 0.8579E-01 0.8928E-01 0.9277E-01 0.9626E-01 0.9975E-01 0.1032E+00 0.1067E+00 0.1102E+00 0.1137E+00 0.1173E+00 0.1208E+00 0.1243E+00 0.1278E+00 0.1314E+00 0.1349E+00 0.1384E+00 0.1419E+00 0.1453E+00 0.1487E+00 0.1520E+00 0.1552E+00 0.1583E+00 0.1613E+00 0.1642E+00 0.1670E+00 0.1696E+00 0.1720E+00 0.1742E+00 0.1762E+00 0.1780E+00 0.1796E+00 0.1809E+00 0.1820E+00 0.1829E+00 0.1834E+00 0.1837E+00 0.1836E+00 0.1832E+00 0.1825E+00 0.1813E+00 0.1798E+00 0.1779E+00 0.1755E+00 0.1726E+00 0.1692E+00 0.1652E+00 0.1606E+00 0.1553E+00 0.1493E+00 0.1426E+00 0.1351E+00 0.1269E+00 0.1180E+00 0.1085E+00 0.9878E-01 0.8900E-01 0.7947E-01 0.7039E-01 0.6188E-01 0.5397E-01 0.4668E-01 0.3996E-01 0.3376E-01 0.2804E-01 0.2273E-01 0.1777E-01 0.1308E-01 0.8605E-02 0.4268E-02 Values of PDE component i = 2 0.8292E-05 0.1926E-04 0.3621E-04 0.6381E-04 0.1089E-03 0.1817E-03 0.2968E-03 0.4749E-03 0.7439E-03 0.1141E-02 0.1711E-02 0.2511E-02 0.3603E-02 0.5056E-02 0.6939E-02 0.9318E-02 0.1225E-01 0.1578E-01 0.1993E-01 0.2471E-01 0.3012E-01 0.3614E-01 0.4273E-01 0.4984E-01 0.5743E-01 0.6546E-01 0.7385E-01 0.8257E-01 0.9156E-01 0.1008E+00 0.1102E+00 0.1197E+00 0.1294E+00 0.1392E+00 0.1491E+00 0.1591E+00 0.1691E+00 0.1791E+00 0.1891E+00 0.1990E+00 0.2088E+00 0.2183E+00 0.2274E+00 0.2362E+00 0.2444E+00 0.2520E+00 0.2591E+00 0.2655E+00 0.2712E+00 0.2762E+00 0.2806E+00 0.2844E+00 0.2875E+00 0.2902E+00 0.2924E+00 0.2941E+00 0.2955E+00 0.2966E+00 0.2975E+00 0.2981E+00 0.2986E+00 0.2990E+00 0.2993E+00 0.2995E+00 0.2996E+00 0.2998E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2998E+00 0.2995E+00 0.2992E+00 0.2985E+00 0.2972E+00 0.2948E+00 0.2904E+00 0.2820E+00 0.2686E+00 0.2491E+00 0.2219E+00 0.1869E+00 0.1474E+00 0.1085E+00 0.7511E-01 0.4945E-01 0.3134E-01 0.1929E-01 0.1159E-01 0.6825E-02 0.3944E-02 0.2239E-02 0.1249E-02 0.6833E-03 0.3643E-03 0.1837E-03 0.7624E-04 Values of PDE component i = 3 0.6162E-07 0.1377E-06 0.2513E-06 0.4520E-06 0.8610E-06 0.1779E-05 0.3932E-05 0.8985E-05 0.2057E-04 0.4616E-04 0.1003E-03 0.2094E-03 0.4185E-03 0.7992E-03 0.1457E-02 0.2534E-02 0.4208E-02 0.6681E-02 0.1016E-01 0.1485E-01 0.2091E-01 0.2846E-01 0.3758E-01 0.4830E-01 0.6060E-01 0.7441E-01 0.8966E-01 0.1062E+00 0.1240E+00 0.1428E+00 0.1626E+00 0.1831E+00 0.2045E+00 0.2264E+00 0.2488E+00 0.2714E+00 0.2941E+00 0.3166E+00 0.3384E+00 0.3594E+00 0.3791E+00 0.3974E+00 0.4141E+00 0.4290E+00 0.4422E+00 0.4536E+00 0.4633E+00 0.4714E+00 0.4781E+00 0.4834E+00 0.4877E+00 0.4910E+00 0.4935E+00 0.4954E+00 0.4968E+00 0.4978E+00 0.4985E+00 0.4990E+00 0.4993E+00 0.4996E+00 0.4997E+00 0.4998E+00 0.4999E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5001E+00 0.5002E+00 0.5007E+00 0.5022E+00 0.5016E+00 0.4950E+00 0.4878E+00 0.4672E+00 0.4132E+00 0.3339E+00 0.2415E+00 0.1512E+00 0.8188E-01 0.3937E-01 0.1726E-01 0.7067E-02 0.2745E-02 0.1022E-02 0.3663E-03 0.1269E-03 0.4241E-04 0.1362E-04 0.4154E-05 0.1175E-05 0.2869E-06 Final statistics for mf = 21: 30 steps, 45 res, 11 jacobians, rwork size = 5368, iwork size = 317 Final output is correct to within 0.14E+02 times local error tolerance. ****************************************************************************************** Run with rtol = 0.1E-02 atol = 0.1E-02 mf = 22 At time t = 0.10000E+00 current h = 0.17219E-01 current order = 2 current nst = 20 At time t = 0.20000E+00 current h = 0.33216E-01 current order = 2 current nst = 24 At time t = 0.30000E+00 current h = 0.33216E-01 current order = 2 current nst = 27 At time t = 0.40000E+00 current h = 0.33216E-01 current order = 2 current nst = 30 Values of PDE component i = 1 0.1713E-02 0.3441E-02 0.5199E-02 0.7001E-02 0.8861E-02 0.1079E-01 0.1280E-01 0.1490E-01 0.1710E-01 0.1940E-01 0.2181E-01 0.2433E-01 0.2695E-01 0.2968E-01 0.3252E-01 0.3545E-01 0.3846E-01 0.4156E-01 0.4473E-01 0.4797E-01 0.5126E-01 0.5460E-01 0.5797E-01 0.6138E-01 0.6481E-01 0.6826E-01 0.7172E-01 0.7520E-01 0.7868E-01 0.8216E-01 0.8565E-01 0.8915E-01 0.9264E-01 0.9614E-01 0.9964E-01 0.1031E+00 0.1066E+00 0.1102E+00 0.1137E+00 0.1172E+00 0.1207E+00 0.1243E+00 0.1278E+00 0.1314E+00 0.1349E+00 0.1384E+00 0.1419E+00 0.1453E+00 0.1487E+00 0.1520E+00 0.1552E+00 0.1584E+00 0.1614E+00 0.1643E+00 0.1670E+00 0.1696E+00 0.1720E+00 0.1742E+00 0.1762E+00 0.1780E+00 0.1796E+00 0.1810E+00 0.1820E+00 0.1829E+00 0.1834E+00 0.1837E+00 0.1836E+00 0.1832E+00 0.1825E+00 0.1814E+00 0.1798E+00 0.1779E+00 0.1755E+00 0.1726E+00 0.1691E+00 0.1651E+00 0.1605E+00 0.1552E+00 0.1493E+00 0.1426E+00 0.1351E+00 0.1270E+00 0.1181E+00 0.1087E+00 0.9897E-01 0.8921E-01 0.7968E-01 0.7060E-01 0.6208E-01 0.5416E-01 0.4685E-01 0.4011E-01 0.3390E-01 0.2816E-01 0.2283E-01 0.1785E-01 0.1314E-01 0.8644E-02 0.4287E-02 Values of PDE component i = 2 0.8220E-05 0.1912E-04 0.3600E-04 0.6353E-04 0.1086E-03 0.1813E-03 0.2965E-03 0.4746E-03 0.7437E-03 0.1141E-02 0.1711E-02 0.2511E-02 0.3603E-02 0.5056E-02 0.6939E-02 0.9317E-02 0.1225E-01 0.1577E-01 0.1992E-01 0.2470E-01 0.3010E-01 0.3611E-01 0.4268E-01 0.4978E-01 0.5735E-01 0.6535E-01 0.7372E-01 0.8242E-01 0.9140E-01 0.1006E+00 0.1100E+00 0.1196E+00 0.1293E+00 0.1391E+00 0.1490E+00 0.1590E+00 0.1691E+00 0.1792E+00 0.1892E+00 0.1991E+00 0.2088E+00 0.2183E+00 0.2275E+00 0.2362E+00 0.2444E+00 0.2521E+00 0.2591E+00 0.2655E+00 0.2712E+00 0.2763E+00 0.2806E+00 0.2844E+00 0.2876E+00 0.2902E+00 0.2924E+00 0.2941E+00 0.2955E+00 0.2966E+00 0.2975E+00 0.2981E+00 0.2986E+00 0.2990E+00 0.2993E+00 0.2995E+00 0.2997E+00 0.2998E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2999E+00 0.2998E+00 0.2995E+00 0.2992E+00 0.2984E+00 0.2971E+00 0.2946E+00 0.2900E+00 0.2819E+00 0.2688E+00 0.2494E+00 0.2221E+00 0.1872E+00 0.1478E+00 0.1089E+00 0.7549E-01 0.4976E-01 0.3156E-01 0.1943E-01 0.1168E-01 0.6879E-02 0.3975E-02 0.2257E-02 0.1259E-02 0.6888E-03 0.3673E-03 0.1853E-03 0.7689E-04 Values of PDE component i = 3 0.6430E-08 0.2134E-07 0.6312E-07 0.1798E-06 0.4960E-06 0.1323E-05 0.3404E-05 0.8433E-05 0.2008E-04 0.4585E-04 0.1003E-03 0.2099E-03 0.4196E-03 0.8010E-03 0.1459E-02 0.2537E-02 0.4212E-02 0.6685E-02 0.1017E-01 0.1486E-01 0.2092E-01 0.2848E-01 0.3762E-01 0.4836E-01 0.6068E-01 0.7453E-01 0.8981E-01 0.1064E+00 0.1242E+00 0.1431E+00 0.1628E+00 0.1835E+00 0.2048E+00 0.2267E+00 0.2490E+00 0.2716E+00 0.2943E+00 0.3167E+00 0.3385E+00 0.3595E+00 0.3792E+00 0.3975E+00 0.4141E+00 0.4290E+00 0.4422E+00 0.4536E+00 0.4633E+00 0.4714E+00 0.4780E+00 0.4834E+00 0.4877E+00 0.4910E+00 0.4935E+00 0.4954E+00 0.4968E+00 0.4978E+00 0.4985E+00 0.4990E+00 0.4993E+00 0.4996E+00 0.4997E+00 0.4998E+00 0.4999E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5001E+00 0.5002E+00 0.5008E+00 0.5024E+00 0.5018E+00 0.4949E+00 0.4881E+00 0.4679E+00 0.4136E+00 0.3344E+00 0.2425E+00 0.1521E+00 0.8249E-01 0.3969E-01 0.1741E-01 0.7129E-02 0.2769E-02 0.1031E-02 0.3705E-03 0.1289E-03 0.4358E-04 0.1433E-04 0.4592E-05 0.1427E-05 0.4012E-06 Final statistics for mf = 22: 30 steps, 155 res, 11 jacobians, rwork size = 5368, iwork size = 317 Final output is correct to within 0.14E+02 times local error tolerance. ****************************************************************************************** Run with rtol = 0.1E-05 atol = 0.1E-05 mf = 11 At time t = 0.10000E+00 current h = 0.79980E-02 current order = 4 current nst = 60 At time t = 0.20000E+00 current h = 0.10116E-01 current order = 4 current nst = 71 At time t = 0.30000E+00 current h = 0.12090E-01 current order = 4 current nst = 89 At time t = 0.40000E+00 current h = 0.12812E-01 current order = 4 current nst = 106 Values of PDE component i = 1 0.1710E-02 0.3434E-02 0.5188E-02 0.6985E-02 0.8839E-02 0.1076E-01 0.1277E-01 0.1486E-01 0.1705E-01 0.1934E-01 0.2174E-01 0.2425E-01 0.2687E-01 0.2960E-01 0.3243E-01 0.3536E-01 0.3838E-01 0.4149E-01 0.4468E-01 0.4793E-01 0.5124E-01 0.5460E-01 0.5799E-01 0.6143E-01 0.6489E-01 0.6836E-01 0.7186E-01 0.7536E-01 0.7887E-01 0.8239E-01 0.8591E-01 0.8943E-01 0.9295E-01 0.9647E-01 0.9999E-01 0.1035E+00 0.1070E+00 0.1106E+00 0.1141E+00 0.1176E+00 0.1211E+00 0.1246E+00 0.1281E+00 0.1316E+00 0.1350E+00 0.1384E+00 0.1418E+00 0.1452E+00 0.1485E+00 0.1517E+00 0.1549E+00 0.1580E+00 0.1609E+00 0.1638E+00 0.1665E+00 0.1691E+00 0.1715E+00 0.1737E+00 0.1758E+00 0.1776E+00 0.1793E+00 0.1806E+00 0.1818E+00 0.1827E+00 0.1833E+00 0.1836E+00 0.1836E+00 0.1832E+00 0.1825E+00 0.1814E+00 0.1799E+00 0.1780E+00 0.1756E+00 0.1727E+00 0.1693E+00 0.1653E+00 0.1606E+00 0.1554E+00 0.1494E+00 0.1427E+00 0.1353E+00 0.1271E+00 0.1183E+00 0.1089E+00 0.9916E-01 0.8935E-01 0.7978E-01 0.7067E-01 0.6212E-01 0.5420E-01 0.4688E-01 0.4015E-01 0.3394E-01 0.2820E-01 0.2286E-01 0.1788E-01 0.1316E-01 0.8659E-02 0.4295E-02 Values of PDE component i = 2 0.7174E-05 0.1708E-04 0.3312E-04 0.6016E-04 0.1053E-03 0.1792E-03 0.2967E-03 0.4789E-03 0.7536E-03 0.1157E-02 0.1734E-02 0.2538E-02 0.3631E-02 0.5078E-02 0.6948E-02 0.9306E-02 0.1221E-01 0.1572E-01 0.1985E-01 0.2462E-01 0.3004E-01 0.3608E-01 0.4271E-01 0.4988E-01 0.5755E-01 0.6566E-01 0.7415E-01 0.8298E-01 0.9208E-01 0.1014E+00 0.1109E+00 0.1206E+00 0.1303E+00 0.1401E+00 0.1500E+00 0.1599E+00 0.1697E+00 0.1794E+00 0.1891E+00 0.1986E+00 0.2079E+00 0.2170E+00 0.2258E+00 0.2343E+00 0.2425E+00 0.2501E+00 0.2573E+00 0.2639E+00 0.2700E+00 0.2754E+00 0.2802E+00 0.2843E+00 0.2878E+00 0.2907E+00 0.2931E+00 0.2949E+00 0.2964E+00 0.2975E+00 0.2983E+00 0.2988E+00 0.2992E+00 0.2995E+00 0.2997E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.2999E+00 0.2999E+00 0.2997E+00 0.2994E+00 0.2989E+00 0.2980E+00 0.2964E+00 0.2937E+00 0.2890E+00 0.2814E+00 0.2693E+00 0.2509E+00 0.2245E+00 0.1899E+00 0.1499E+00 0.1099E+00 0.7540E-01 0.4903E-01 0.3061E-01 0.1852E-01 0.1091E-01 0.6277E-02 0.3530E-02 0.1941E-02 0.1042E-02 0.5461E-03 0.2774E-03 0.1334E-03 0.5328E-04 Values of PDE component i = 3 0.1970E-09 0.1996E-08 0.1198E-07 0.5563E-07 0.2186E-06 0.7563E-06 0.2357E-05 0.6709E-05 0.1762E-04 0.4304E-04 0.9826E-04 0.2107E-03 0.4262E-03 0.8157E-03 0.1482E-02 0.2562E-02 0.4229E-02 0.6681E-02 0.1013E-01 0.1479E-01 0.2084E-01 0.2843E-01 0.3766E-01 0.4855E-01 0.6109E-01 0.7522E-01 0.9082E-01 0.1078E+00 0.1259E+00 0.1450E+00 0.1650E+00 0.1857E+00 0.2069E+00 0.2284E+00 0.2500E+00 0.2717E+00 0.2932E+00 0.3143E+00 0.3350E+00 0.3550E+00 0.3741E+00 0.3922E+00 0.4092E+00 0.4248E+00 0.4389E+00 0.4514E+00 0.4623E+00 0.4715E+00 0.4791E+00 0.4852E+00 0.4899E+00 0.4933E+00 0.4958E+00 0.4974E+00 0.4985E+00 0.4992E+00 0.4996E+00 0.4998E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.4998E+00 0.4994E+00 0.4982E+00 0.4949E+00 0.4861E+00 0.4652E+00 0.4219E+00 0.3479E+00 0.2496E+00 0.1516E+00 0.7801E-01 0.3480E-01 0.1387E-01 0.5059E-02 0.1712E-02 0.5396E-03 0.1582E-03 0.4287E-04 0.1061E-04 0.2351E-05 0.4522E-06 0.7093E-07 0.7792E-08 Final statistics for mf = 11: 106 steps, 148 res, 23 jacobians, rwork size = 7447, iwork size = 317 Final output is correct to within 0.67E+01 times local error tolerance. ****************************************************************************************** Run with rtol = 0.1E-05 atol = 0.1E-05 mf = 12 At time t = 0.10000E+00 current h = 0.64185E-02 current order = 3 current nst = 61 At time t = 0.20000E+00 current h = 0.13037E-01 current order = 4 current nst = 71 At time t = 0.30000E+00 current h = 0.68710E-02 current order = 3 current nst = 89 At time t = 0.40000E+00 current h = 0.84070E-02 current order = 3 current nst = 100 Values of PDE component i = 1 0.1710E-02 0.3434E-02 0.5188E-02 0.6985E-02 0.8839E-02 0.1076E-01 0.1276E-01 0.1486E-01 0.1705E-01 0.1934E-01 0.2174E-01 0.2425E-01 0.2687E-01 0.2960E-01 0.3243E-01 0.3536E-01 0.3838E-01 0.4149E-01 0.4467E-01 0.4793E-01 0.5124E-01 0.5460E-01 0.5800E-01 0.6142E-01 0.6489E-01 0.6836E-01 0.7186E-01 0.7535E-01 0.7888E-01 0.8238E-01 0.8591E-01 0.8942E-01 0.9295E-01 0.9646E-01 0.1000E+00 0.1035E+00 0.1070E+00 0.1105E+00 0.1141E+00 0.1176E+00 0.1211E+00 0.1246E+00 0.1281E+00 0.1316E+00 0.1350E+00 0.1384E+00 0.1418E+00 0.1452E+00 0.1485E+00 0.1517E+00 0.1549E+00 0.1580E+00 0.1609E+00 0.1638E+00 0.1665E+00 0.1691E+00 0.1715E+00 0.1737E+00 0.1758E+00 0.1776E+00 0.1793E+00 0.1806E+00 0.1818E+00 0.1827E+00 0.1833E+00 0.1836E+00 0.1836E+00 0.1832E+00 0.1825E+00 0.1814E+00 0.1799E+00 0.1780E+00 0.1756E+00 0.1727E+00 0.1693E+00 0.1653E+00 0.1606E+00 0.1554E+00 0.1494E+00 0.1427E+00 0.1353E+00 0.1271E+00 0.1183E+00 0.1089E+00 0.9916E-01 0.8935E-01 0.7978E-01 0.7067E-01 0.6212E-01 0.5420E-01 0.4688E-01 0.4015E-01 0.3394E-01 0.2820E-01 0.2286E-01 0.1788E-01 0.1316E-01 0.8659E-02 0.4295E-02 Values of PDE component i = 2 0.7174E-05 0.1708E-04 0.3312E-04 0.6016E-04 0.1053E-03 0.1792E-03 0.2967E-03 0.4789E-03 0.7536E-03 0.1157E-02 0.1734E-02 0.2538E-02 0.3631E-02 0.5078E-02 0.6948E-02 0.9306E-02 0.1221E-01 0.1572E-01 0.1985E-01 0.2462E-01 0.3004E-01 0.3608E-01 0.4271E-01 0.4988E-01 0.5755E-01 0.6566E-01 0.7415E-01 0.8298E-01 0.9208E-01 0.1014E+00 0.1109E+00 0.1206E+00 0.1303E+00 0.1401E+00 0.1500E+00 0.1599E+00 0.1697E+00 0.1794E+00 0.1891E+00 0.1986E+00 0.2079E+00 0.2170E+00 0.2258E+00 0.2343E+00 0.2425E+00 0.2501E+00 0.2573E+00 0.2639E+00 0.2700E+00 0.2754E+00 0.2802E+00 0.2843E+00 0.2878E+00 0.2907E+00 0.2931E+00 0.2949E+00 0.2964E+00 0.2975E+00 0.2983E+00 0.2988E+00 0.2992E+00 0.2995E+00 0.2997E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.2999E+00 0.2999E+00 0.2997E+00 0.2994E+00 0.2989E+00 0.2980E+00 0.2964E+00 0.2937E+00 0.2890E+00 0.2814E+00 0.2693E+00 0.2509E+00 0.2245E+00 0.1899E+00 0.1499E+00 0.1099E+00 0.7540E-01 0.4902E-01 0.3061E-01 0.1852E-01 0.1091E-01 0.6278E-02 0.3530E-02 0.1941E-02 0.1042E-02 0.5461E-03 0.2775E-03 0.1334E-03 0.5329E-04 Values of PDE component i = 3 0.2003E-09 0.2003E-08 0.1198E-07 0.5562E-07 0.2185E-06 0.7561E-06 0.2356E-05 0.6708E-05 0.1762E-04 0.4303E-04 0.9826E-04 0.2107E-03 0.4262E-03 0.8157E-03 0.1482E-02 0.2562E-02 0.4229E-02 0.6681E-02 0.1013E-01 0.1479E-01 0.2084E-01 0.2843E-01 0.3766E-01 0.4855E-01 0.6109E-01 0.7522E-01 0.9082E-01 0.1078E+00 0.1259E+00 0.1450E+00 0.1650E+00 0.1857E+00 0.2069E+00 0.2284E+00 0.2500E+00 0.2717E+00 0.2932E+00 0.3143E+00 0.3350E+00 0.3550E+00 0.3741E+00 0.3922E+00 0.4092E+00 0.4248E+00 0.4389E+00 0.4514E+00 0.4623E+00 0.4715E+00 0.4791E+00 0.4852E+00 0.4899E+00 0.4933E+00 0.4958E+00 0.4974E+00 0.4985E+00 0.4992E+00 0.4996E+00 0.4998E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.4998E+00 0.4994E+00 0.4982E+00 0.4949E+00 0.4861E+00 0.4652E+00 0.4219E+00 0.3479E+00 0.2496E+00 0.1516E+00 0.7801E-01 0.3480E-01 0.1387E-01 0.5061E-02 0.1713E-02 0.5401E-03 0.1584E-03 0.4295E-04 0.1063E-04 0.2361E-05 0.4556E-06 0.7202E-07 0.8103E-08 Final statistics for mf = 12: 100 steps, 353 res, 22 jacobians, rwork size = 7447, iwork size = 317 Final output is correct to within 0.12E+02 times local error tolerance. ****************************************************************************************** Run with rtol = 0.1E-05 atol = 0.1E-05 mf = 21 At time t = 0.10000E+00 current h = 0.57520E-02 current order = 4 current nst = 72 At time t = 0.20000E+00 current h = 0.92622E-02 current order = 5 current nst = 85 At time t = 0.30000E+00 current h = 0.12274E-01 current order = 5 current nst = 95 At time t = 0.40000E+00 current h = 0.12274E-01 current order = 5 current nst = 103 Values of PDE component i = 1 0.1710E-02 0.3434E-02 0.5188E-02 0.6985E-02 0.8839E-02 0.1076E-01 0.1277E-01 0.1486E-01 0.1705E-01 0.1934E-01 0.2174E-01 0.2425E-01 0.2687E-01 0.2960E-01 0.3243E-01 0.3536E-01 0.3838E-01 0.4149E-01 0.4468E-01 0.4793E-01 0.5124E-01 0.5460E-01 0.5799E-01 0.6143E-01 0.6489E-01 0.6836E-01 0.7186E-01 0.7536E-01 0.7887E-01 0.8239E-01 0.8591E-01 0.8943E-01 0.9295E-01 0.9647E-01 0.9999E-01 0.1035E+00 0.1070E+00 0.1106E+00 0.1141E+00 0.1176E+00 0.1211E+00 0.1246E+00 0.1281E+00 0.1316E+00 0.1350E+00 0.1384E+00 0.1418E+00 0.1452E+00 0.1485E+00 0.1517E+00 0.1549E+00 0.1580E+00 0.1609E+00 0.1638E+00 0.1665E+00 0.1691E+00 0.1715E+00 0.1737E+00 0.1758E+00 0.1776E+00 0.1793E+00 0.1806E+00 0.1818E+00 0.1827E+00 0.1833E+00 0.1836E+00 0.1836E+00 0.1832E+00 0.1825E+00 0.1814E+00 0.1799E+00 0.1780E+00 0.1756E+00 0.1727E+00 0.1693E+00 0.1653E+00 0.1606E+00 0.1554E+00 0.1494E+00 0.1427E+00 0.1353E+00 0.1271E+00 0.1183E+00 0.1089E+00 0.9916E-01 0.8935E-01 0.7978E-01 0.7067E-01 0.6212E-01 0.5420E-01 0.4688E-01 0.4015E-01 0.3394E-01 0.2820E-01 0.2286E-01 0.1788E-01 0.1316E-01 0.8659E-02 0.4295E-02 Values of PDE component i = 2 0.7174E-05 0.1708E-04 0.3312E-04 0.6016E-04 0.1053E-03 0.1792E-03 0.2967E-03 0.4789E-03 0.7536E-03 0.1157E-02 0.1734E-02 0.2538E-02 0.3631E-02 0.5078E-02 0.6948E-02 0.9306E-02 0.1221E-01 0.1572E-01 0.1985E-01 0.2462E-01 0.3004E-01 0.3608E-01 0.4271E-01 0.4988E-01 0.5755E-01 0.6566E-01 0.7415E-01 0.8298E-01 0.9208E-01 0.1014E+00 0.1109E+00 0.1206E+00 0.1303E+00 0.1401E+00 0.1500E+00 0.1599E+00 0.1697E+00 0.1794E+00 0.1891E+00 0.1986E+00 0.2079E+00 0.2170E+00 0.2258E+00 0.2343E+00 0.2425E+00 0.2501E+00 0.2573E+00 0.2639E+00 0.2700E+00 0.2754E+00 0.2802E+00 0.2843E+00 0.2878E+00 0.2907E+00 0.2931E+00 0.2949E+00 0.2964E+00 0.2975E+00 0.2983E+00 0.2988E+00 0.2993E+00 0.2995E+00 0.2997E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.2999E+00 0.2999E+00 0.2997E+00 0.2994E+00 0.2989E+00 0.2980E+00 0.2964E+00 0.2937E+00 0.2890E+00 0.2814E+00 0.2693E+00 0.2509E+00 0.2245E+00 0.1899E+00 0.1499E+00 0.1099E+00 0.7540E-01 0.4903E-01 0.3061E-01 0.1852E-01 0.1091E-01 0.6277E-02 0.3530E-02 0.1941E-02 0.1042E-02 0.5460E-03 0.2774E-03 0.1333E-03 0.5327E-04 Values of PDE component i = 3 0.1721E-09 0.1938E-08 0.1186E-07 0.5543E-07 0.2183E-06 0.7559E-06 0.2356E-05 0.6708E-05 0.1762E-04 0.4304E-04 0.9826E-04 0.2107E-03 0.4262E-03 0.8157E-03 0.1482E-02 0.2562E-02 0.4228E-02 0.6681E-02 0.1013E-01 0.1479E-01 0.2084E-01 0.2843E-01 0.3766E-01 0.4855E-01 0.6109E-01 0.7522E-01 0.9082E-01 0.1078E+00 0.1259E+00 0.1450E+00 0.1650E+00 0.1857E+00 0.2069E+00 0.2284E+00 0.2500E+00 0.2717E+00 0.2932E+00 0.3143E+00 0.3350E+00 0.3550E+00 0.3741E+00 0.3922E+00 0.4092E+00 0.4248E+00 0.4389E+00 0.4514E+00 0.4623E+00 0.4715E+00 0.4791E+00 0.4852E+00 0.4899E+00 0.4933E+00 0.4958E+00 0.4974E+00 0.4985E+00 0.4992E+00 0.4996E+00 0.4998E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.4998E+00 0.4994E+00 0.4982E+00 0.4949E+00 0.4861E+00 0.4652E+00 0.4218E+00 0.3479E+00 0.2497E+00 0.1517E+00 0.7801E-01 0.3480E-01 0.1387E-01 0.5057E-02 0.1711E-02 0.5392E-03 0.1581E-03 0.4283E-04 0.1060E-04 0.2350E-05 0.4530E-06 0.7165E-07 0.8130E-08 Final statistics for mf = 21: 103 steps, 122 res, 17 jacobians, rwork size = 5368, iwork size = 317 Final output is correct to within 0.11E+02 times local error tolerance. ****************************************************************************************** Run with rtol = 0.1E-05 atol = 0.1E-05 mf = 22 At time t = 0.10000E+00 current h = 0.57538E-02 current order = 4 current nst = 72 At time t = 0.20000E+00 current h = 0.92575E-02 current order = 5 current nst = 85 At time t = 0.30000E+00 current h = 0.12294E-01 current order = 5 current nst = 95 At time t = 0.40000E+00 current h = 0.12294E-01 current order = 5 current nst = 103 Values of PDE component i = 1 0.1710E-02 0.3434E-02 0.5188E-02 0.6985E-02 0.8839E-02 0.1076E-01 0.1277E-01 0.1486E-01 0.1705E-01 0.1934E-01 0.2174E-01 0.2425E-01 0.2687E-01 0.2960E-01 0.3243E-01 0.3536E-01 0.3838E-01 0.4149E-01 0.4468E-01 0.4793E-01 0.5124E-01 0.5460E-01 0.5799E-01 0.6143E-01 0.6489E-01 0.6836E-01 0.7186E-01 0.7536E-01 0.7887E-01 0.8239E-01 0.8591E-01 0.8943E-01 0.9295E-01 0.9647E-01 0.9999E-01 0.1035E+00 0.1070E+00 0.1106E+00 0.1141E+00 0.1176E+00 0.1211E+00 0.1246E+00 0.1281E+00 0.1316E+00 0.1350E+00 0.1384E+00 0.1418E+00 0.1452E+00 0.1485E+00 0.1517E+00 0.1549E+00 0.1580E+00 0.1609E+00 0.1638E+00 0.1665E+00 0.1691E+00 0.1715E+00 0.1737E+00 0.1758E+00 0.1776E+00 0.1793E+00 0.1806E+00 0.1818E+00 0.1827E+00 0.1833E+00 0.1836E+00 0.1836E+00 0.1832E+00 0.1825E+00 0.1814E+00 0.1799E+00 0.1780E+00 0.1756E+00 0.1727E+00 0.1693E+00 0.1653E+00 0.1606E+00 0.1554E+00 0.1494E+00 0.1427E+00 0.1353E+00 0.1271E+00 0.1183E+00 0.1089E+00 0.9916E-01 0.8935E-01 0.7978E-01 0.7067E-01 0.6212E-01 0.5420E-01 0.4688E-01 0.4015E-01 0.3394E-01 0.2820E-01 0.2286E-01 0.1788E-01 0.1316E-01 0.8659E-02 0.4295E-02 Values of PDE component i = 2 0.7174E-05 0.1708E-04 0.3312E-04 0.6016E-04 0.1053E-03 0.1792E-03 0.2967E-03 0.4789E-03 0.7536E-03 0.1157E-02 0.1734E-02 0.2538E-02 0.3631E-02 0.5078E-02 0.6948E-02 0.9306E-02 0.1221E-01 0.1572E-01 0.1985E-01 0.2462E-01 0.3004E-01 0.3608E-01 0.4271E-01 0.4988E-01 0.5755E-01 0.6566E-01 0.7415E-01 0.8298E-01 0.9208E-01 0.1014E+00 0.1109E+00 0.1206E+00 0.1303E+00 0.1401E+00 0.1500E+00 0.1599E+00 0.1697E+00 0.1794E+00 0.1891E+00 0.1986E+00 0.2079E+00 0.2170E+00 0.2258E+00 0.2343E+00 0.2425E+00 0.2501E+00 0.2573E+00 0.2639E+00 0.2700E+00 0.2754E+00 0.2802E+00 0.2843E+00 0.2878E+00 0.2907E+00 0.2931E+00 0.2949E+00 0.2964E+00 0.2975E+00 0.2983E+00 0.2988E+00 0.2993E+00 0.2995E+00 0.2997E+00 0.2998E+00 0.2999E+00 0.2999E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.3000E+00 0.2999E+00 0.2999E+00 0.2997E+00 0.2994E+00 0.2989E+00 0.2980E+00 0.2964E+00 0.2937E+00 0.2890E+00 0.2814E+00 0.2693E+00 0.2509E+00 0.2245E+00 0.1899E+00 0.1499E+00 0.1099E+00 0.7540E-01 0.4903E-01 0.3061E-01 0.1852E-01 0.1091E-01 0.6277E-02 0.3530E-02 0.1941E-02 0.1042E-02 0.5460E-03 0.2774E-03 0.1333E-03 0.5327E-04 Values of PDE component i = 3 0.1916E-09 0.1973E-08 0.1191E-07 0.5547E-07 0.2183E-06 0.7558E-06 0.2356E-05 0.6708E-05 0.1762E-04 0.4304E-04 0.9826E-04 0.2107E-03 0.4262E-03 0.8157E-03 0.1482E-02 0.2562E-02 0.4228E-02 0.6681E-02 0.1013E-01 0.1479E-01 0.2084E-01 0.2843E-01 0.3766E-01 0.4855E-01 0.6109E-01 0.7522E-01 0.9082E-01 0.1078E+00 0.1259E+00 0.1450E+00 0.1650E+00 0.1857E+00 0.2069E+00 0.2284E+00 0.2500E+00 0.2717E+00 0.2932E+00 0.3143E+00 0.3350E+00 0.3550E+00 0.3741E+00 0.3922E+00 0.4092E+00 0.4248E+00 0.4389E+00 0.4514E+00 0.4623E+00 0.4715E+00 0.4791E+00 0.4852E+00 0.4899E+00 0.4933E+00 0.4958E+00 0.4974E+00 0.4985E+00 0.4992E+00 0.4996E+00 0.4998E+00 0.4999E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.5000E+00 0.4998E+00 0.4994E+00 0.4982E+00 0.4949E+00 0.4861E+00 0.4652E+00 0.4218E+00 0.3479E+00 0.2497E+00 0.1517E+00 0.7801E-01 0.3480E-01 0.1387E-01 0.5057E-02 0.1711E-02 0.5392E-03 0.1581E-03 0.4283E-04 0.1060E-04 0.2350E-05 0.4530E-06 0.7161E-07 0.8102E-08 Final statistics for mf = 22: 103 steps, 292 res, 17 jacobians, rwork size = 5368, iwork size = 317 Final output is correct to within 0.12E+02 times local error tolerance. ****************************************************************************************** Run completed: 0 errors encountered