Zss `##
dy(1)=y(2); 这行报错,不知什么原因。 m<DiYxK
function dy=thermallens(t,y) "&,Gn#'FG
global n2 ZM~kc|&
n0=1.5897; 8<VO>WA>E
A=0.5^2; GT)63|
n2=n0*A;%thermaload(1000);% 5:o$]LkOWC
dy=zeros(2,1); *nPB+@f
dy(1)=y(2); A* =r~T5B
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);