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)