跪求这个程序怎么改

daoying

hold on; box on;xlim([0,6]);
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=[0.01,0.1,1];
    [t,x]=ode45('hp',[0,200],x0,options,r,k,a1,b1,a2,b2,d1,d2);
    x0=x(end,:);
    for n=1:N;
        [t,x]=ode45('hp',[0,2],x0,options,r,k,a1,b1,a2,b2,d1,d2);
        x0=x(end,:);
        xd=x(:,2);
        [M,IK]=min(xd);
        if 1<IK &IK<length(xd);
            Tn=[Tn;M];
        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);