################################
# Diagram information
################################

name = "debox"

edges = [[1, 4], [1, 2], [2, 3], [3, 4], [3, 4]]
nodes = [1, 2, 3, 4]
internal_masses = [m2, m2, m2, m2, m2]
external_masses = [0, 0, 0, 0]

U = x[1]*x[4] + x[1]*x[5] + x[2]*x[4] + x[2]*x[5] + x[3]*x[4] + x[3]*x[5] + x[4]*x[5]
F = -m2*x[1]^2*x[4] - m2*x[1]^2*x[5] - 2*m2*x[1]*x[2]*x[4] - 2*m2*x[1]*x[2]*x[5] + (-2*m2 + s)*x[1]*x[3]*x[4] + (-2*m2 + s)*x[1]*x[3]*x[5] - m2*x[1]*x[4]^2 - 3*m2*x[1]*x[4]*x[5] - m2*x[1]*x[5]^2 - m2*x[2]^2*x[4] - m2*x[2]^2*x[5] - 2*m2*x[2]*x[3]*x[4] - 2*m2*x[2]*x[3]*x[5] - m2*x[2]*x[4]^2 + (-3*m2 + t)*x[2]*x[4]*x[5] - m2*x[2]*x[5]^2 - m2*x[3]^2*x[4] - m2*x[3]^2*x[5] - m2*x[3]*x[4]^2 - 3*m2*x[3]*x[4]*x[5] - m2*x[3]*x[5]^2 - m2*x[4]^2*x[5] - m2*x[4]*x[5]^2
parameters = [m2, s, t]
variables = [x[1], x[2], x[3], x[4], x[5]]

χ_generic = 25
f_vector = [21, 57, 64, 36, 10]

################################
# Component 1
################################

D[1] = m2
χ[1] = 3
weights[1] = [[-1, 0, 0, 0, 1], [-1, 0, 0, 1, 0], [0, 1, 1, -1, 1], [0, 1, 1, 1, -1], [0, -1, 0, 0, 1], [0, -1, 0, 1, 0], [1, 0, 1, -1, 1], [1, 0, 1, 1, -1], [0, 0, -1, 0, 1], [0, 0, -1, 1, 0], [1, 1, 0, -1, 1], [1, 1, 0, 1, -1], [1, 1, 1, -1, 0], [1, 1, 1, 0, -1], [-1, 0, 0, 0, 0], [-1, 0, 0, -1, 0], [-1, -1, 0, 0, 1], [-1, 0, 0, 0, -1], [-1, -1, 0, 1, 0], [0, 0, 1, -1, 1], [0, 1, 0, -1, 1], [0, 1, 1, -1, 0], [0, 0, 1, 1, -1], [0, 1, 0, 1, -1], [0, 1, 1, 0, -1], [0, -1, 0, 0, 0], [0, -1, 0, -1, 0], [0, -1, -1, 0, 1], [0, -1, 0, 0, -1], [0, -1, -1, 1, 0], [1, 0, 0, -1, 1], [1, 0, 1, -1, 0], [1, 0, 0, 1, -1], [1, 0, 1, 0, -1], [0, 0, -1, 0, 0], [0, 0, -1, -1, 0], [0, 0, -1, 0, -1], [1, 1, 0, -1, 0], [1, 1, 0, 0, -1], [0, 0, 0, -1, -1], [-1, 0, 0, -1, -1], [-1, -1, 0, 0, 0], [-1, -1, 0, -1, 0], [-1, -1, 0, 0, -1], [0, 0, 0, -1, 1], [0, 0, 1, -1, 0], [0, 1, 0, -1, 0], [0, 0, 0, 1, -1], [0, 0, 1, 0, -1], [0, 1, 0, 0, -1], [0, -1, -1, 0, 0], [0, -1, -1, -1, 0], [0, -1, -1, 0, -1], [1, 0, 0, -1, 0], [1, 0, 0, 0, -1], [0, 0, -1, -1, -1], [0, 0, 0, -1, 0], [0, 0, 0, 0, -1]]
computed_with[1] = ["PLD_sym"]

################################
# Component 2
################################

D[2] = m2 + 4//9*s
χ[2] = 24
weights[2] = [[-1, 0, -1, -1, -1]]
computed_with[2] = ["PLD_sym", "HyperInt"]

################################
# Component 3
################################

D[3] = m2 - 1//4*s
χ[3] = 22
weights[3] = [[-1, 0, -1, 0, 1], [-1, 0, -1, 1, 0], [-1, 0, -1, 0, 0]]
computed_with[3] = ["PLD_sym", "HyperInt"]

################################
# Component 4
################################

D[4] = m2 - 1//9*t
χ[4] = 24
weights[4] = [[0, -1, 0, -1, -1]]
computed_with[4] = ["PLD_sym", "HyperInt"]

################################
# Component 5
################################

D[5] = m2 - t
χ[5] = 22
weights[5] = [[0, -1, 0, -1, -1], [0, -1, 0, -1, -1], [0, -1, 0, -1, -1]]
computed_with[5] = ["PLD_sym", "HyperInt"]

################################
# Component 6
################################

D[6] = m2^2*s - 10//9*m2*s*t - 4//9*m2*t^2 + 1//9*s*t^2
χ[6] = 24
weights[6] = [[-1, -1, -1, -1, -1]]
computed_with[6] = ["PLD_sym"]

################################
# Component 7
################################

D[7] = s
χ[7] = 10
weights[7] = [[-1, 0, -1, 0, 1], [-1, 0, -1, 1, 0], [-1, 0, -1, 0, 0], [-1, -1, -1, 0, 1], [-1, -1, -1, 1, 0], [-1, 0, -1, -1, -1], [-1, -1, -1, 0, 0], [0, 1, 0, 0, 1], [0, 1, 0, 1, 0], [0, 1, 0, 0, 0], [0, 0, 0, 0, 1], [0, 0, 0, 1, 0]]
computed_with[7] = ["PLD_sym"]

################################
# Component 8
################################

D[8] = s + t
χ[8] = 22
weights[8] = [[-1, -1, -1, -1, -1], [0, 0, 0, 0, 0]]
computed_with[8] = ["PLD_sym", "HyperInt"]

################################
# Component 9
################################

D[9] = t
χ[9] = 13
weights[9] = [[0, -1, 0, -1, -1], [-1, -1, 0, -1, -1], [1, 0, 1, 0, 0], [0, -1, -1, -1, -1], [0, 0, 1, 0, 0], [1, 0, 0, 0, 0]]
computed_with[9] = ["PLD_sym"]