zha.txt

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

#include <stdio.h>
#include <malloc.h>
void jsqr(float a[],float b[],float c[],int n);
void mulqr(float a[],float b[],float c[],int n);
void pd(float a[],float b[],int n);
void findr(float x0[],int n,float r0[]) 
{    
   float *r,*q,*a;
   r=(float*)calloc(n*n,4);
   q=(float*)calloc(n*n,4);
   a=(float*)calloc(n*n,4);
   jsqr(x0,q,r,n);
   mulqr(q,r,a,n);
   pd(a,r0,n);
   return ;
   free(r);
   free(q);
   free(a);