Gauss Siedel Iteration Method.py

# Gauss Seidel Iteration


f1 = lambda x, y, z: (7 + y) / 2  # x
f2 = lambda x, y, z: (1 + z + x) / 2  # y
f3 = lambda x, y, z: (1 + y) / 2  # z

# Initial setup
x0 = 0
y0 = 0
z0 = 0
count = 1

# Reading tolerable error
e = float(input("Enter Tolerable Error: "))

# Implementation of Gauss Seidel Iteration
print("\nIteration\t\tx \ty \tz \n")

condition = True

while condition:
    x1 = f1(x0, y0, z0)
    y1 = f2(x1, y0, z0)
    z1 = f3(x1, y1, z0)
    print("%d\t\t %0.4f\t %0.4f\t %0.4f \n" % (count, x1, y1, z1))
    e1 = abs(x0 - x1)
    e2 = abs(y0 - y1)
    e3 = abs(z0 - z1)

    count += 1
    x0 = x1
    y0 = y1
    z0 = z1
    condition = e1 > e and e2 > e and e3 > e

print("\nSolution: x=%0.3f, y=%0.3f and z = %0.3f\n" % (x1, y1, z1))