跪求这个程序怎么改
hold on; box on;xlim();N=100;
Tn=[];
a1=4;
a2=4;
b1=3;
b2=3;
d1=0.4;
d2=0.6;
k=1.5;
options=odeset('RelTol',1e-9);
for r=0:0.01:6;
x0=;
=ode45('hp',,x0,options,r,k,a1,b1,a2,b2,d1,d2);
x0=x(end,:);
for n=1:N;
=ode45('hp',,x0,options,r,k,a1,b1,a2,b2,d1,d2);
x0=x(end,:);
xd=x(:,2);
=min(xd);
if 1<IK &IK<length(xd);
Tn=;
end
end
plot(r,Tn,'k.','markersize',1);
Tn=[];
pause(0.01);
end
xlabel('\itc','Fontsize',16,'Fontname','Times new roman');
xlabel('\itx','Fontsize',16,'Fontname','Times new roman','Rotation',0);
运行时为
Error in ==> hhhh at 25
plot(r,Tn,'k.','markersize',1);
页:
[1]