matlab在激光器中的应用,请教关于MATLAB 中bvp4c函数 用于光纤激光器多模传输的仿真....

function Pout= DCFlaseroutpower_LP01    global  sigma_ap sigma_ep sigma_as sigma_es gamma_s01 gamma_s11 gamma_s21 gamma_s02 gamma_p N ...

alpha_p alpha_s   nu_s nu_p tau  a  b   i X1 X2 Y1 Y2 R1 R2 Ppl Ppr

%光纤激光器参数设置

lambda_s=1064*1e-9;

lambda_p=975*1e-9;

tau=0.8e-3;

sigma_ap=25e-21*1e-4;

sigma_ep=25e-21*1e-4;

sigma_as=6.4e-23*1e-4;

sigma_es=3.2e-21*1e-4;

N=4.6e+019*1e6;

alpha_p=2e-5*1e2;

alpha_s=4e-6*1e2;

a=14e-6;

b=250e-6;

NAcore=0.06;

R1=0.99;

R2=0.04;

L=10;

V=2*pi*a*NAcore/lambda_s;

U01=1.9912;

W01=sqrt(V^2-U01^2);

K=2*pi/lambda_s;

%物理常数

c=3e8;

h=6.626e-34;

nu_s=c/lambda_s;

nu_p=c/lambda_p;

%