home *** CD-ROM | disk | FTP | other *** search
/ The Fred Fish Collection 1.5 / ffcollection-1-5-1992-11.iso / ff_progs / libs / matlab.lzh / MATLAB / MATLAB.LZH / Source / MatLab / WATAN.FOR < prev    next >
Encoding:
Text File  |  1991-04-13  |  552 b   |  29 lines
  1.       SUBROUTINE WATAN (XR, XI, YR, YI)
  2.       IMPLICIT NONE
  3. C
  4. C     Y = ATAN(X) = (I/2)*LOG ((I+X)/(I-X))
  5. C
  6.       DOUBLE PRECISION XR, XI, YR, YI
  7. C
  8.       DOUBLE PRECISION TR, TI
  9. C
  10. C
  11.       IF (XI.NE.0.0D0) GO TO 10
  12.       YR = DATAN2 (XR, 1.0D0)
  13.       YI = 0.0D0
  14.       RETURN
  15. C
  16. 10    CONTINUE
  17.       IF (XR.NE.0.0D0 .OR. DABS (XI).NE.1.0D0) GO TO 20
  18.       CALL ERROR (32)
  19.       RETURN
  20. C
  21. 20    CONTINUE
  22.       CALL WDIV (XR, 1.0D0+XI, -XR, 1.0D0-XI, TR, TI)
  23.       CALL WLOG (TR, TI, TR, TI)
  24.       YR = -TI/2.0D0
  25.       YI = TR/2.0D0
  26. C
  27.       RETURN
  28.       END
  29.