Monte Carlo methods: Difference between revisions
Content added Content deleted
No edit summary |
|||
Line 598:
100000000 3.14147
</pre>
{{works with|Fortran|2008 and later}}
<lang fortran>
program mc
integer :: n,i
real(8) :: pi
n=10000
do i=1,5
print*,n,pi(n)
n = n * 10
end do
end program
function pi(nsamples)
integer :: nsamples
real(8) :: x(2,nsamples)
real(8) :: pi
call random_number(x)
x = 2.d0*x - 1.d0
pi = 4.d0 * dble( count( hypot(x(1,:),x(2,:)) <= 1.d0 ) ) / nsamples
end function
</lang>
=={{header|Futhark}}==
| |||