Linear Algebra: How to solve a system of equations algebraically using SymPy and Python
Import the Sympy library.
Initiate printing format
Define symbols
code:
from sympy import *
init_printing()
x,y,z = symbols("x y z")
eqs = [3*x +2*y -z -4,2*x -y +2*z -10, x -3*y -4*z -5]
A,b = linear_eq_to_matrix(eqs,[x,y,z])
A
b
sol = A.solve(b)
sol
soln = A.LUsolve(b)
soln
ans = linsolve((A,b),[x,y,z])
ans
M = A.row_join(b)
M
M.rref()
AM = A.row_join(eye(3))
AM
AM.rref()
A.inv()
A.inv()*b
b.equals(A*sol)
