outx.cpp

用QR法求矩阵全部特征值,qr算法是幕法的推广和变形

#include <Stdio.h>
void outx(float B[],int N)
{
  int i;
  FILE *fp;
  fp=fopen("Outx.dat","w");
  for (i=0;i<N*2;i++)
   {
    fprintf(fp,"x(%d)=%f \n",i,B[i]); 
    printf("r(%d)=%f ",i/2,B[i]); 
    i++;
    printf("  +(%f) j\n",B[i]);
   }
  fclose(fp);
}