ode45 matlab仿真程序,MATLAB的ode45函数的程序运行问题？

以下程序存在以Laser2为名的.m文件中。

function dy=Laser2(t,y)

Delta=2.8e-23;%激光增益介质受激发射截面面积

Deltags=1.1e-22;%饱和吸收体的基态吸收截面面积

Deltaes=1.2e-23;%饱和吸收体的激发态吸收截面面积

tg=2.3e-4;%饱和吸收体的恢复时间

ts=3.4e-6;

C=3e8;%真空中光速

P=0.2;%泵浦光的泵浦功率

L=0.047;%谐振腔的光学长度

lg=1.2e-3;%激光增益介质长度

ls=8e-4;

WL=808e-9;%泵浦光波长

V=C/WL;%泵浦光频率

h=6.63e-34;%普朗克常数

T0=0.9;%可饱和吸收体的初始透过率

tr=2*L/C;%光在谐振腔中往返一周所需时间

asc=0.2;%激光器的散射损耗，典型值为0.002cm-1

R=0.95;%输出镜反射率

Wp=5e-5;%泵浦光的光斑半径

Rp=P/(h*V*pi*(Wp^2)*(lg+ls));

dy=zeros(3,1);

dy(1)=y(1)/tr*(2*Delta*y(2)*lg-2*Deltaes*((log(1/T0))/(Deltags*(lg+ls))-y(3))*ls-(log(1/R)+2*asc*(lg+ls)));

dy(2)=Rp-Delta*C*y(1)*y(2)-y(2)/tg;

dy(3)=-Deltags*C*y(1)*y(3)+((log(1/T0))/(Deltags*(lg+ls))-y(3))/ts;

end

以下依次程序在command命令窗口运行。

[t,y]=ode45('Laser2',[0 5],[0,0,0]);

figure,plot(t,y(:,1),'r');

figure,plot(t,y(:,2),'b');

figure,plot(t,y(:,3),'g');

我的电脑性能不是太好，貌似运行好久也出不来，不知道哪位大侠能帮忙跑一下，把结果给我啦！！感激不尽啊，急等！谢谢了！