含外腔的半导体激光的理论研究通常也基于Lang-Kobayashi速率方程。Lang-Lobayashi速率方程由复光强方程和激光腔内的载流子密度方程组成，请问如何用如何用matlab语言的四阶-龙格库塔算法 解Lang-Kobayashi速率方程，有现成仿真程序最好，谢谢