home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / eispack-1.0-src.tgz / tar.out / contrib / eispack / cdiv.f < prev    next >
Text File  |  1996-09-28  |  390b  |  17 lines
  1.       subroutine cdiv(ar,ai,br,bi,cr,ci)
  2.       double precision ar,ai,br,bi,cr,ci
  3. c
  4. c     complex division, (cr,ci) = (ar,ai)/(br,bi)
  5. c
  6.       double precision s,ars,ais,brs,bis
  7.       s = dabs(br) + dabs(bi)
  8.       ars = ar/s
  9.       ais = ai/s
  10.       brs = br/s
  11.       bis = bi/s
  12.       s = brs**2 + bis**2
  13.       cr = (ars*brs + ais*bis)/s
  14.       ci = (ais*brs - ars*bis)/s
  15.       return
  16.       end
  17.