附录六 求解模型的Matlab程序