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