main.cpp

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

#include <stdio.h>
#include <malloc.h>
#define num 100
void inarry(float x[], int *n);
void findr(float x[],int n,float r[]);
void outx(float x[],int n);
void outf(float r[],int n);
void main()
{
 repeat:  float x0[num];
   int n, m;
   inarry(x0, &n);
   float *r0;
   r0=(float*)calloc(2*n,4);
  
   printf("是否用QR法求全部特征值？输入1，是；输入0，不是 \n");
      scanf("%d",&m);
   if(m==1)
   findr(x0,n,r0);
   outx(r0,n);
   outf(r0,n);
   printf("是否继续输入1，是；输入0，不是 \n");
      scanf("%d",&m);
   if(m==1)
   goto repeat;
   else 
   printf("over");
  }