C-----------------------------------------------------------------------
C IMSL Name:  CWPL/ZWPL (Single/Double precision version)
C
C Purpose:    Evaluate the Weierstrass' P function in the lemniscatic
C             case for complex argument with unit period
C             parallelogram.
C
C Usage:      CWPL(Z)
C
C Example 1:
C                                  Declare variables
      INTEGER    NOUT
      COMPLEX    CWPL, VALUE, Z
      EXTERNAL   CWPL, UMACH
C                                  Compute
      Z     = (0.25, 0.25)
      VALUE = CWPL(Z)
C                                  Print the results
      CALL UMACH (2, NOUT)
      WRITE (NOUT,99999) Z, VALUE
99999 FORMAT (' CWPL(', F6.3, ',', F6.3, ') = (',
     &     F6.3, ',', F6.3, ')')
      END
C  CWPL( 0.250, 0.250) = ( 0.000,-6.875)