complex function catanh (z)
c april 1977 version.  w. fullerton, c3, los alamos scientific lab.
      complex z, ci, catan
      external catan
      data ci /(0.,1.)/
c
      catanh = -ci*catan(ci*z)
c
      return
      end