SUBROUTINE EXPELAST ( Y, X, XBOX, Z, P, V, SQRTFZ, BETA, DELTA, VLY, VLX, VLZ, SIGW2, N, NX, NZ, NOBS, NULLX, NB, NVAZ, NLAMY, NUMLX, NUMLZ, NUMXEQZ, EXPY, EXPDERX, EXPELAX, EXPDERZ, EXPELAZ, R2E, PLIM, EYMEAN )
Argument Definitions (+ indicates altered content)
DOUBLE PRECISION   BETA(*),    DELTA(*),  +EXPDERX(*)
DOUBLE PRECISION  +EXPDERZ(*),            +EXPELAX(*)
DOUBLE PRECISION  +EXPELAZ(*),            +EXPY(*),   +EYMEAN
DOUBLE PRECISION   P(*),      +PLIM,      +R2E,        SIGW2
DOUBLE PRECISION   SQRTFZ(*),  V(*),       VLX(*),     VLY
DOUBLE PRECISION   VLZ(*),     X(N,NX),    XBOX(N,NX)
DOUBLE PRECISION   Y(N),       Z(N,NZ)
INTEGER            N,          NB,         NLAMY,      NOBS
INTEGER            NULLX(*),   NUMLX(*),   NUMLZ(*),   NUMXEQZ(*)
INTEGER            NVAZ,       NX,         NZ
Description
Version with the formulas involving P and V. 
Input arguments : (y,...,numxeqz) (see Definitions in Subroutine FCN). 
Output arguments: (expy,...,eymean) (see Definitions below).
Source file:expelast.f
Intrinsic Functions Called
DOUBLE PRECISION   DABS,       DEXP,       DSQRT
External Functions and Subroutines Called
DOUBLE PRECISION   ZLAM
SUBROUTINE         NDTR
Local Variables (+ indicates altered content)
DOUBLE PRECISION  +AVLY,      +BINF,      +BIS,       +BSS
DOUBLE PRECISION  +BSUP,      +BUS,       +CONST,      DENS
DOUBLE PRECISION  +DERX,      +DERZ,      +DEVEY,     +DEVY
DOUBLE PRECISION  +ELAX,      +ELAZ,      +FACTORX,   +FACTORZ
DOUBLE PRECISION  +FCT,       +PROB,      +PSF,       +PV
DOUBLE PRECISION  +RJAC,      +SDY2,      +SIGW,      +SIGY
DOUBLE PRECISION  +SQ2PI,     +SQSW,      +SUMEY2,    +SUMY2
DOUBLE PRECISION  +SUMYEY,    +TGAUSS(32),            +UPLIMY
DOUBLE PRECISION  +W,         +WGAUSS(32),            +WLOWER
DOUBLE PRECISION  +WPV,       +WUPPER,    +XB,        +XBW
DOUBLE PRECISION  +YMAX,      +YMEAN,     +YT,        +YTPOWER
INTEGER           +I,         +IK,        +IX,        +J
INTEGER           +K,         +NGROUP