我已经定义了函数
function dzdt=odefun2222(t,z)
beta=1+exp(-t); delta=1+exp(-t); rho=1+exp(-t);
gama=exp(-t); theta=exp(-t);
f1=1/2+1/(t+1); f2=1/2+1/(t+1);
dzdt=zeros(6,1);
dzdt(1)=z(2)-f1*z(1);
dzdt(2)=(1/4+1/(t+1)-beta)*z(1)+gama*z(3)-f1*z(2)-f1*(z(1)^2);
dzdt(3)=z(4);
dzdt(4)=gama*z(1)-delta*z(3)+theta*z(5);
dzdt(5)=z(6)-f2*z(5);
dzdt(6)=(1/4+1/(t+1)-rho)*z(5)+theta*z(3)-f2*z(6)-f2*(z(5)^2);
end
我已经使用了命令
tspan = [0 20];
z0 = [0.001 0.001 0.001 0.001 0.001 0.001];
[t,z] = ode45(@(t,z) odefun2222(t,z), tspan, z0);
plot(t,z(:,1),'b',t,z(:,3),'r',t,z(:,5),'g')
我只得到了 z(3) 和 z(5) 的绘图，z(1) 根本没有绘图。为了在同一输出中绘制所有三个函数 z(1)、z(3)、z(5)，我该怎么做？