euler2a integrator ncase = 1 neqn = 22 nsteps = 1000 t u(1) u(im) u(n-1) 0.0 298.00 298.00 298.00 33.3 1987.96 1002.15 394.17 66.7 1991.05 1218.09 622.84 100.0 1992.55 1345.04 834.65 133.3 1993.72 1447.20 1015.31 166.7 1994.70 1533.03 1168.10 euler2b integrator ncase = 2 neqn = 22 nsteps = 1000 t u(1) u(im) u(n-1) 0.0 298.00 298.00 298.00 33.3 1987.93 1002.11 394.14 66.7 1991.00 1217.92 622.58 100.0 1992.54 1344.68 834.02 133.3 1993.68 1446.70 1014.43 166.7 1994.65 1532.29 1166.78 rkc4a integrator ncase = 3 neqn = 22 nsteps = 1000 t u(1) u(im) u(n-1) 0.0 298.00 298.00 298.00 33.3 1987.96 1002.16 394.17 66.7 1991.05 1218.10 622.85 100.0 1992.55 1345.05 834.66 133.3 1993.72 1447.20 1015.32 166.7 1994.70 1533.03 1168.10 rkc4b integrator ncase = 4 neqn = 22 nsteps = 1000 t u(1) u(im) u(n-1) 0.0 298.00 298.00 298.00 33.3 1987.94 1001.97 394.06 66.7 1990.90 1217.84 622.46 100.0 1992.54 1344.63 833.92 133.3 1993.70 1446.55 1014.16 166.7 1994.65 1532.15 1166.53 rkc45a integrator ncase = 5 neqn = 22 nsteps = 1000 t u(1) u(im) u(n-1) 0.0 298.00 298.00 298.00 33.3 1987.96 1002.16 394.17 66.7 1991.05 1218.10 622.85 100.0 1992.55 1345.05 834.66 133.3 1993.72 1447.20 1015.32 166.7 1994.70 1533.03 1168.10 rkc45b integrator ncase = 6 neqn = 22 nsteps = 1000 t u(1) u(im) u(n-1) 0.0 298.00 298.00 298.00 33.3 1988.01 1002.16 394.17 66.7 1991.05 1217.81 622.40 100.0 1992.54 1344.61 833.89 133.3 1993.73 1446.47 1014.02 166.7 1994.70 1532.10 1166.45