!@>:k3DC&
dy(1)=y(2); 这行报错,不知什么原因。 Vb?_RE_H
function dy=thermallens(t,y) g5
y*-t
global n2 'f6!a5qC
n0=1.5897; Ex{;&UWm
A=0.5^2; e{.P2rnh
n2=n0*A;%thermaload(1000);% 6lwWFR+k
dy=zeros(2,1); H+Q_%%[N
dy(1)=y(2); t<cWMx5ra
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);