5strq0.c

用C语言编写的关于矩阵特征值和特征向量的计算

  #include "stdio.h"
  #include "5strq.c"
  main()
  { int i,j;
    static double b[5],c[5],q[5][5];
    static double a[5][5]={ {10.0,1.0,2.0,3.0,4.0},
        {1.0,9.0,-1.0,2.0,-3.0},{2.0,-1.0,7.0,3.0,-5.0},
        {3.0,2.0,3.0,12.0,-1.0},{4.0,-3.0,-5.0,-1.0,15.0}};
    strq(a,5,q,b,c);
    printf("MAT A IS:\n");
    for (i=0; i<=4; i++)
      { for (j=0; j<=4; j++)
          printf("%13.7e ",a[i][j]);
        printf("\n");
      }
    printf("\n");
    printf("MAT Q IS:\n");
    for (i=0; i<=4; i++)
      { for (j=0; j<=4; j++)
          printf("%13.7e ",q[i][j]);
        printf("\n");
      }
    printf("\n");
    printf("MAT B IS:\n");
    for (i=0; i<=4; i++)
      printf("%13.7e ",b[i]);
    printf("\n\n");
    printf("MAT C IS:\n");
    for (i=0; i<=4; i++)
      printf("%13.7e ",c[i]);
    printf("\n\n");
  }