aar_gram.c

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

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

/* --------------------------------- INCLUDES ------------------------------- */
#include "v_ahfile.h"
/* --------------------------------- MAIN ----------------------------------- */
int main(void)
{
   int n;
   int m;
   int i=0;
double f;

double pbu[5][5]=
{
     3,          7,          3,          3,          3,
     1,          1,          3,          3,          1,
    -2,          1,         -2,         -2,         -2,
    -2,          3,          3,          1,          1,
     3,          3,         -2,          1,          3,
};
mR mu={5,5,&pbu[0][0]};

double pbv  [5][5];mR mv  ={5,5,&pbv [0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
   clrscrn();

   printf(" Consider the vector space R**5 with the Euclidean inner product.\n");
   printf(" A Matrix U with linear independent rows vectors, a basis of R**5.\n");
   p_mR(&mu);

   printf("\n");
   p_innerprod_mR(&mu,'u');
   getchar();

   clrscrn();
   printf("\n If we apply the Gram-SchidmRt process on the vectors\n");
   printf(" on the matrixU, we obtain a orthogonal basis.\n");
   Gram_Schimdt_mR(&mu,&mv);
   p_mR(&mv);

   printf("\n");
   p_innerprod_mR(&mv,'v');

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

 return 0;
}