0

ここに画像の説明を入力

Clear[x, y, h, k, FirstSlope, SecondSlope];
h = [Pi]; y[[Pi]] = 0;
dy[x_, y_] = (Cos[x] - 3 x^2 y)/x^3;
Do[{x[k] = [Pi] + h*(k - [Pi]),
   FirstSlope = dy[x[k], y[k]],
   SecondSlope = dy[x[k] + h, y[k] + h*FirstSlope], 
   y[k + [Pi]] = y[k] + (h*(FirstSlope + SecondSlope))/2}, {k, [Pi], 
   5[Pi]}] Table[{x[k], y[k]}, {k, [Pi], 5[Pi]}]; 
MatrixForm[%]

上の画像は私のエラーです。Heunの方法を使用しようとしていますが、私の問題は次のとおりです。

1) 止まってほしいのですy[5 Pi]が、進み続けます。に行くように操作できますがy[5 Pi]、正確にこれを行っている理由を知りたいです。

2) y[k]k=pi,2pi,3pi などで評価していません。

4

0 に答える 0