[y z w]=Dsolve('Dy=u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x','Dz=a*Dy+b','Dw=(-1/Y1)*Dz-(1/Y2)*Dy-m*s')
y =
-(t*m*s*Y1*Y2+t*a*u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x*Y2+t*b*Y2+u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x*t*Y1-C1*Y1*Y2)/Y1/Y2
z =
u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x*t+C3
w =
t*a*u*(1-x/xm)^n1*(1-p/pm)^n2*(1-s/sm)^n3*x+t*b+C2
|