8.1 p164 (プログラムリスト1)

164ページのプログラムで、166-167ページの計算実習1をやってみよう。

% ./p164
 LINEAR SYSTEM        : NFUNC=1
 VAN DER POL EQUATION :      =2
 EXERCISE-2           :      =3
 INPUT : NFUNC=
1
 EULER       : METHOD=1
 HEUN        :       =2
 RUNGE-KUTTA :       =3
 INPUT : METHOD=
1
 INPUT : N=
2
 INPUT : T0=
0
 INPUT : X0(  1)=
1
 INPUT : X0(  2)=
2
 INPUT : H=
0.1
 INPUT : TMAX=
2
 T=  1.0000E-01
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  1.0000E+00   2  1.9000E+00
 T=  2.0000E-01
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  9.8000E-01   2  1.7900E+00
(中略)
 T=  1.9000E+00
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  2.5576E-01   2  3.9084E-01
 T=  2.0000E+00
   I     X(I)      I     X(I)      I     X(I)      I     X(I)      I     X(I)
   1  2.3162E-01   2  3.5320E-01
%

これは $x_1(0)=1$ , $x_2(0)=2$ の場合に、 Euler法で計算したものだが、167ページの数表と見比べてどうだろうか？