jM07&o]D
dy(1)=y(2); 这行报错,不知什么原因。 MpCK/eiC
function dy=thermallens(t,y) drkY~!a
global n2 lKT<aYX
n0=1.5897; _X/`7!f
A=0.5^2; r0'6\MS13
n2=n0*A;%thermaload(1000);% xp68-&
dy=zeros(2,1); 8p0ZIrD%
dy(1)=y(2); l2|[
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);