home *** CD-ROM | disk | FTP | other *** search
/ minnie.tuhs.org / unixen.tar / unixen / PDP-11 / Trees / V6 / usr / source / fort / rt1 / csqrt.f < prev    next >
Encoding:
Text File  |  1975-07-17  |  217 b   |  15 lines
  1.     function csqrt(z1)
  2.     complex csqrt, z1, z2
  3.     dimension dummy(2)
  4.     equivalence (a,z2,dummy(1)), (b,dummy(2))
  5.  
  6.     z2 = z1
  7.     c = cabs(z2)
  8.     c = sqrt(c)
  9.     b = atan2(b, a)/2.0
  10.     a = c*cos(b)
  11.     b = c*sin(b)
  12.     csqrt = z2
  13.     return
  14.     end
  15.