fox-li迭代法，用matlab编程解决折叠腔的的输出模式问题。

想请教一个问题，用matlab编程解决：求出四镜折叠腔经过多次迭代达到稳态输出时在输出镜上的振幅和相位分布，编程时用fox-li迭代法，而不是矩阵法。也就是最后画出这两个图来，一个是相位分布，另一个是振幅分布。急！！！我只会最简单的两个镜的平行平面腔，还是对称腔，我不明白在四镜折叠腔中，折叠角和腔镜的曲率半径这些参数怎么在amtlab编程中体现出来，也就是说怎么加进去。请愿意帮我解答的人最好给个现成的程序例子。