Monte Carlo methods: Difference between revisions
Content added Content deleted
No edit summary |
|||
Line 598: | Line 598: | ||
100000000 3.14147 |
100000000 3.14147 |
||
</pre> |
</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}}== |
=={{header|Futhark}}== |