bief_erf.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\bief_erf.f
00002 !
00041                     DOUBLE PRECISION FUNCTION BIEF_ERF
00042 !                   **********************************
00043 !
00044      &(X)
00045 !
00046 !***********************************************************************
00047 ! BIEF   V6P3                                   07/03/2013
00048 !***********************************************************************
00049 !
00050 !
00051 !
00052 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00053 !| X              |-->| OPERAND
00054 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00055 !
00056 !     USE BIEF
00057 !
00058       IMPLICIT NONE
00059       INTEGER LNG,LU
00060       COMMON/INFO/LNG,LU
00061 !
00062 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00063 !
00064       DOUBLE PRECISION, INTENT(IN) :: X
00065 !
00066 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00067 !
00068       DOUBLE PRECISION BIEF_GAMMP
00069       EXTERNAL         BIEF_GAMMP
00070 !
00071 !-----------------------------------------------------------------------
00072 !
00073       IF(X.LT.0.D0) THEN
00074         BIEF_ERF = -BIEF_GAMMP(0.5D0,X**2)
00075       ELSE
00076         BIEF_ERF =  BIEF_GAMMP(0.5D0,X**2)
00077       ENDIF
00078 !
00079 !-----------------------------------------------------------------------
00080 !
00081       RETURN
00082       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0