bief_gammp.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\bief_gammp.f
00002 !
00046                    DOUBLE PRECISION FUNCTION BIEF_GAMMP
00047 !                  ************************************
00048 !
00049      &(A,X)
00050 !
00051 !***********************************************************************
00052 ! BIEF   V6P3                                   07/03/2013
00053 !***********************************************************************
00054 !
00055 !
00056 !
00057 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00058 !| A              |-->| PARAMETER
00059 !| X              |-->| OPERAND
00060 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00061 !
00062 !     USE BIEF
00063 !
00064       IMPLICIT NONE
00065       INTEGER LNG,LU
00066       COMMON/INFO/LNG,LU
00067 !
00068 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00069 !
00070       DOUBLE PRECISION, INTENT(IN) :: X,A
00071 !
00072 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00073 !
00074       DOUBLE PRECISION GAMMCF,GLN
00075 !
00076 !-----------------------------------------------------------------------
00077 !
00078       IF(X.LT.A+1.D0) THEN
00079 !       Use the series representation
00080         CALL BIEF_GSER(BIEF_GAMMP,A,X,GLN)
00081       ELSE
00082 !       Use the continued fraction representation and take its complement
00083         CALL BIEF_GCF(GAMMCF,A,X,GLN)
00084         BIEF_GAMMP = 1.D0 - GAMMCF
00085       ENDIF
00086 !
00087 !-----------------------------------------------------------------------
00088 !
00089       RETURN
00090       END

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