+o ;}*
dy(1)=y(2); 这行报错,不知什么原因。 6T`F'Fk[
function dy=thermallens(t,y) @5=2+ M
global n2 )j_Y9`R
n0=1.5897; 8kRqF?rbj
A=0.5^2; q{c/TRp7
n2=n0*A;%thermaload(1000);% 0#/N ZO
dy=zeros(2,1); u)hr
dy(1)=y(2); mw!EDJ;'
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);