matlab求解复制动态方程代码

** matlab求解复制动态方程代码**
syms x y z o s a C F1 r2 r3 V a P1 P2 D1 D2 Y1 Y2 e B
f1=x(1-x)((-s-F1)y+eyz+(s-C+F1))
f2=y*(1-y)((r2+s)x+(V+a)z+(P1+o-D1-P2))
f3=z(1-z)(r3x*y+(Y1-Y2)*y+B-D1)

f11=diff(f1,x)
f21=diff(f1,y)
f31=diff(f1,z)
f12=diff(f2,x)
f22=diff(f2,y)
f32=diff(f2,z)
f13=diff(f3,x)
f23=diff(f3,y)
f33=diff(f3,z)
x=0,y=0,z=0
eval(f11),eval(f21),eval(f31),eval(f12),eval(f22),eval(f32),eval(f13),eval(f23),eval(f33)
x=0,y=0,z=1
eval(f11),eval(f21),eval(f31),eval(f12),eval(f22),eval(f32),eval(f13),eval(f23),eval(f33)
x=0,y=1,z=0
eval(f11),eval(f21),eval(f31),eval(f12),eval(f22),eval(f32),eval(f13),eval(f23),eval(f33)
x=1,y=0,z=0
eval(f11),eval(f21),eval(f31),eval(f12),eval(f22),eval