inarry.cpp

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

#include <stdio.h>
void inarry(float x0[],int *p)  
{
 int i,n,j,k;
 FILE *fp;
 fp=fopen("inarry.dat","r");
   fscanf(fp,"%d",&n);
   *p=n;
   printf(" N=%d \n",n);
   for (i=0;i<n;i++)
    { for(j=0;j<n;j++)
    {
    k=i*n+j;
    fscanf(fp,"%f",&x0[k]);
    printf(" X(%d)(%d)=%f    ",i,j,x0[i]);
    
    }
     printf("\n");
    }
  fclose(fp);
  printf("\n");
}