关于参数变化的常微分方程求解

spriner

需要模拟一个脉冲微分方程组,无脉冲方程简要为: dx/dt=r(K-x)-axy, dy/dt=qaxy-(m+d)y..
现在需要进行模拟，其中的r,K,a,q,d均为常数，而m的情况比较复杂，当x/y小于1时，m=0.1*(y/x-1)；当x/y大于等于1时，m=0. 求助如何在编写DHTI时候体现参数m的变化，我写了一个，但有些问题，如下：
function DHTI=Model20120326000(t,x)
global r m a K q d;
if x(1)/x(2)<1;
  m=0.1*((x(2)/x(1))-1);
else
    m=0;
end
DHTI=[r*(K-x(1))-a*x(1)*x(2);
           q*a*x(1)*x(2)-x(2)*(m+d);