#Diff eq.
dM/dt=kappa*M^2*P-M
dP/dt= -kappa*M^2*P-P+1


# model parameters
param kappa=.5

done