aar_pow_mr.c

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

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

/* ------------------------------- INCLUDES --------------------------------- */
#include "v_ahfile.h"
/* --------------------------------- MAIN ----------------------------------- */
int main(void)
{
double pbA[5][5]=
{
        0.0,        7.0,        0.0,       -6.0,        0.0,
        7.0,       -2.0,       -2.0,       -2.0,        3.0,
        1.0,        9.0,        1.0,        3.0,       -2.0,
        7.0,        1.0,        0.0,        5.0,        7.0,
        9.0,       -2.0,        1.0,        5.0,        9.0,
};mR mA={5,5,&pbA[0][0]};

double pbA_pn[5][5];mR mA_pn={5,5,&pbA_pn[0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
 clrscrn();
 printf(" A \n");   p_mR(&mA);   printf("\n");
 pow_mR(&mA, 0,&mA_pn);
 printf(" A**0 \n");p_mR(&mA_pn);printf("\n");
 getchar();

 clrscrn();
 printf(" A \n");   p_mR(&mA);   printf("\n");
 pow_mR(&mA, 1,&mA_pn);
 printf(" A**1 \n");p_mR(&mA_pn);printf("\n");
 getchar();

 clrscrn();
 printf(" A \n");   p_mR(&mA);   printf("\n");
 pow_mR(&mA, 2,&mA_pn);
 printf(" A**2 \n");p_mR(&mA_pn);printf("\n");
 getchar();

 clrscrn();
 printf(" A \n");   p_mR(&mA);   printf("\n");
 pow_mR(&mA, 5,&mA_pn);
 printf(" A**5 \n");p_mR(&mA_pn);printf("\n");

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

 return 0;
}