/* MENGHITUNG INVERS MATRIKS DENGAN METODE GAUSS-JORDAN */

#include "stdio.h"
#include "conio.h"
void main()
{ float p[20], a[20][20], t;
int m, i, j, k, x;
clrscr();
printf("\nMasukkan ukuran matriks : \n");
scanf("%d", &m);
printf("\nMasukkan nilai elemen matriks yang akan diinvers”);
printf("\ nsecara baris per baris\n");
/* Membaca matriks asli */
for(i=1; i<=m; i++)
{ printf("\n");
for(j=1; j<=m; j++)
{ printf("A(%d,%d)= ",i, j);
scanf("%f", &a[i][j]);
}
}
/* Mencetak Matriks asli */
printf("\nMatriks asli : ");
for(i=1; i<=m; i++)
{ printf("\n");
for(j=1; j<=m; j++)
printf(" %.2f", a[i][j]);
}
/* Proses inversi */
for(i=1; i<=m; i++)
{ p[i] = a[i][j];
a[i][j] = 1;
for(j=1; j<=m; j++)
{ a[i][j] = a[i][j]/p[i];
}
for(k=1; k<=m; k++)
{ if(k != i)
{ t = a[k][i];
a[k][i] = 0;
for(x=1; x<=m; x++)
a[k][x] = a[k][x] - a[i][x] * t;
}
}
}
/* Mencetak matriks hasil inversi*/
printf("\n\nMatriks inversi : \n");
for(i =1; i <=m; i++)
{ for(j=1; j<=m; j++)
printf(" %.1f", a[i][j]);
printf(" \n");
}
getch();
}