5 * Revision 1.1.1.1 1995/10/24 10:21:06 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.40 by S.Giani
12 FUNCTION DPDNUC(A,B,C)
15 C ORIGIN : H.FESEFELDT
17 #if !defined(CERNLIB_SINGLE)
18 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
21 C DPDK = SQRT(A*A+(B*B-C*C)**2/(A*A) - 2.0*(B*B+C*C))/2.0
27 ARG=A2+(B2-C2)**2/A2-2.0*(B2+C2)
32 DPDNUC = 0.5*SQRT(A2 + (B2-C2)**2/A2 - 2.0*(B2+C2))