[交流] 微分方程组就解，求大神帮忙 已有2人参与

function dy=rigid(t,y) dy=zeros(4,1); ey=11.65e-25;Ny=2.01e+27;k1=5e-22; k2=5e-22;ry=1000;c=3e+8;n1=1.5;ee=8e-25; re=125;C=6e-25;n2=1.75;l1=0.0031;l2=0.0019;L=0.005;egsa=3.5e-25;eesa=1e-26;R=0.85; Ns0=1.29e24;l=(n1-1).*l1+(n2-1).*l2+L; tr=2*l/(c*(log(1/R)+0.04)); ta=3.4e-7; Pin=200;wp=0.1e-3;h=6.62e-34; vp=c/(1535e-9); Wp=2.*Pin.*exp(-280.*0.0031)./(pi.*h.*vp.*(wp).^2); dy(1)=ey.*Ny.*Wp-k1.*y(1).*(9.88e25-y(2))-k2.*y(1).*y(2)-y(1).*ry; dy(2)=k1.*y(1).*(9.88e25-y(2))-(2*y(2)-9.88e25).*(c./n1).*ee.*y(3)-y(2).*re-C.*(y(2)).^2; dy(3)=y(3)./tr.* (2.*ee.*(2.*y(2)-9.88e+25).*l1-2.*egsa.*y(4).*l2-2.*eesa.*(Ns0-y(4)).*l2-log(1./R)-0.04); dy(4)=-egsa.*c./n2.*y(3).*y(4)+(Ns0-y(4))./ta; end [T,Y] = ode45(@rigid,[0 5e-08],[0 0 0 1.29e24] ,5e-09) 结果中y（3）y(4).值一直不变，和理论不一致。图片是结果。上边是源程序。求大神看看有什么问题，谢谢了~ 发自小木虫Android客户端

回复此楼