aar_gssjd.c

一些矩阵运算的C代码实现：矩阵一般运算

/*         .h         freeware           http://groups.yahoo.com/group/mathc/ */

/* --------------------------------- INCLUDES ------------------------------- */
#include "v_ahfile.h"

/* --------------------------------- MAIN ----------------------------------- */
int main(void)
{
double pb[6][7]=
{
1, 3,  3,  1, -4, -1, -3,
1, 1,  7,  1,  4, -2, -2,
2, 8,  1, -6,  4, -3, -1,
1, 1, -6,  1,  1, -1, -3,
1, 4, -1,  5,  1, -2, -2,
1, 4, -1,  1,  1, -2, -1,
};
mR m={6,7,&pb[0][0]};

/*-------------------------------- PROGRAM ---------------------------------- */

  printf(" M \n");
  p_mR(&m);

  gaussjordan_mR(&m);
            p_mR(&m);

 printf("\n Press Return to continue");
 getchar();

 return 0;
}