matlab sol函数,sol=bvp4c(@f,@fsbc,solinit)解决有上下限的ODE函数

function Pout=fiberlaser_twoend

global R1 R2 Pp1 sigma_ap sigma_ep sigma_as sigma_es sigma_s ...

gamma_p N alpha_p alpga_s Pssat Ppsat Ppsp Ppsv mu kappa elta

%参数设置

lambda_s=1100*1e-9;

lambda_p=974*1e-9;

tau=0.9e-3;

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

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

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

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

A_c=3.1416e-10;

N=5.5351e+025;

alpha_p=2e-5*1e2;

alpha_s=4e-6*1e2;

gamma_s=0.82;

gamma_p=0.0024;

R1=0.99;

R2=.035;

L=40;

%物理常数及中间过程计算

c=3e8;

h=6.626e-34;

nu_s=c/lambda_s;

nu_p=c/lambda_p;

Pssat=h * nu_s * A_c/(gamma_s * (sigma_es+sigma_as) * tau);

Ppsat=h * nu_p * A_c/(gamma_p * (sigma_ep+sigma_ap) * tau);

%抽运光功率