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 / WDIV.FOR < prev    next >
Encoding:
Text File  |  1991-04-13  |  532 b   |  27 lines
  1.       SUBROUTINE WDIV (AR, AI, BR, BI, CR, CI)
  2.       IMPLICIT NONE
  3. C
  4. C C = A/B
  5. C
  6.       DOUBLE PRECISION AR, AI, BR, BI, CR, CI
  7. C
  8.       DOUBLE PRECISION S, D, ARS, AIS, BRS, BIS
  9. C
  10.       DOUBLE PRECISION FLOP
  11. C
  12. C
  13.       S = DABS (BR)+DABS (BI)
  14.       IF (S.EQ.0.0D0) CALL ERROR (27)
  15.       IF (S.EQ.0.0D0) RETURN
  16.       ARS = AR/S
  17.       AIS = AI/S
  18.       BRS = BR/S
  19.       BIS = BI/S
  20.       D = BRS**2+BIS**2
  21.       CR = FLOP ((ARS*BRS+AIS*BIS)/D)
  22.       CI = (AIS*BRS-ARS*BIS)/D
  23.       IF (CI.NE.0.0D0) CI = FLOP (CI)
  24. C
  25.       RETURN
  26.       END
  27.