I have two matrices, A (a 4 x 4 matrix) and b (a 4 x 1) matrix, as shown:
A = ([[2, 3, 5, 2],
[2, 4, 1, 1],
[1, 7, 4, 2],
[3, 8, 2, 3]])
b = ([[2],
[2],
[1],
[5]])
I am trying to solve for x in Ax = b. This means that x = A^-1 * b. I tried the built-in matrix solver from numpy but I want to know if there is another way to do so without using this. Perhaps a something with a more lines of code than numpy's solver?
x = np.linalg.solve(A,b)