e
RjpR?!\
dy(1)=y(2); 这行报错,不知什么原因。 AyMMr_q
function dy=thermallens(t,y) Tfl4MDZb
global n2 DOm5 azO!>
n0=1.5897; V OViOD
A=0.5^2;
$dLPvN
n2=n0*A;%thermaload(1000);% G&Sp }
dy=zeros(2,1); Y+tXWN"8
dy(1)=y(2); A@0%7xm
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);