qrtt0.c

来自「qr方法求赫申伯格矩阵的全部特征值」· C语言 代码 · 共 27 行

  #include "stdio.h"
  #include "qrtt.c"
  #include "hhbg.c"
  main()
  { int i,j,jt=60;
    double eps=0.000001;
    static double u[5],v[5];
    static double a[5][5]={ {1.0,6.0,-3.0,-1.0,7.0},
       {8.0,-15.0,18.0,5.0,4.0},{-2.0,11.0,9.0,15.0,20.0},
       {-13.0,2.0,21.0,30.0,-6.0},{17.0,22.0,-5.0,3.0,6.0}};
    hhbg(5,a);
    printf("MAT H IS:\n");
    for (i=0; i<=4; i++)
      { for (j=0; j<=4; j++)
          printf("%13.7e ",a[i][j]);
        printf("\n");
      }
    printf("\n");
    i=qrtt(5,a,u,v,eps,jt);
    if (i>0)
      for (i=0; i<=4; i++)
        printf("%13.7e +J %13.7e\n",u[i],v[i]);
    printf("\n");
  }