bief_erf.f

Go to the documentation of this file.

!                   **********************************
                    DOUBLE PRECISION FUNCTION bief_erf
!                   **********************************
!
     &(x)
!
!***********************************************************************
! BIEF   V6P3                                   07/03/2013
!***********************************************************************
!
!brief    Function erf (until all compilers have it).
!+        From 'Numerical Recipes' Chapter 6, and Tideway model.
!
!history  D J Evans-Roberts (HRW)
!+        23/11/1993
!+        V6P3
!+   First version sent by Michiel Knaapen (HRW) on 07/03/2013.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| X              |-->| OPERAND
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
!     USE BIEF
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION, INTENT(IN) :: X
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION BIEF_GAMMP
      EXTERNAL         bief_gammp
!
!-----------------------------------------------------------------------
!
      IF(x.LT.0.d0) THEN
        bief_erf = -bief_gammp(0.5d0,x**2)
      ELSE
        bief_erf =  bief_gammp(0.5d0,x**2)
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END