function??Pout?=?fiberlaser_twoend
????global?R1?R2?Ppl?Ppr?sigma_ap?sigma_ep?sigma_as?sigma_es??...
????????gamma_p?gamma_s?N?alpha_p?alpha_s?Pssat?Ppsat?Ppsp?Ppsv?mu?kappa?elta
????
????%參數設置
????lambda_s?=?1100?*?1e-9;
????lambda_p?=?974?*?1e-9;
????tau?=?0.8e-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?=.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）;
????
????%抽運光功率設置
????Ppl?=?50;
????Ppr?=?50;
????
???%%?%‘端面抽運的光纖激光器邊值問題數值求解
????OPTION?=?bvpset（‘Stats‘‘ON‘）;
????solinit?=?bvpinit（linspace（0L10）[Ppl?Ppr?30?Ppr]）;%%%%沒看懂
????sol?=?bvp4c（@f@fsbcsolinit）;%給出速率方程?邊值條件?
????%%
????
????%數值計算結果分析和顯示nz就是N2/N
????x?=?[sol.x];
????y?=?[sol.y];
????nz?=?[（sigma_ap/（sigma_ap+sigma_ep）*（y（1:）+y（2:））/Ppsat+...
????????sigma_as/（sigma_as+sigma_es）*（y（3:）+y（4:））/Pssat）./...
????????（（y（1:）+y（2: