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 / WLOG.FOR < prev    next >
Encoding:
Text File  |  1991-04-13  |  411 b   |  23 lines
  1.       SUBROUTINE WLOG (XR, XI, YR, YI)
  2.       IMPLICIT NONE
  3. C
  4. C Y = LOG(X)
  5. C
  6.       DOUBLE PRECISION XR, XI, YR, YI
  7. C
  8.       DOUBLE PRECISION T, R
  9. C
  10.       DOUBLE PRECISION PYTHAG
  11. C
  12. C
  13.       R = PYTHAG (XR, XI)
  14.       IF (R.EQ.0.0D0) CALL ERROR (32)
  15.       IF (R.EQ.0.0D0) RETURN
  16.       T = DATAN2 (XI, XR)
  17.       IF (XI.EQ.0.0D0 .AND. XR.LT.0.0D0) T = DABS (T)
  18.       YR = DLOG (R)
  19.       YI = T
  20. C
  21.       RETURN
  22.       END
  23.