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