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

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

/* --------------------------------- MAIN ----------------------------------- */
 main()
{

double pb[4][5*TWOCOL]=
{
 11,1,  14,1, -16,1,  19,1,  11,1,
 12,1,  13,1,  13,1,  11,1, -14,1,
  0,1,   0,1,  17,1,  11,1,  12,1,
 12,1,  18,1,   1,1, -16,1,  14,1,
};
matrix m={4,5*TWOCOL,&pb[0][0]};

double pbT   [4][5*TWOCOL]; matrix mT    ={4,5*TWOCOL,&pbT   [0][0]};
double pbID  [4][4*TWOCOL]; matrix mID   ={4,4*TWOCOL,&pbID  [0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
 copymF(&m,&mT);

 printf("\n");
 printf(" The matrix  mA \n");
 printmF(&m);
  getchar();

 printf(" Find mID in one step \n");
 idF4(&mID,&m);
 printmF(&mID);
  getchar();

 printf("  mID * mA  \n");
 multmF(&mID,&mT,&m);
 printmF(&m);

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

 return 0;
}