SED navigation bar go to SED home page go to Dataplot home page go to NIST home page SED Home Page SED Staff SED Projects SED Products and Publications Search SED Pages
Dataplot Vol 2 Auxiliary Chapter

PSIFN

Name:
    PSIFN (LET)
Type:
    Library Function
Purpose:
    Compute the scaled kth derivative of the digamma (or psi) function.
Description:
    The psifn function is defined as:

      psifn(x,k) = ((-1)**(K+1)/GAMMA(K+1))*PSI(X,K)

    where GAMMA is the gamma function and PSI(X,K) is the kth derivative of the digamma function. Enter HELP DIGAMMA for details on the digamma function.

    Note that this is the kth derivative of the digamma function, not the log gamma function. That is, for k=1, the trigamma function is computed, not the digamma function.

Syntax:
    LET <y> = PSIFN(<x>)             <SUBSET/EXCEPT/FOR qualification>
    where <x> is a number, variable or a parameter;
                  <y> is a variable or a parameter (depending on what <x> is) where the computed psifn values are stored;
    and where the <SUBSET/EXCEPT/FOR qualification> is optional.
Examples:
    LET A = PSIFN(1,2)
    LET X2 = PSIFN(X1,K)
Note:
    DATAPLOT uses the routine PSIFN from the SLATEC Common Mathematical Library to compute this function. SLATEC is a large set of high quality, portable, public domain Fortran routines for various mathematical capabilities maintained by seven federal laboratories.
Default:
    None
Synonyms:
    None
Related Commands:
    GAMMA = Compute the gamma function.
    DIGAMMA = Compute the digamma function.
    LOGGAMMA = Compute the log (to base e) gamma function.
    GAMMAI = Compute the incomplete Gamma function.
    GAMMAIP = Compute an alternate form of the incomplete gamma function.
    GAMMAIC = Compute the complementary incomplete Gamma function.
    GAMMAR = Compute the reciprocal gamma function.
    TRICOMI = Compute Tricomi's incomplete gamma function.
    BETA = Compute the Beta function.
Reference:
    "Handbook of Mathematical Functions, Applied Mathematics Series, Vol. 55", Abramowitz and Stegun, National Bureau of Standards, 1964 (chapter 6).
Applications:
    Special Functions
Implementation Date:
    1998/5
Program:
    YLIMITS 0 100
    MULTIPLOT 2 2
    MULTIPLOT CORNER COORDINATES 5 5 95 95
    TITLE AUTOMATIC
    PLOT PSIFN(X,1) FOR X = 0.1 0.1 3
    PLOT PSIFN(X,2) FOR X = 0.1 0.1 3
    PLOT PSIFN(X,3) FOR X = 0.1 0.1 3
    PLOT PSIFN(X,4) FOR X = 0.1 0.1 3
    END OF MULTIPLOT
    MOVE 50 97
    JUSTIFICATION CENTER
    TEXT SCALED PSI FUNCTIONS

    plot generated by sample program

Date created: 6/5/2001
Last updated: 4/4/2003
Please email comments on this WWW page to alan.heckert@nist.gov.