- 积分
- 13
- 注册时间
- 2008-11-28
- 仿真币
-
- 最后登录
- 1970-1-1
|
大家请看图,In[24]运行后没有输出任何结果,但当再运行OUt[24]之后却有内容显示,这是怎么回事?附上代码,麻烦各位帮忙解释一下。- bf[a_, b_, h_, x0_: - 1, y0_: - 1, z0_: - 1] :=
- Module[{t, k1, v1, m1, k2, v2, m2, k3, v3, m3, k4, v4, m4, x = x0,
- y = y0, z = z0},
- For[data = {}; t = a, t <= b, t += h, k1 = -10*(x - y);
- v1 = 28*x - y - x*z; m1 = (-8/3)*z + x*y;
- k2 = -10*((x + h/2.0*k1) - (y + h/2.0*v1));
- v2 = 28*(x + h/2.0*k1) - (y +
- h/2.0*v1) - (x + h/2.0*k1)*(z + h/2.0*m1);
- m2 = (-8/3)*(z + h/2.0*m1) + (x + h/2.0*k1)*(y + h/2.0*v1);
- k3 = -10*((x + h/2.0*k2) - (y + h/2.0*v2));
- v3 = 28*(x + h/2.0*k2) - (y +
- h/2.0*v2) - (x + h/2.0*k2)*(z + h/2.0*m2);
- m3 = (-8/3)*(z + h/2.0*m2) + (x + h/2.0*k2)*(y + h/2.0*v2);
- k4 = -10*((x + h*k3) - (y + h*v3));
- v4 = 28*(x + h*k3) - (y + h*v3) - (x + h*k3)*(z + h*m3);
- m4 = (-8/3)*(z + h*m3) + (x + h*k3)*(y + h*v3);
- x = x + h/6.0*(k1 + 2*k2 + 2*k3 + k4);
- y = y + h/6.0*(v1 + 2*v2 + 2*v3 + v4);
- z = z + h/6.0*(m1 + 2*m2 + 2*m3 + m4);
- data = Append[data, {t, x, y, z}]]]
- bf[0, 1, 0.01]
复制代码 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|