def main():
print("Start working")
A = numpy.array([[3, -1, 1],
[1, 2, 4],
[5, 1, 2]], "f")
b = numpy.array([12, 6, 3], "f")
x = gauss.num_gauss(A, b)
A = sympy.Matrix([i for i in A])
print("Triangle matrix A:")
sympy.pprint(A)
print("\nVector b: {0}".format(b))
print("\nResults: {0}".format(x))
def main():
# A = numpy.array([[3, -2, 4, 5],
# [1, -1, -1, 1],
# [2, -2, 3, -3],
# [4, 2, -6, -1]], "f")
# b = numpy.array([5, 0, 1, 0], "f")
A = numpy.array([[1, 2, 3, -2],
[2, -1, -2, -3],
[3, 2, -1, 2],
[2, -3, 2, 1]], "f")
b = numpy.array([6, 8, 4, -8], "f")
print("Initial matrix A:")
sympy.pprint(sympy.Matrix([i for i in A]))
x = gauss.num_gauss(A, b)
A = sympy.Matrix([i for i in A])
print("\nTriangle matrix A:")
sympy.pprint(A)
print("\nVector b: {0}".format(b))
print("\nResults: {0}".format(x))
