matlab中ode指令,matlab中使用ode45函数老是报错，请帮忙看下程序 - 仿真模拟 - 小木虫 - 学术 科研 互动社区...

function laser_New                %铒镱共掺光纤激光器的传输方程

global   sigma_12 sigma_21 sigma_13 sigma_56 sigma_65 Ner....

Nyb  alpha_p alpha_s gamma_p gamma_s h c nu_p nu_s   Ccr Acore tau_er tau_yb  lambda_s Es  Ep Pp0  Ps0  lambda_p ps0

%参量设置

lambda_p=980*1e-9;     %泵浦光波长1

lambda_s=1550*1e-9;    %激光波长

tau_er=10*1e-3;           %Er上能级粒子寿命

tau_yb=1.35*1e-3;          %Yb上能级粒子寿命

sigma_12=1.75*1e-25;    %铒离子的吸收截面面积

sigma_21=2.4*1e-25;    %铒离子的发射截面面积

sigma_13=1.97*1e-27;

sigma_56=2.19*1e-25;    %镱离子的吸收截面面积

sigma_65=3.14*1e-26;    %镱离子的发射截面面积

Ner=4.33*1e25;          %铒离子的浓度

Nyb=3.61*1e26;         %镱离子的浓度

alpha_p=0.15;        %双包层光纤对泵浦光的损耗

alpha_s=0.1;      %双包层光纤对激光的损耗

gamma_p=0.0058;       %泵浦光重叠因子

gamma_s=1;          %激光重叠因子

Ccr=1e-2