/* invid03.c                          freeware                 [[Email Removed]]
*/

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

/* --------------------------------- MAIN ----------------------------------- */
main()
{
double pbA[5][5*TWOCOL]=
{
 1,1,  1,1,  2,1,  1,1,  5,1,
 2,1, -1,1,  4,1,  2,1,  4,1,
-2,1,  3,1, -0,1,  0,1,  3,1,
 2,1,  3,1,  5,1, -3,1,  2,1,
-2,1,  2,1,  1,1,  4,1, -1,1,
};
matrix mA={5,5*TWOCOL,&pbA[0][0]};

double pbID  [5][ 5*TWOCOL];matrix mID   = {5, 5*TWOCOL,&pbID  [0][0]};
double pbinvA[5][ 5*TWOCOL];matrix minvA = {5, 5*TWOCOL,&pbinvA[0][0]};
double pbAID [5][10*TWOCOL];matrix mAID  = {5,10*TWOCOL,&pbAID [0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
 printf("\n The matrix mA \n");
 printmF(&mA);

 printf("\n");
 invidF2(&mA,&mAID);
 printf("The inverse matrix is : \n");
 copyinv(&mAID,&minvA);
 printmF(      &minvA);

 printf("\n mA * minvA give \n");
 multmF(&mA, &minvA,&mID);
 printmF(           &mID);

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

 return 0;
}