zs0hXxTY:
dy(1)=y(2); 这行报错,不知什么原因。 >GmN~"iJ
function dy=thermallens(t,y) 0lBat_<8
global n2 Z~_8P
n0=1.5897; >J>>\Y(p
A=0.5^2; oCbpK
n2=n0*A;%thermaload(1000);% O<>#>[
dy=zeros(2,1); 6W$rY] h!
dy(1)=y(2); !D7"=G}HD
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);