C
FUNCTION G(Y1,Y2,PT2)
IMPLICIT REAL*8 (A-H,O-Z)
+#define BLANKET_SAVE
#include "hiparnt.inc"
DIMENSION F(2,7)
SAVE
XT=2.0*DSQRT(PT2)/HINT1(1)
X1=0.50*XT*(DEXP(Y1)+DEXP(Y2))
X2=0.50*XT*(DEXP(-Y1)+DEXP(-Y2))
- Z=DSQRT(1.D0-XT**2/X1/X2)
+ DUMMY = 1.D0-XT**2/X1/X2
+ if (DUMMY.LT.0.D0) DUMMY = 0.D0
+ Z=DSQRT(DUMMY)
SS=X1*X2*HINT1(1)**2
T=-(1.0-Z)/2.0
U=-(1.0+Z)/2.0