康博虎的博客

    蒙特卡罗方法又称随机模拟方法或统计试验方法。它是通过不断产生随机数序列来模拟过程。

蒙特卡罗方法首先被美国科学家应用于原子弹的研制，目前这一方法已经广泛应用到物理学许多领域，像系统工程，科学管理，生物遗传，社会科学等这样一些学科领域也采用了这种研究方法。

    最简单的一个例子是用蒙特卡罗方法计算圆周率，方法是在正方形里做内切圆，在正方形内随机取两个数构成一个坐标，看它在圆内还是圆外，最终把在圆内的点数除以总点数乘以4得到结果。下面是用fortran语言编制的程序。

    program main
   real:: x,y,w,u,n,p,q
   p=0
   q=0
   read *,n
   do i=1,n
   call random(x)
   call random(y)
   u=x**2+y**2
   print *,u
   if(u>1) then
    p=p+1
   else
    q=q+1
   endif
   end do
   print *,p
   print *,q
   w=4*q/n
   write(*,*) w
end program main