home *** CD-ROM | disk | FTP | other *** search
/ minnie.tuhs.org / unixen.tar / unixen / PDP-11 / Trees / V6 / usr / source / fort / rt1 / dcsqrt.f < prev    next >
Encoding:
Text File  |  1975-07-17  |  250 b   |  15 lines
  1.     function dcsqrt(z1)
  2.     double complex dcsqrt, z1, z2
  3.     double precision dummy(2), a, b, c
  4.     equivalence (a,z2,dummy(1)), (b,dummy(2))
  5.  
  6.     z2 = z1
  7.     c = dcabs(z2)
  8.     c = dsqrt(c)
  9.     b = datan2(b, a)/2.0d0
  10.     a = c*dcos(b)
  11.     b = c*dsin(b)
  12.     dcsqrt = z2
  13.     return
  14.     end
  15.