请帮忙看下这个程序

gwei

r=ones(1,10);
[t,yp]=ode45('funyp',[0:10],1);
[t,ym]=ode45('funym',[0:10],1);
e=yp-ym;
[t,c]=ode45('fundc',[0:10],1);
[t,d]=ode45('fundd',[0:10],1);
u=c*r+d*yp;

function dy=funyp(t,y)
a=0.005;b=0.005;
r=ones(1,100);
dy=-a*y+b*r;

function dy=funym(t,y)
r=ones(1,10);
dy=-0.1*y+0.1*r;


function dy=fundc(t,y)
dy=-r*e;

function dy=fundd(t,y)
dy=-yp*e;
刚开始学matlab，请帮小弟看下哪有问题，多谢了！