program square integer, parameter :: count=100000 integer :: i,hit=0 real(kind(1d0)) :: x(2) do i=1,count call random_number(x) if ((x(1)**2+x(2)**2)<1d0) hit=hit+1 enddo write(*,101) hit/real(count,kind(1d0)) 101 format ('Ratio: ',f13.9) end program