o=@0Bd8
dy(1)=y(2); 这行报错,不知什么原因。 &D0suK#
function dy=thermallens(t,y) exT
O#*o
global n2 K=T]@ix$
n0=1.5897; d*Q:[RUf,
A=0.5^2; >oSNKE
n2=n0*A;%thermaload(1000);% golr,+LSo
dy=zeros(2,1); "h?;)Ye
dy(1)=y(2); 2NHuZ.af
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);