c CONVEX DEFINE



      SUBROUTINE DIGRESP(OM,NZ,Z,NP,P,RESP1,RESP2)
      real*4 om,resp1,resp2	
      COMPLEX Z(20),P(30),T,B,O
      T=CMPLX(1.,0.)
      B=CMPLX(1.,0.)
      O=CMPLX(0.,OM)
      DO 10 I=1,NZ
         T=T*(O-Z(I))
   10 continue
      DO 20 I=1,NP
         B=B*(O-P(I))
   20 continue
      RESP1=REAL(T/B)
      RESP2=AIMAG(T/B)
      RETURN
      END