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