如题，程序已经编好了
这是jixiao.m里面的
function dy=设计(t,y,R,M,A,E,S)
dy=zeros(2,1);
dy(1)=y(1)*(1-y(1))*(y(2)*(S-A-M)-(R-M));
dy(2)=y(2)*(1-y(2))*(y(1)*(A+M)-E);
end
这是主程序.m里面的
p0=0.2,q0=0.5,M=0.2,R=0.25,A=0.32,E=0.48,S=0.8;
[T,y]=ode45('jixiao',[0:40],[p0,q0],[],M,R,A,E,S);
y1=y(:,1);
plot(T,y1,'b*');
xlabel('t');
ylabel('p');
hold on;
p0=0.4,q0=0.5,M=0.2,R=0.25,A=0.32,E=0.48,S=0.8;
[T,y]=ode45('jixiao',[0:40],[p0,q0],[],M,R,A,E,S);
y1=y(:,1);
plot(T,y1,'b0');
xlabel('t');
ylabel('p');
hold on;
title('演化示意图');
legend('p0=0.2','p0=0.4','p0=0.6','p0=0.8');
但是运行主程序的时候总是出现Error: Missing variable or function.不知道是哪里的问题，求助诸位大神啊，急求