激光器速率方程matlab程序

% 微分方程  
function f=laser18（zp）
global input
global pp
wavelength1=975e-9;    
wavelength2=1086e-9;   
t=1e-3;                
aap=2.5e-20;         
aep=2.5e-20;        
aes=2.5e-21;         
aas=1.6e-23;      
A=pi*（29.7/2）*1e-8;            
n=8.93e19 ;          %掺杂浓度
     ap=5e-5; 
tp=0.0043 ;      %重叠因子
ts=0.8;