作业帮 > 数学 > 作业

maple作图问题现在有一个微分方程组如下:diff(x,t) = x*(z/(0.7857142858+z)-0.11

来源:学生作业帮 编辑:大师作文网作业帮 分类:数学作业 时间:2024/10/01 10:29:30
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.11
estarte;
> 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]