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

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

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

double pbID  [3][3*TWOCOL];  matrix mID   = {3,3*TWOCOL,&pbID  [0][0]};
double pbinvA[3][3*TWOCOL];  matrix minvA = {3,3*TWOCOL,&pbinvA[0][0]};
double pbAID [3][6*TWOCOL];  matrix mAID  = {3,6*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;
}