//Definition of the MAC differential equations. X(1)=x and X(2)=y.

function [dif]=mac(t,X,a,b,m,HD,WN,r,p)
    dif(1)=a*b*m*X(2)-a*b*m*X(1)*X(2)-r*X(1)
    dif(2)=a*X(1)-a*X(1)*X(2)+log(p)*X(2)
endfunction