数学建模社区-数学中国

标题: 请教这个mathematica程序什么问题？ [打印本页]

作者: bancal 时间: 2006-4-20 17:17
标题: 请教这个mathematica程序什么问题？

请教一下：这个实现欧拉方法mathematica程序怎么图形显示不出来？谢谢！

euler[f_, {t_, t0_, t1_, h_}, {x_, x0_}] =
Module[{ti, xi, graph}, ti = t0; xi = x0; graph = {{t0, x0}};
While[ti < t1, xi = N[xi + hf /. {t -> ti, x -> xi}]; ti = ti + h;
AppendTo[graph, {ti, xi}]]; ListPlot[graph] ]
pic2 = euler[f[t, p], {t, 0, 15, 0.5}, {p, 20}]

作者: madio 时间: 2006-4-20 23:56

估计你没有注意错误提示，里面显然告诉你hf不是一个浮点数。你的函数参数里面没有它，应该改成h，下面是修改过的代码，可以显示出图！

[attach]2820[/attach]

请教这个mathematica程序什么问题？.jpg (22.44 KB, 下载次数: 187)

请教这个mathematica程序什么问题？

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)