#include int main(int argc, char *argv[]) { void MatrixInvert(float*, int, float*); float A[][3] = { {2, 5, 9}, {7, 3, 7}, {8, 3, 1} }; float A1[3][3]; float B[][4] = { {3, 2, 5, 9}, {1, 7, 3, 7}, {0, -9, 1, 8}, {9, 8, 3, 1} }; float B1[4][4]; int i, j; cout << "A:\n\n"; for (i = 0; i < 3; i = i + 1) { for (j = 0; j < 3; j = j + 1) cout << A[i][j] << "\t"; cout << "\n"; } MatrixInvert((float *)A, 3, (float *)A1); cout << "\nA inverse:\n\n"; for (i = 0; i < 3; i = i + 1) { for (j = 0; j < 3; j = j + 1) cout << A1[i][j] << "\t"; cout << "\n"; } cout << "\n\nB:\n\n"; for (i = 0; i < 4; i = i + 1) { for (j = 0; j < 4; j = j + 1) cout << B[i][j] << "\t"; cout << "\n"; } MatrixInvert((float *)B, 4, (float *)B1); cout << "\nB inverse:\n\n"; for (i = 0; i < 4; i = i + 1) { for (j = 0; j < 4; j = j + 1) cout << B1[i][j] << "\t"; cout << "\n"; } }