各位大佬帮忙看一下哪里出问题了？
clc;clear;
tspan=[0 3];
no_lines=10;
a=1;
y0=linspace(0,a*0.95,no_lines); x0=sqrt(1-y0.^2);
x1=ones(1,length(x0))*2*a;y1=y0;
line([x0;x1],[y0;y1],'color','g','linewidth',2);hold on;
phi=linspace(0,2*pi,50);
x=a*cos(phi);y=a*sin(phi);
plot(x,y,'r','linewidth',3);
for lp=1:no_lines Y0=[x0(lp) y0(lp) -1 0]; options=odeset('RelTol',le-7,'AbsTol',[le-8]); [T,Y]=ode45(@ham,tspan,Y0,options); IN=inpolygon(Y(:,1),Y(:,2),[a*cos(phi)],[a*sin(phi)]); tmpx=Y(:,1);tmpy=Y(:,2); plot(tmpx(IN),tmpy(IN),'b','linewidth',2);hold on;
end
axis equal;
y0=-y0;x0=x0;
x1=ones(1,length(x0))*2*a;y1=y0;
line([x0;x1],[y0;y1],'color','g','linewidth',2);hold on;