code to solve equations of the form Ax=b, where A is non-singular