program randr
	real::random@
    do i=1,100
      r1=random@()
      r2=random@()
      z=sqrt(-2.*log(r1))*cos(2.*3.1415*r2)
      print *,z
    enddo
end program