SSYE&
dy(1)=y(2); 这行报错,不知什么原因。 |lm
function dy=thermallens(t,y) Aa[p7{e
global n2 UBs'3M
n0=1.5897; wQgW9546
A=0.5^2; "1$OPt5
n2=n0*A;%thermaload(1000);% 8G`fSac`
dy=zeros(2,1); 51W\ %aB
dy(1)=y(2); }i!hzkK#
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);