maple作图问题现在有一个微分方程组如下:diff(x,t) = x*(z/(0.7857142858+z)-0.1169737892-0.30599022e-1*y),diff(y,t) = y*(4*z/(0.75+z)-0.5),diff(z,t) = z*(1-z-x/(0.7857142858+z)-4*y/(0.75+z)).给定初值x0 = 0.4,y0 =0.1,z0 =0.1,请问
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/24 20:37:49
maple作图问题现在有一个微分方程组如下:diff(x,t) = x*(z/(0.7857142858+z)-0.1169737892-0.30599022e-1*y),diff(y,t) = y*(4*z/(0.75+z)-0.5),diff(z,t) = z*(1-z-x/(0.7857142858+z)-4*y/(0.75+z)).给定初值x0 = 0.4,y0 =0.1,z0 =0.1,请问
maple作图问题
现在有一个微分方程组如下:diff(x,t) = x*(z/(0.7857142858+z)-0.1169737892-0.30599022e-1*y),diff(y,t) = y*(4*z/(0.75+z)-0.5),diff(z,t) = z*(1-z-x/(0.7857142858+z)-4*y/(0.75+z)).给定初值x0 = 0.4,y0 =0.1,z0 =0.1,请问作它的相图,最好用maple自带的命令.
maple作图问题现在有一个微分方程组如下:diff(x,t) = x*(z/(0.7857142858+z)-0.1169737892-0.30599022e-1*y),diff(y,t) = y*(4*z/(0.75+z)-0.5),diff(z,t) = z*(1-z-x/(0.7857142858+z)-4*y/(0.75+z)).给定初值x0 = 0.4,y0 =0.1,z0 =0.1,请问
restarte;
> eq:=
diff(x(t),t) = x(t)*(z(t)/(0.7857142858+z(t))-0.1169737892-0.30599022e-1*y(t)),
diff(y(t),t) = y(t)*(4*z(t)/(0.75+z(t))-0.5),
diff(z(t),t) = z(t)*(1-z(t)-x(t)/(0.7857142858+z(t))-4*y(t)/(0.75+z(t))),x(0) = 0.4,y(0) =0.1,z(0) =0.1;
> sss:= dsolve([eq],numeric):
> with(plots):odeplot(sss,[x(t),y(t),z(t)],t=0..30);
想要具体的t=0.01点数值,可以这样得到:
> sss(0.01);
[t = 0.01,x(t) = 0.399971436602128,y(t) = 0.0999705465997950,
z(t) = 0.0999778859911688]