运行以上的程序，当Nstep增大时，计算得到的π值会更接近于真实值，但由于每次开始计算时用的“种子”不同，Nstep取相同值时得到的结果也会有所波动。

　　1.9 参考文献
　　[1] 计算物理，杨清建，上海科学技术出版社，1988.5 第一版，p133-149
　　[2] 蒙特卡罗方法及应用，尹增谦等，物理与工程，Vo1.12，No.3，2002
　　[3] Compytational Physics, J.M.Thijssen, Cambridge University Press, First Published 1999, p271
　　[4] 计算物理讲义，罗成林，南京师范大学物科院，2002.9版，p61-72
　　[5] 大量粒子系统的课程教学体系探讨，王彬，物理与工程，Vol.11，No.1，2001
　　[6] 热学，秦允豪，高等教育出版社，2004.6，第2版，p261-264
　　[7] 数值分析，孙志忠、袁慰平等编，东南大学出版社，2002.2，第2版，p183-232