home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / CXSQRT.C < prev    next >
Encoding:
Text File  |  1984-08-01  |  352 b   |  17 lines
  1.     cxsqrt(a,b)
  2.  
  3.       /*subroutine finds the square root of a complex number.*/
  4.  
  5.       float a[],b[];
  6.  
  7.     {
  8.       float c,d;
  9.       extern double fabs(),sqrt();
  10.  
  11.       d = a[1];
  12.       c = sqrt(a[0]*a[0] + a[1]*a[1]);
  13.       b[1] = sqrt(fabs(.5*(-a[0] + c)));
  14.       b[0] = sqrt(fabs(.5*(a[0] + c)));
  15.       if(d < 0.0) b[1]=-b[1];
  16.     }
  17.