DOUBLE PRECISION FUNCTION KSGENG(PW, VP, MVP)
#endif
#if !defined(CLEO_DOUBLE)
      REAL FUNCTION KSGENG(PW, VP, MVP)
#endif
 
C  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
C  Returns the error in energy for a track in the W parameterization.
C       sig = sqrt(Vx0(4,4))
C
C  PW        real array (read)
C            Track parameters in W form (at least 4-vector part)
C
C  VP        real array (read)
C            Covariance matrix
C
C  MVP       integer variable (read)
C            First dimension of VP (must be at least 4)
C  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>