6ggje0.c

来自「C语言描述的常用数值算法程序集」· C语言 代码 · 共 20 行

  #include "stdio.h"
  #include "6ggje.c"
  main()
  { int i;
    static double a[8][8]={ {0.0,0.0,-1.0,0.0,0.0,0.0,2.0,0.0},
                         {0.0,6.0,0.0,0.0,0.0,-6.0,0.0,0.0},
                         {0.0,0.0,0.0,2.0,0.0,0.0,0.0,-4.0},
                         {3.0,0.0,0.0,0.0,-2.0,0.0,1.0,0.0},
                         {0.0,0.0,6.0,0.0,0.0,0.0,5.0,0.0},
                         {1.0,0.0,0.0,0.0,-3.0,0.0,0.0,2.0},
                         {0.0,4.0,0.0,-1.0,0.0,0.0,0.0,0.0},
                         {0.0,0.0,1.0,0.0,-1.0,0.0,0.0,-2.0}};
    static double b[8]={4.0,6.0,-8.0,-2.0,27.0,-9.0,2.0,-4.0};
    if (ggje(a,8,b)!=0)
      for (i=0; i<=7; i++)
        printf("x(%d)=%13.7e\n",i,b[i]);
  }