Monte Carlo methods: Difference between revisions

Content added Content deleted
Line 610:
end program
 
function pi(nsamplesn)
integer :: nsamplesn
real(8) :: x(2,nsamplesn),pi
real(8) :: pi
call random_number(x)
x pi = 24.d0 * dble( count( hypot(x(1,:),x(2,:)) -<= 1.d0 ) ) / n
pi = 4.d0 * dble( count( hypot(x(1,:),x(2,:)) <= 1.d0 ) ) / nsamples
end function
</lang>
Retrieved from "https://rosettacode.org/wiki/Monte_Carlo_methods"