From: pallav pal on
I want to get the isocline surface of the 3D system of differential equations
dx/dt=r*(1-(x+y)/K)*x-lambda*x*y-theta1*x*z/(a+x);
dy/dt=lambda*x*y-m1*y*z-mu*y;
dz/dt=theta2*x*z/(a+x)-m2*y*z-delta*z;
where r=3;k=45;lambda=.006;m1=.05;
a=15;mu=.24;delta=.09;theta1=.5;theta2=.2;m2=.04;

Please supply me the detailed matlab code for finding the all isocline surfaces in one fig.

Please note:
from 1st equation-
z=(a+x)*(r*(1-(x+y)/K)-lambda*y)/theta1
from 2nd equation-
z=(lambda*x-mu)/m1;
from 3rd equation-
y=(theta2*x/(a+x)-delta)/m2;
Please give me the matlab code