哪位高手知道这篇文章里面模型的编程方法呀，在国外一篇文献上也看到这一模型，但都没有程序，直接给出结果了。本人也想用用这个模型的改版形式。

请高手指点！ 不甚感激。