!-------------------------------
!       get gaussian probability
!       for given sigma:  P(x>sigma) =1-erf(x/sqrt(2.))
Program GaussProb

Real*8 x,y
x =3.
do
  write(*,'(a,$)') ' enter x='
  read(*,*) x
  y = 1.-erf(x/sqrt(2.d0))
  write (*,*) 'x=',x,' probability =',y
enddo

stop
end Program GaussProb