G0 = (0, 4)
x0: x(G0)
Z0: (x0, 0)
G3 = (3, 3)
x3: x(G3)
U3 = Point(Line[(x3, 0), (x3, 1)])
G2 = (2, 3)
x2: x(G2)
Z2: (x2, 0)
G1 = (1, 2)
x1: x(G1)
Z1: (x1, 0)
U2 = Point(Line[(x2, 0), (x2, 1)])
Z3: (x3, 0)
U1 = Point(Line[(x1, 0), (x1, 1)])
U0 = Point(Line[(x0, 0), (x0, 1)])
ρ_1 = Polynomial({Z0, U1})
ρ_2 = Polynomial({Z0, Z1, U2})
ρ_3 = Polynomial({Z0, Z1, Z2, U3})
ρ_0(x) = y(U0)
G4 = (4, 2)
G5 = (5, 5)
x4: x(G4)
x5: x(G5)
Z4: (x4, 0)
Z5: (x5, 0)
U4 = Point(Line[(x4, 0), (x4, 1)])
U5 = Point(Line[(x5, 0), (x5, 1)])
ρ_4 = Polynomial({Z0, Z1, Z2, Z3, U4})
ρ_5 = Polynomial({Z0, Z1, Z2, Z3, Z4, U5})
φ(x) = ρ_0(x) + ρ_1(x) + ρ_2(x) + ρ_3(x) + ρ_4(x) + ρ_5(x)