function dx=diffeq(t,x,eta,kappa,nu,omega,phi,alpha,beta,gamma,delta)

n=x(1);
m=x(2);
f=x(3);
c=x(4);

dn=0;
df=alpha*eta*(m-f);
dm=beta*kappa*n-f*c+gamma*f-m;
dc=nu*f*m-omega*n-delta*phi*c;

dx=zeros(4,1);
dx(1)=dn;
dx(2)=dm;
dx(3)=df;
dx(4)=dc;