seirs模型微分方程求解

function dy= ey1(t,y) %UNTITLED2 此处显示有关此函数的摘要 % 此处显示详细说明 b=1; w=0.3; r=0.35; a=0.2; dy=zeros(4,1); dy(1)=-b*y(1)*y(3)+a*y(4); dy(2)=b*y(1)*y(3)-w*y(2); dy(3)=w*y(2)-r*y(3); dy(4)=r*y(3)-a*y(4); end