IooNb:(
dy(1)=y(2); 这行报错,不知什么原因。 |~<N -~.C
function dy=thermallens(t,y) xE-7P|2
global n2 *{o UWt
n0=1.5897; ^JR;epVJ
A=0.5^2; /b;K
n2=n0*A;%thermaload(1000);% PNeh#PI6)
dy=zeros(2,1); _#_
E^!
dy(1)=y(2); _e!F~V.
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);