参考张军教授的方法，计算结果有细微差异，不知道问题出在哪里。