Fortran - vypocet pi

Dobrý den,
pokouším se vypočítat pí pomocí metody Monte Carlo, ale nějak jsem se do toho zamotal a nevím kudy kam. Zkoušek jsem to takto:

program vypocet_pi
integer :: n0, i
real xp, yp, y

call random_seed()
n0 = 0
do i =0, 1000000
call random_number(xp)
call random_number(yp)
y=sqrt (1.0 -xp**2)
if (yp<=) n0 = n0+1 ----------------------v tomto řádku dělám asi chybu, ale co tam je špatně?
if(mod(1,1000) == 0) write (*,*) i, 4.0*real(n0)/real(i),pi
end do

end program vypocet_pi

Budu rád za každou reakci, děkuji