HbUadPr
dy(1)=y(2); 这行报错,不知什么原因。 Ts~)0
function dy=thermallens(t,y) 320g!r
global n2 p#T^o]+
n0=1.5897; F,l%SQCyj
A=0.5^2; (qky&}H
n2=n0*A;%thermaload(1000);% ~l@-gAyw
dy=zeros(2,1); 5(}Qg9%
dy(1)=y(2); Wt8=j1>
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);