DOUBLE PRECISION FUNCTION DFUN( N, xy )
C ----------------------------------------------------------------------------
C           Top level function call: integrated by DIVON4.
C           What function is actually called depends upon the value of
C           Function_type in the SUSFUN common block.
C
C           1 = Cross section for photino production
C
C         101 = Decay matrix element/phase space for decay to neutrino
C         102 = Decay matrix element/phase space for decay to electron
C ----------------------------------------------------------------------------
C
      IMPLICIT NONE
C
      INCLUDE 'SUSANA.INC'
C
      INTEGER N
      DOUBLE PRECISION xy(2)
      DOUBLE PRECISION suxsec, suswd1, suswd2
C
      IF (Function_type .EQ. 1) THEN
C
C --    Cross section, process 1.
        DFUN = SUXSEC(N,xy)
C
      ELSEIF (Function_type .EQ. 101) THEN
C
C --    Phase space for photino -> neutrino decay.
        DFUN = SUSWD1(N,xy)
C
      ELSEIF (Function_type .EQ. 102) THEN
C
C --    Phase space for photino -> electron decay.
        DFUN = SUSWD2(N,xy)
C
      ELSE
C
C --    Illegal call.
        WRITE (isustr(1),*) 'Illegal value of Function_type!'
        STOP
C
      ENDIF
C
      RETURN
      END