求助ode45求解问题！

sindy209

clear all
lb=inline('[xa(2)+sin(2*xa(1))*xa(1)^2+2sin(xa(1))*xa(3)+0.5*xa(1)*sin(t);gama+5*xa(1)+(6*xa(2)^2)/(1+2*xa(1)^2)+(1-cos(xa(2)))*xa(3)+0.5*xa(1)*cos(t);-xa(3)+xa(1)^2+0.5;0.5*q;-xa(5)+xa(1)^4+0.625]','t','x');

p=-1.5*xa(:,1)-xa(:,4).*xa(:,1).*(xa(:,1).^4+2*xa(:,5)+1);
a=-1.5-xa(:,4).*(5*xa(:,1)^4+2*xa(:,5)+1);
b=-xa(:,1)^5-2*xa(:,1).*xa(:,1);
c=-2*xa(:,1).*xa(:,4);
r1=xa(:,1).^2;
r2=xa(:,1).^2+xa(:,2).^2;
q1=xa(:,1).^4+xa(:,1).^2.*2xa(:,5)+xa(:,1).^2-0.001*xa(:,4);
q2=q1+(xa(:,2)-p).^2*(1+a^2)*(2*xa(:,5)+1)+p^2*(r2^2+(a*r1)^2)
gama=-2*p+(-1.5-xa(:,4).*(5*xa(:,1)^4+2*xa(:,5)+1))*xa(:,2)+0.5*(-xa(:,1)^5-2*xa(:,1).*xa(:,1)).*q2+(-2*xa(:,1).*xa(:,4)).*xa(:,5)-xa(:,4).*p*(1+a^2)*(2*xa(:,5)+1)-xa(:,4).*q2*(r2^2+(a*r1)^2);
options=odeset('RelTol',1e-4,'AbsTol',1e-4);
[t1,xa]=ode45(lb,[0 2],[1,1,1,1,1,1],options);
figure(1)
plot(t2,xa(:,1),'k:')

高手帮忙看下错在那！谢谢了！

[ 本帖最后由 无水1324 于 2008-4-14 18:25 编辑 ]

无水1324

提示的错误是什么？