A = [1 2 3 4; 5 2 5 1; 12 3 1 4; 6 2 15 3];

C = cofactorMatrix(A);
D = determinant(A);
I = (1/D)*C';
