Hey StudentMCCS.

Here is some output from Octave for your matrix:

>> A

A =

2 5

-4 10

>> inv(A)

ans =

0.250000 -0.125000

0.100000 0.050000

>> det(A)

ans = 40

>> rref(A)

ans =

1 0

0 1

All of this (although you only need one of them) says that the matrix is of full rank so you will never ever get any inconsistent situations for solving x = A^(-1)*b since a unique solution exists for all b in R^2.

The only time you get a chance of inconsistent solutions is when your matrix is not full-rank (i.e. has linearly dependent vectors in the rows of the matrix) and if your matrix is square, the easiest way to check this computationally is to check det(A). If you want to find out what rows are linearly independent, do a row reduction and see what rows go to zero and what rows don't and the number of non-zero rows in rref(A) give the number of linearly independent rows and you can get your rank and your nullity (nullity(A) = dim(space) - rank(A)).

When you have nullity, then and only then do have a chance of inconsistency and this can be determined by doing the normal reduction with a variable matrix and seeing when you get the 0 = 1 situation and those situations are the cases where you get inconsistencies).