请教各位高手问题，不胜感激！

我正做edfa,，想用matlab模拟gils方程。刚接触不久，现有几点疑问，求各位大哥大姐指点：

1。其pk正和负包括前向和后向的光，是不是对每个频率点都是这样，如波长间隔取1nm,从1520-1560划分。泵浦光是不是在求平均粒子束反转度的时候也要加上一起计算；

2。使用龙格-库塔法进行计算可行不？因其为一阶线性偏微分方程，从上一个点的值迭代出下一点。步长取多大合适？

3。在带入参数进行计算时alfa（吸收系数）和g（增益系数 ）是否要转化为倍数关系。若alfa取为db/m,在计算时用10^(0.1*alfa)是否正确？ alfa和g怎样取呢，是否在厂家给出的曲线频率点上大概取值？

4。基于以上假设，步长h取为0.2m ，我在计算时得到的第一个值在减小，很困惑。如在1550nm处，pk+=0, pk(1550+)=pin, pp(泵浦+)=pp;

恳请高手给我指点，非常感谢！

大家帮帮我这个菜鸟吧，是不是问题太ben了，有没有程序给我参考一下，邮箱lai-fan@163.com