jcbi0.c

数值分析算法描述与习题解答

  #include "stdio.h"
  #include "jcbi.c"
  main()
  { int i,j;
    double eps;
    static double v[3][3];
    static double a[3][3]={ {2.0,-1.0,0.0},{-1.0,2.0,-1.0},
                            {0.0,-1.0,2.0}};
    eps=0.000001;
    i=jcbi(3,a,v,eps,100);
    if (i>0)
      { for (i=0; i<=2; i++)
          printf("%13.7e  ",a[i][i]);
        printf("\n\n");
        for (i=0; i<=2; i++)
          { for (j=0; j<=2; j++)
              printf("%13.7e  ",v[i][j]);
            printf("\n");
          }
        printf("\n");
      }
  }