gausd0.c

来自「数值分析算法描述与习题解答」· C语言 代码 · 共 16 行

  #include "stdio.h"
  #include "gausdl.c"
  main()
  { int i;
    double eps;
    static double a[4][4]={ {7.0,2.0,1.0,-2.0},{9.0,15.0,3.0,-2.0},
                          {-2.0,-2.0,11.0,5.0},{1.0,3.0,2.0,13.0}};
    static double x[5],b[4]={4.0,7.0,-1.0,0.0};
    eps=0.000001;
    if (gausdl(4,a,b,x,eps)>0)
      for (i=0;i<=3; i++)
        printf("x(%d)=%13.7e\n",i,x[i]);
  }