[THOUGHT] I will use Sympy to solve the given system for natural numbers a, b, c, then compute a^4 + b^4 + c^4 from the solutions. I'll set up the equations exactly as given and let Sympy find integer solutions, then verify they are natural numbers, and finally compute the required sum of fourth powers. [/THOUGHT] [PYTHON] from sympy import symbols, Eq, solve

a, b, c = symbols('a b c', integer=True)

solutions = solve([ Eq(a + b + c, 151), Eq(a2 + b2 + c2, 10939), Eq(a3 + b3 + c3, 957871) ], (a, b, c), domain=None)

print("Solutions:", solutions)

Compute a^4 + b^4 + c^4 for each solution

results = [] for sol in solutions: av, bv, cv = sol val = av4 + bv4 + cv**4 results.append((sol, val))

print("Fourth power sums:", results) [/PYTHON]

[THOUGHT] The system has solutions that are permutations of (21, 33, 97). For any of these, the value of a^4 + b^4 + c^4 is the same, namely 89909683 from the Python output. [/THOUGHT] \boxed{89909683}