6bint0.c

c 算法第三版 徐士良的 源代码 非常好

  #include "stdio.h"
  #include "6bint.c"
  main()
  { int k;
    double eps=0.000001,x[4];
    static double a[4][4]=
                 {{3.4336,-0.5238,0.67105,-0.15272},
                  {-0.5238,3.28326,-0.73051,-0.2689},
                  {0.67105,-0.73051,4.02612,0.01835},
                  {-0.15272,-0.2689,0.01835,2.75702}};
    static double b[4]={-1.0,1.5,2.5,-2.0};
    if (bint(a,4,b,eps,x)!=0)
      { printf("\n");
        for (k=0; k<=3; k++)
        printf("x[%d]=%13.7e\n",k,x[k]);
        printf("\n");
      }
  }