+++ /dev/null
-*
-* $Id$
-*
-* $Log$
-* Revision 1.1.1.1 1996/04/01 15:01:48 mclareni
-* Mathlib gen
-*
-*
-#include "gen/pilot.h"
-#if defined(CERNLIB_DOUBLE)
- FUNCTION DWIG9J(A,B,C,P,Q,R,X,Y,Z)
- IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-
- PARAMETER (R1 = 1, HF = R1/2)
-
- IA=NINT(2*A)
- IB=NINT(2*B)
- IC=NINT(2*C)
- IP=NINT(2*P)
- IQ=NINT(2*Q)
- IR=NINT(2*R)
- IX=NINT(2*X)
- IY=NINT(2*Y)
- IZ=NINT(2*Z)
-
- H=0
- IF(IA .LT. 0 .OR. IB .LT. 0 .OR. IC .LT. 0 .OR.
- 1 IP .LT. 0 .OR. IQ .LT. 0 .OR. IR .LT. 0 .OR.
- 2 IX .LT. 0 .OR. IY .LT. 0 .OR. IZ .LT. 0) GO TO 99
- J0=MAX(ABS(IA-IZ),ABS(IY-IP),ABS(IB-IR))
- J1=MIN(IA+IZ,IY+IP,IB+IR)
- S=0
- V=(-1)**J0
- DO 1 J = J0,J1
- AJ=HF*J
- H=DWIG6J(A,P,X,Y,Z,AJ)
- IF(H .NE. 0) H=H*DWIG6J(B,Q,Y,P,AJ,R)
- IF(H .NE. 0) H=H*DWIG6J(C,R,Z,AJ,A,B)
- S=S+V*(AJ+HF)*H
- 1 V=-V
- H=2*S
- 99 DWIG9J=H
- RETURN
- END
-#endif