def main():
A = sympy.Matrix([[3, -1],
[-6, 2]])
b = sympy.Matrix([0, 0])
print("\nInitial matrix:")
sympy.pprint(A)
x = gauss.sym_gauss(A, b)
sympy.pprint(x)
def main():
a = sympy.Symbol("a")
A = sympy.Matrix([[3, 4, a],
[1, -5, 2],
[4, a, 1]])
print("\nInitial matrix:")
sympy.pprint(A)
aValue = sympy.solve(sympy.Eq(A.det(), 0), a)
print("\na value: {0}".format(aValue))
for i in aValue:
Asubs = A.subs(a, i)
b = sympy.zeros(1, A.rows)
x = gauss.sym_gauss(Asubs, b)
print("\na = {0}".format(i), )
sympy.pprint(x)
