# MatveyPanchenko DSAI-03
import sympy as sp, math

def rref(M, e=1e-12):
    A, m, n, r, c, piv = [row[:] for row in M], len(M), len(M[0]), 0, 0, []
    while r < m and c < n:
        k = next((i for i in range(r, m) if abs(A[i][c]) > e), None)
        if k is None: c += 1; continue
        A[r], A[k] = A[k], A[r]
        d = A[r][c]; A[r] = [x / d for x in A[r]]
        for i in range(m):
            if i != r and abs(A[i][c]) > e:
                f = A[i][c]; A[i] = [u - f * v for u, v in zip(A[i], A[r])]
        piv.append(c); r += 1; c += 1
    return A, piv

def null_space(M):
    R, piv = rref(M); n = len(R[0]); free = [j for j in range(n) if j not in piv]
    if not free: return [[0]*n]
    B = []
    for f in free:
        v = [0]*n; v[f] = 1
        for i, p in enumerate(piv): v[p] = -R[i][f]
        B.append(v)
    return B

si = lambda A, l: [[A[i][j] - (l if i == j else 0) for j in range(3)] for i in range(3)]

def char_coeffs(A):
    a, b, c = A[0]; d, e, f = A[1]; g, h, i = A[2]
    tr = a + e + i
    s2 = a*e + a*i + e*i - b*d - c*g - f*h
    det = a*(e*i - f*h) - b*(d*i - f*g) + c*(d*h - e*g)
    return 1, -tr, s2, -det

A = [list(map(float, input().split())) for _ in range(3)]
p = char_coeffs(A)
λ = sp.symbols('λ')
poly = p[0]*λ**3 + p[1]*λ**2 + p[2]*λ + p[3]
vals = [complex(sp.N(v)) for v in sp.solve(poly, λ)]
vals = [v.real if abs(v.imag) < 1e-10 else v for v in vals]

print("Eigenvalues:")
for v in vals:
    print(v)

print("Eigenvectors:")
for l in vals:
    v = next((u for u in null_space(si(A, l)) if any(abs(x) > 1e-12 for x in u)), [0]*3)
    n = math.sqrt(sum(abs(x)**2 for x in v))
    v = [x / n for x in v] if n > 1e-12 else v
    print(v)