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 などで評価していません。