.!}hhiF,Z
dy(1)=y(2); 这行报错,不知什么原因。 )n=ARDd^e
function dy=thermallens(t,y) x>~p;z#VX
global n2 U4Nh
n0=1.5897; !eJCM`cp
A=0.5^2; MbXq`%
n2=n0*A;%thermaload(1000);% t!_x(u
dy=zeros(2,1); s^T+5E&}
dy(1)=y(2); y>\S@I
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);