/kd6Yq(y
dy(1)=y(2); 这行报错,不知什么原因。 s\1c.
function dy=thermallens(t,y) ]M9r<x*
global n2 EtvYIfemr
n0=1.5897; #>\8m+h 9
A=0.5^2; {B6tGLt#bf
n2=n0*A;%thermaload(1000);% G`R2=bb8
dy=zeros(2,1); RT=(vq @
dy(1)=y(2); .[:*bo3
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);