各位大佬帮忙看一下哪里出问题了？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;