1 SUBROUTINE zs1s2(ZRR, ZRI, S1R, S1I, S2R, S2I, NZ, ASCLE, ALIM,
17 DOUBLE PRECISION AA, ALIM, ALN, ASCLE, AS1, AS2, C1I, C1R, S1DI,
18 * s1dr, s1i, s1r, s2i, s2r, zeroi, zeror, zri, zrr, xzabs
20 DATA zeror,zeroi / 0.0d0 , 0.0d0 /
24 IF (s1r.EQ.0.0d0 .AND. s1i.EQ.0.0d0)
GO TO 10
25 IF (as1.EQ.0.0d0)
GO TO 10
26 aln = -zrr - zrr + dlog(as1)
32 IF (aln.LT.(-alim))
GO TO 10
33 CALL xzlog(s1dr, s1di, c1r, c1i, idum)
36 CALL xzexp(c1r, c1i, s1r, s1i)
41 IF (aa.GT.ascle)
RETURN