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 / PYTHAG.FOR < prev    next >
Encoding:
Text File  |  1991-04-13  |  671 b   |  34 lines
  1.       DOUBLE PRECISION FUNCTION PYTHAG (A, B)
  2.       IMPLICIT NONE
  3. C
  4.       DOUBLE PRECISION A, B
  5. C
  6.       INCLUDE MATLAB$KOM:IOP.KOM
  7. C
  8.       DOUBLE PRECISION P, Q, R, S, T
  9. C
  10. C
  11.       P = DMAX1 (DABS (A), DABS (B))
  12.       Q = DMIN1 (DABS (A), DABS (B))
  13.       IF (Q.EQ.0.0D0) GO TO 20
  14.       IF (DDT.EQ.25) WRITE (WTE, 1)
  15. 1     FORMAT (' PYTHAG', 1P2D23.15)
  16.       IF (DDT.EQ.25) WRITE (WTE, 2) P, Q
  17. 2     FORMAT (1X, 1P2D23.15)
  18. C
  19. 10    CONTINUE
  20.       R = (Q/P)**2
  21.       T = 4.0D0+R
  22.       IF (T.EQ.4.0D0) GO TO 20
  23.       S = R/T
  24.       P = P+2.0D0*P*S
  25.       Q = Q*S
  26.       IF (DDT.EQ.25) WRITE (WTE, 2) P, Q
  27.       GO TO 10
  28. C
  29. 20    CONTINUE
  30.       PYTHAG = P
  31. C
  32.       RETURN
  33.       END
  34.