8.2 p170 (プログラムリスト2)

170ページのプログラムで - ページの計算実習2をやってみよう。

% ./p170
 EULER       : METHOD=1
 HEUN        :       =2
 RUNGE-KUTTA :       =3
 INPUT : METHOD=
1
 INPUT : N=
2
 INPUT : A(  1,  1)=
0
 INPUT : A(  1,  2)=
1
 INPUT : A(  2,  1)=
-1
 INPUT : A(  2,  2)=
0
 INPUT : T0=
0
 INPUT : X0(  1)=
1
 INPUT : X0(  2)=
1
 INPUT : H=
0.05
 INPUT : TMAX=
20
 T=  5.0000E-02
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  1.0500E+00   2  9.5000E-01
 T=  1.0000E-01
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  1.0975E+00   2  8.9750E-01
 T=  1.5000E-01
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  1.1424E+00   2  8.4262E-01
(中略)
 T=  1.9900E+01
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  2.2530E+00   2 -5.7165E-01
 T=  1.9950E+01
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  2.2244E+00   2 -6.8430E-01
 T=  2.0000E+01
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  2.1902E+00   2 -7.9552E-01
%

これは Euler 法で解いたものだが、テキストには結果が数値で与えられていないので、 グラフを描いて、172ページの図と比べてみよう。