请教一个关于KTP参量允许角的问题

KTP的参量允许角是怎么算出来的啊，以及在Matlab里是怎么编程表示出图像。
我看见一个资料是采用姚建铨院士的理论模型，首先计算光波与光轴的夹角在最佳匹配方向处将失配量展开为Taylor级数并保留到二阶项再计算，可是我觉得编程里的未知项太多，没有办法编下去了，一个上午都在忙这个忙得焦头烂额的，哪位高手能帮一下我啊？谢谢！！！

你可以用SNLO软件啊