Wolfram mathematica Mathematica 8使用Heun'；s方法/改进的Euler'；s法

上图是我的错误。我尝试使用Heun的方法，我的问题是：

1） 我想让它停在

y[5pi]

但它一直在运行。我可以操纵它，使它转到

y[5pi]

，但我想知道它为什么这样做

---

2）

y[k]

不是在k=pi、2pi、3pi等位置进行计算。

请以文本形式发布您的代码，因此任何想要回答的人都不需要重新键入。谢谢你。。。你忘了定义函数

y[x]

对不起，我是新来的！下面是我的代码清晰[x，y，h，k，FirstSlope，SecondSlope]；h=[Pi]；y[[Pi]]=0；dy[x，y][=（Cos[x]-3x^2y）/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]}表[{x[k]，y[k]}，{k[Pi]，5[Pi]}；MatrixForm[%]我不认为定义y[x_]是个问题——我认为它与π有关。清除[x，y，h，k，第一坡度，第一坡度，第二坡度]；h=1；y[1]=0.5；dy[x_，y_]=（yLog[x]）/Log[y]；Do[{x[k]=1+h*（k-1），FirstSlope=dy[x[k]，y[k]]，SecondSlope=dy[x[k]+h，y[k]+hFirstSlope]，y[k+1]=y[k]+（h*（FirstSlope+SecondSlope））/2}，{k，1，5}表[{x[k]，y[k]}，{k，1，5}]；例如，MatrixForm[%]，它工作非常好，不在注释中发布代码；没人会去读的。将其编辑到您的问题中，并使用提供的工具设置其格式。

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[%]