MINICERN/mathlib/gen/d/rfrdh264.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:22  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10 #if defined(CERNLIB_DOUBLE)
  11       FUNCTION DFRDH2(F,G,X,WS,IDIM,N)
  12 C
  13 #include "gen/imp64.inc"
  14 C
  15       CHARACTER*(*) NAME
  16       PARAMETER(NAME='DFRDH2')
  17 #endif
  18 #if !defined(CERNLIB_DOUBLE)
  19       FUNCTION RFRDH2(F,G,X,WS,IDIM,N)
  20 C
  21       CHARACTER*(*) NAME
  22       PARAMETER(NAME='RFRDH2')
  23 #endif
  24 C
  25       DIMENSION WS(IDIM,*)
  26
  27       SUM=F(X)
  28       DO 2 J = 1,N
  29     2 SUM=SUM+WS(J,N+4)*G(X,WS(J,N+3))*WS(J,N+1)
  30       GO TO 9
  31
  32 #if defined(CERNLIB_DOUBLE)
  33       ENTRY DFRDH3(H,WS,IDIM,N)
  34 #endif
  35 #if !defined(CERNLIB_DOUBLE)
  36       ENTRY RFRDH3(H,WS,IDIM,N)
  37 #endif
  38
  39       SUM=0
  40       DO 3 J = 1,N
  41     3 SUM=SUM+WS(J,N+4)*H(WS(J,N+3))*WS(J,N+1)
  42 #if defined(CERNLIB_DOUBLE)
  43     9 DFRDH2=SUM
  44 #endif
  45 #if !defined(CERNLIB_DOUBLE)
  46     9 RFRDH2=SUM
  47 #endif
  48       RETURN
  49       END