求解以下微分方程

dx/dt=a(y-h(x))
dy/dt=x-y+z
dz/dt=-by
其中a,b为常数，当x大于1时，h(x)=2,当小于1时，h(x)=4

例如数值解法[t,y]=ode45(@qwe,[0 20],[2 0 1]);

function q=qwe(t,y)
q=zeros(3,1);
if y(1)>1
q(1)=y(2)-2
else
q(1)=y(2)-4
end

q(2)=y(1)-y(2)+y(3)
q(3)=-y(2)

T98yIUCo.txt (125 Bytes) [求助]求解微分方程组时如果方程内含有分段函数怎么办