'v.i' 6
dy(1)=y(2); 這行報錯,不知什么原因。 r3oAP[+n
function dy=thermallens(t,y) 9^Q:l0|
global n2 5n,?>>p$
n0=1.5897; `QV}je
A=0.5^2; L;=:OX0
n2=n0*A;%thermaload(1000);% TIbqUR
dy=zeros(2,1); c#M'Mye
dy(1)=y(2); ]q DhGt
dy(2)=(-n2*y(1)*(1+y(2)^2))/(n0-n2*y(1)^2/2);