InvertMatrix code does not work I added DO row = 1, n DO col = 1, n res = 0.0 DO k = 1, n res = res + Matrix(row, k) * InvMatrix(k, col) ENDDO print *, "Res ", res ENDDO ENDDO