simple_ssss.py

def f(x):
    return 1234 + 166 * x + 94 * x * x

x = [2,4,5]

f0 = 0
for j in range(3):
    prod = 1
    for m in range(3):
        if m != j:
            prod *= (x[m] / (x[m] - x[j]))
    f0 += f(x[j]) * prod
        
print(f0)