- 01
- 02
- 03
- 04
- 05
- 06
- 07
- 08
- 09
- 10
- 11
- 12
- 13
- 14
- 15
- 16
- 17
- 18
- 19
- 20
- 21
- 22
- 23
- 24
- 25
- 26
- 27
- 28
- 29
- 30
- 31
- 32
- 33
- 34
- 35
- 36
- 37
- 38
- 39
- 40
- 41
- 42
- 43
Matrix Matrix::inverse ()
{
Matrix Temp;
Matrix A;
Matrix AL(4,4);
int c1,c2;
for (int i = 0; i < 5; i++)
for (int j = 0; j < 5; j++)
{
c1 = 0;
for (int k = 0; k < 5; k++)
{
if (k == i) continue;
c2 = 0;
for (int l = 0; l < 5; l++)
{
if (l == j) continue;
AL[c1][c2] = array[k][l];
c2++;
}
c1++;
}
A[i][j] = pow(-1,static_cast<double>(i + j + 2))*AL.determinant();
}
for (int i = 0; i < 5; i++)
for (int j = 0; j < 5; j++)
{
Temp[i][j] = 1/determinant()*A[i][j];
}
return Temp;
}