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

name = "debox"

edges = [[1, 4], [1, 2], [2, 3], [3, 4], [3, 4]]
nodes = [1, 2, 3, 4]
internal_masses = [m[1], m[2], m[3], m[4], m[5]]
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 = -m[1]*x[1]^2*x[4] - m[1]*x[1]^2*x[5] + (-m[1] - m[2])*x[1]*x[2]*x[4] + (-m[1] - m[2])*x[1]*x[2]*x[5] + (-m[1] - m[3] + s)*x[1]*x[3]*x[4] + (-m[1] - m[3] + s)*x[1]*x[3]*x[5] - m[4]*x[1]*x[4]^2 + (-m[1] - m[4] - m[5])*x[1]*x[4]*x[5] - m[5]*x[1]*x[5]^2 - m[2]*x[2]^2*x[4] - m[2]*x[2]^2*x[5] + (-m[2] - m[3])*x[2]*x[3]*x[4] + (-m[2] - m[3])*x[2]*x[3]*x[5] - m[4]*x[2]*x[4]^2 + (-m[2] - m[4] - m[5] + t)*x[2]*x[4]*x[5] - m[5]*x[2]*x[5]^2 - m[3]*x[3]^2*x[4] - m[3]*x[3]^2*x[5] - m[4]*x[3]*x[4]^2 + (-m[3] - m[4] - m[5])*x[3]*x[4]*x[5] - m[5]*x[3]*x[5]^2 - m[4]*x[4]^2*x[5] - m[5]*x[4]*x[5]^2
parameters = [m[1], m[2], m[3], m[4], m[5], s, t]
variables = [x[1], x[2], x[3], x[4], x[5]]

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

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

D[1] = m[1]
χ[1] = 27
weights[1] = [[-1, 0, 0, 0, 1], [-1, 0, 0, 1, 0], [-1, 0, 0, 0, 0]]
computed_with[1] = ["PLD_sym"]

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

D[2] = m[1] - m[2]
χ[2] = 28
weights[2] = [[-1, -1, 0, 0, 1], [-1, -1, 0, 1, 0], [-1, -1, 0, 0, 0]]
computed_with[2] = ["PLD_sym"]

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

D[3] = m[1]*m[2] - m[1]*m[3] - m[2]^2 + m[2]*m[3] - m[2]*s
χ[3] = 28
weights[3] = [[-1, -1, -1, 0, 1], [-1, -1, -1, 1, 0], [-1, -1, -1, 0, 0]]
computed_with[3] = ["PLD_num"]

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

D[4] = m[1]^2 - 2*m[1]*m[3] - 2*m[1]*s + m[3]^2 - 2*m[3]*s + s^2
χ[4] = 28
weights[4] = [[-1, 0, -1, 0, 1], [-1, 0, -1, 1, 0], [-1, 0, -1, 0, 0]]
computed_with[4] = ["PLD_sym"]

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

D[5] = m[1]^2*m[3]^2 - 2*m[1]^2*m[3]*m[4] - 2*m[1]^2*m[3]*m[5] + m[1]^2*m[4]^2 - 2*m[1]^2*m[4]*m[5] + m[1]^2*m[5]^2 - 2*m[1]*m[3]^2*m[4] - 2*m[1]*m[3]^2*m[5] + 4*m[1]*m[3]*m[4]^2 + 8*m[1]*m[3]*m[4]*m[5] + 2*m[1]*m[3]*m[4]*s + 4*m[1]*m[3]*m[5]^2 + 2*m[1]*m[3]*m[5]*s - 2*m[1]*m[4]^3 + 2*m[1]*m[4]^2*m[5] - 2*m[1]*m[4]^2*s + 2*m[1]*m[4]*m[5]^2 + 4*m[1]*m[4]*m[5]*s - 2*m[1]*m[5]^3 - 2*m[1]*m[5]^2*s + m[3]^2*m[4]^2 - 2*m[3]^2*m[4]*m[5] + m[3]^2*m[5]^2 - 2*m[3]*m[4]^3 + 2*m[3]*m[4]^2*m[5] - 2*m[3]*m[4]^2*s + 2*m[3]*m[4]*m[5]^2 + 4*m[3]*m[4]*m[5]*s - 2*m[3]*m[5]^3 - 2*m[3]*m[5]^2*s + m[4]^4 - 4*m[4]^3*m[5] + 2*m[4]^3*s + 6*m[4]^2*m[5]^2 - 2*m[4]^2*m[5]*s + m[4]^2*s^2 - 4*m[4]*m[5]^3 - 2*m[4]*m[5]^2*s - 2*m[4]*m[5]*s^2 + m[5]^4 + 2*m[5]^3*s + m[5]^2*s^2
χ[5] = 30
weights[5] = [[-1, 0, -1, -1, -1]]
computed_with[5] = ["PLD_num"]

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

D[6] = m[1]^4 - 2*m[1]^3*m[2] - 2*m[1]^3*m[4] - 2*m[1]^3*m[5] + 2*m[1]^3*t + m[1]^2*m[2]^2 + 4*m[1]^2*m[2]*m[4] + 4*m[1]^2*m[2]*m[5] - 2*m[1]^2*m[2]*t + m[1]^2*m[4]^2 - 2*m[1]^2*m[4]*m[5] - 2*m[1]^2*m[4]*t + m[1]^2*m[5]^2 - 2*m[1]^2*m[5]*t + m[1]^2*t^2 - 2*m[1]*m[2]^2*m[4] - 2*m[1]*m[2]^2*m[5] - 2*m[1]*m[2]*m[4]^2 + 4*m[1]*m[2]*m[4]*m[5] + 2*m[1]*m[2]*m[4]*t - 2*m[1]*m[2]*m[5]^2 + 2*m[1]*m[2]*m[5]*t + m[2]^2*m[4]^2 - 2*m[2]^2*m[4]*m[5] + m[2]^2*m[5]^2
χ[6] = 30
weights[6] = [[-1, -1, 0, -1, -1]]
computed_with[6] = ["PLD_num"]

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

D[7] = m[1]^4*t^4 + 4*m[1]^3*m[2]*s*t^3 - 8*m[1]^3*m[3]*s*t^3 - 4*m[1]^3*m[3]*t^4 + 4*m[1]^3*m[4]*s*t^3 + 4*m[1]^3*m[5]*s*t^3 - 4*m[1]^3*s*t^4 + 6*m[1]^2*m[2]^2*s^2*t^2 - 16*m[1]^2*m[2]*m[3]*s^2*t^2 - 4*m[1]^2*m[2]*m[3]*s*t^3 + 4*m[1]^2*m[2]*m[4]*s^2*t^2 - 8*m[1]^2*m[2]*m[4]*s*t^3 + 4*m[1]^2*m[2]*m[5]*s^2*t^2 - 8*m[1]^2*m[2]*m[5]*s*t^3 - 12*m[1]^2*m[2]*s^2*t^3 + 16*m[1]^2*m[3]^2*s^2*t^2 + 16*m[1]^2*m[3]^2*s*t^3 + 6*m[1]^2*m[3]^2*t^4 - 16*m[1]^2*m[3]*m[4]*s^2*t^2 - 4*m[1]^2*m[3]*m[4]*s*t^3 - 16*m[1]^2*m[3]*m[5]*s^2*t^2 - 4*m[1]^2*m[3]*m[5]*s*t^3 + 16*m[1]^2*m[3]*s^2*t^3 + 4*m[1]^2*m[3]*s*t^4 + 6*m[1]^2*m[4]^2*s^2*t^2 + 4*m[1]^2*m[4]*m[5]*s^2*t^2 - 12*m[1]^2*m[4]*s^2*t^3 + 6*m[1]^2*m[5]^2*s^2*t^2 - 12*m[1]^2*m[5]*s^2*t^3 + 6*m[1]^2*s^2*t^4 + 4*m[1]*m[2]^3*s^3*t - 8*m[1]*m[2]^2*m[3]*s^3*t + 4*m[1]*m[2]^2*m[3]*s^2*t^2 - 4*m[1]*m[2]^2*m[4]*s^3*t - 16*m[1]*m[2]^2*m[4]*s^2*t^2 - 4*m[1]*m[2]^2*m[5]*s^3*t - 16*m[1]*m[2]^2*m[5]*s^2*t^2 - 12*m[1]*m[2]^2*s^3*t^2 - 16*m[1]*m[2]*m[3]^2*s^2*t^2 - 4*m[1]*m[2]*m[3]^2*s*t^3 + 16*m[1]*m[2]*m[3]*m[4]*s^3*t + 56*m[1]*m[2]*m[3]*m[4]*s^2*t^2 + 16*m[1]*m[2]*m[3]*m[4]*s*t^3 + 16*m[1]*m[2]*m[3]*m[5]*s^3*t + 56*m[1]*m[2]*m[3]*m[5]*s^2*t^2 + 16*m[1]*m[2]*m[3]*m[5]*s*t^3 + 16*m[1]*m[2]*m[3]*s^3*t^2 - 8*m[1]*m[2]*m[3]*s^2*t^3 - 4*m[1]*m[2]*m[4]^2*s^3*t - 16*m[1]*m[2]*m[4]^2*s^2*t^2 + 40*m[1]*m[2]*m[4]*m[5]*s^3*t + 32*m[1]*m[2]*m[4]*m[5]*s^2*t^2 - 8*m[1]*m[2]*m[4]*s^3*t^2 + 16*m[1]*m[2]*m[4]*s^2*t^3 - 4*m[1]*m[2]*m[5]^2*s^3*t - 16*m[1]*m[2]*m[5]^2*s^2*t^2 - 8*m[1]*m[2]*m[5]*s^3*t^2 + 16*m[1]*m[2]*m[5]*s^2*t^3 + 12*m[1]*m[2]*s^3*t^3 - 8*m[1]*m[3]^3*s*t^3 - 4*m[1]*m[3]^3*t^4 - 16*m[1]*m[3]^2*m[4]*s^2*t^2 - 4*m[1]*m[3]^2*m[4]*s*t^3 - 16*m[1]*m[3]^2*m[5]*s^2*t^2 - 4*m[1]*m[3]^2*m[5]*s*t^3 + 16*m[1]*m[3]^2*s^2*t^3 + 4*m[1]*m[3]^2*s*t^4 - 8*m[1]*m[3]*m[4]^2*s^3*t + 4*m[1]*m[3]*m[4]^2*s^2*t^2 - 48*m[1]*m[3]*m[4]*m[5]*s^3*t - 40*m[1]*m[3]*m[4]*m[5]*s^2*t^2 + 16*m[1]*m[3]*m[4]*s^3*t^2 - 8*m[1]*m[3]*m[4]*s^2*t^3 - 8*m[1]*m[3]*m[5]^2*s^3*t + 4*m[1]*m[3]*m[5]^2*s^2*t^2 + 16*m[1]*m[3]*m[5]*s^3*t^2 - 8*m[1]*m[3]*m[5]*s^2*t^3 - 8*m[1]*m[3]*s^3*t^3 + 4*m[1]*m[3]*s^2*t^4 + 4*m[1]*m[4]^3*s^3*t - 4*m[1]*m[4]^2*m[5]*s^3*t - 12*m[1]*m[4]^2*s^3*t^2 - 4*m[1]*m[4]*m[5]^2*s^3*t - 8*m[1]*m[4]*m[5]*s^3*t^2 + 12*m[1]*m[4]*s^3*t^3 + 4*m[1]*m[5]^3*s^3*t - 12*m[1]*m[5]^2*s^3*t^2 + 12*m[1]*m[5]*s^3*t^3 - 4*m[1]*s^3*t^4 + m[2]^4*s^4 + 4*m[2]^3*m[3]*s^3*t - 4*m[2]^3*m[4]*s^4 - 8*m[2]^3*m[4]*s^3*t - 4*m[2]^3*m[5]*s^4 - 8*m[2]^3*m[5]*s^3*t - 4*m[2]^3*s^4*t + 6*m[2]^2*m[3]^2*s^2*t^2 - 4*m[2]^2*m[3]*m[4]*s^3*t - 16*m[2]^2*m[3]*m[4]*s^2*t^2 - 4*m[2]^2*m[3]*m[5]*s^3*t - 16*m[2]^2*m[3]*m[5]*s^2*t^2 - 12*m[2]^2*m[3]*s^3*t^2 + 6*m[2]^2*m[4]^2*s^4 + 16*m[2]^2*m[4]^2*s^3*t + 16*m[2]^2*m[4]^2*s^2*t^2 + 4*m[2]^2*m[4]*m[5]*s^4 - 32*m[2]^2*m[4]*m[5]*s^3*t - 32*m[2]^2*m[4]*m[5]*s^2*t^2 + 4*m[2]^2*m[4]*s^4*t + 16*m[2]^2*m[4]*s^3*t^2 + 6*m[2]^2*m[5]^2*s^4 + 16*m[2]^2*m[5]^2*s^3*t + 16*m[2]^2*m[5]^2*s^2*t^2 + 4*m[2]^2*m[5]*s^4*t + 16*m[2]^2*m[5]*s^3*t^2 + 6*m[2]^2*s^4*t^2 + 4*m[2]*m[3]^3*s*t^3 + 4*m[2]*m[3]^2*m[4]*s^2*t^2 - 8*m[2]*m[3]^2*m[4]*s*t^3 + 4*m[2]*m[3]^2*m[5]*s^2*t^2 - 8*m[2]*m[3]^2*m[5]*s*t^3 - 12*m[2]*m[3]^2*s^2*t^3 - 4*m[2]*m[3]*m[4]^2*s^3*t - 16*m[2]*m[3]*m[4]^2*s^2*t^2 + 40*m[2]*m[3]*m[4]*m[5]*s^3*t + 32*m[2]*m[3]*m[4]*m[5]*s^2*t^2 - 8*m[2]*m[3]*m[4]*s^3*t^2 + 16*m[2]*m[3]*m[4]*s^2*t^3 - 4*m[2]*m[3]*m[5]^2*s^3*t - 16*m[2]*m[3]*m[5]^2*s^2*t^2 - 8*m[2]*m[3]*m[5]*s^3*t^2 + 16*m[2]*m[3]*m[5]*s^2*t^3 + 12*m[2]*m[3]*s^3*t^3 - 4*m[2]*m[4]^3*s^4 - 8*m[2]*m[4]^3*s^3*t + 4*m[2]*m[4]^2*m[5]*s^4 + 8*m[2]*m[4]^2*m[5]*s^3*t + 4*m[2]*m[4]^2*s^4*t + 16*m[2]*m[4]^2*s^3*t^2 + 4*m[2]*m[4]*m[5]^2*s^4 + 8*m[2]*m[4]*m[5]^2*s^3*t - 40*m[2]*m[4]*m[5]*s^4*t - 32*m[2]*m[4]*m[5]*s^3*t^2 + 4*m[2]*m[4]*s^4*t^2 - 8*m[2]*m[4]*s^3*t^3 - 4*m[2]*m[5]^3*s^4 - 8*m[2]*m[5]^3*s^3*t + 4*m[2]*m[5]^2*s^4*t + 16*m[2]*m[5]^2*s^3*t^2 + 4*m[2]*m[5]*s^4*t^2 - 8*m[2]*m[5]*s^3*t^3 - 4*m[2]*s^4*t^3 + m[3]^4*t^4 + 4*m[3]^3*m[4]*s*t^3 + 4*m[3]^3*m[5]*s*t^3 - 4*m[3]^3*s*t^4 + 6*m[3]^2*m[4]^2*s^2*t^2 + 4*m[3]^2*m[4]*m[5]*s^2*t^2 - 12*m[3]^2*m[4]*s^2*t^3 + 6*m[3]^2*m[5]^2*s^2*t^2 - 12*m[3]^2*m[5]*s^2*t^3 + 6*m[3]^2*s^2*t^4 + 4*m[3]*m[4]^3*s^3*t - 4*m[3]*m[4]^2*m[5]*s^3*t - 12*m[3]*m[4]^2*s^3*t^2 - 4*m[3]*m[4]*m[5]^2*s^3*t - 8*m[3]*m[4]*m[5]*s^3*t^2 + 12*m[3]*m[4]*s^3*t^3 + 4*m[3]*m[5]^3*s^3*t - 12*m[3]*m[5]^2*s^3*t^2 + 12*m[3]*m[5]*s^3*t^3 - 4*m[3]*s^3*t^4 + m[4]^4*s^4 - 4*m[4]^3*m[5]*s^4 - 4*m[4]^3*s^4*t + 6*m[4]^2*m[5]^2*s^4 + 4*m[4]^2*m[5]*s^4*t + 6*m[4]^2*s^4*t^2 - 4*m[4]*m[5]^3*s^4 + 4*m[4]*m[5]^2*s^4*t + 4*m[4]*m[5]*s^4*t^2 - 4*m[4]*s^4*t^3 + m[5]^4*s^4 - 4*m[5]^3*s^4*t + 6*m[5]^2*s^4*t^2 - 4*m[5]*s^4*t^3 + s^4*t^4
χ[7] = 30
weights[7] = [[-1, -1, -1, -1, -1]]
computed_with[7] = ["PLD_num"]

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

D[8] = m[2]
χ[8] = 28
weights[8] = [[0, -1, 0, 0, 1], [0, -1, 0, 1, 0], [0, -1, 0, 0, 0]]
computed_with[8] = ["PLD_sym"]

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

D[9] = m[2] - m[3]
χ[9] = 28
weights[9] = [[0, -1, -1, 0, 1], [0, -1, -1, 1, 0]]
computed_with[9] = ["PLD_sym"]

################################
# Component 10
################################

D[10] = m[2]^2*m[3]^2 - 2*m[2]^2*m[3]*m[4] - 2*m[2]^2*m[3]*m[5] + m[2]^2*m[4]^2 - 2*m[2]^2*m[4]*m[5] + m[2]^2*m[5]^2 - 2*m[2]*m[3]^3 + 4*m[2]*m[3]^2*m[4] + 4*m[2]*m[3]^2*m[5] - 2*m[2]*m[3]^2*t - 2*m[2]*m[3]*m[4]^2 + 4*m[2]*m[3]*m[4]*m[5] + 2*m[2]*m[3]*m[4]*t - 2*m[2]*m[3]*m[5]^2 + 2*m[2]*m[3]*m[5]*t + m[3]^4 - 2*m[3]^3*m[4] - 2*m[3]^3*m[5] + 2*m[3]^3*t + m[3]^2*m[4]^2 - 2*m[3]^2*m[4]*m[5] - 2*m[3]^2*m[4]*t + m[3]^2*m[5]^2 - 2*m[3]^2*m[5]*t + m[3]^2*t^2
χ[10] = 30
weights[10] = [[0, -1, -1, -1, -1]]
computed_with[10] = ["PLD_num"]

################################
# Component 11
################################

D[11] = m[2]^4 - 4*m[2]^3*m[4] - 4*m[2]^3*m[5] - 4*m[2]^3*t + 6*m[2]^2*m[4]^2 + 4*m[2]^2*m[4]*m[5] + 4*m[2]^2*m[4]*t + 6*m[2]^2*m[5]^2 + 4*m[2]^2*m[5]*t + 6*m[2]^2*t^2 - 4*m[2]*m[4]^3 + 4*m[2]*m[4]^2*m[5] + 4*m[2]*m[4]^2*t + 4*m[2]*m[4]*m[5]^2 - 40*m[2]*m[4]*m[5]*t + 4*m[2]*m[4]*t^2 - 4*m[2]*m[5]^3 + 4*m[2]*m[5]^2*t + 4*m[2]*m[5]*t^2 - 4*m[2]*t^3 + m[4]^4 - 4*m[4]^3*m[5] - 4*m[4]^3*t + 6*m[4]^2*m[5]^2 + 4*m[4]^2*m[5]*t + 6*m[4]^2*t^2 - 4*m[4]*m[5]^3 + 4*m[4]*m[5]^2*t + 4*m[4]*m[5]*t^2 - 4*m[4]*t^3 + m[5]^4 - 4*m[5]^3*t + 6*m[5]^2*t^2 - 4*m[5]*t^3 + t^4
χ[11] = 30
weights[11] = [[0, -1, 0, -1, -1]]
computed_with[11] = ["PLD_num"]

################################
# Component 12
################################

D[12] = m[3]
χ[12] = 27
weights[12] = [[0, 0, -1, 0, 1], [0, 0, -1, 1, 0], [0, 0, -1, 0, 0]]
computed_with[12] = ["PLD_sym"]

################################
# Component 13
################################

D[13] = m[4]
χ[13] = 20
weights[13] = [[0, 1, 1, -1, 1], [1, 0, 1, -1, 1], [1, 1, 0, -1, 1], [1, 1, 1, -1, 0], [0, 0, 1, -1, 1], [0, 1, 0, -1, 1], [0, 1, 1, -1, 0], [1, 0, 0, -1, 1], [1, 0, 1, -1, 0], [1, 1, 0, -1, 0], [0, 0, 0, -1, 1], [0, 0, 1, -1, 0], [0, 1, 0, -1, 0], [1, 0, 0, -1, 0], [0, 0, 0, -1, 0]]
computed_with[13] = ["PLD_sym", "PLD_num"]

################################
# Component 14
################################

D[14] = m[5]
χ[14] = 20
weights[14] = [[0, 1, 1, 1, -1], [1, 0, 1, 1, -1], [1, 1, 0, 1, -1], [1, 1, 1, 0, -1], [0, 0, 1, 1, -1], [0, 1, 0, 1, -1], [0, 1, 1, 0, -1], [1, 0, 0, 1, -1], [1, 0, 1, 0, -1], [1, 1, 0, 0, -1], [0, 0, 0, 1, -1], [0, 0, 1, 0, -1], [0, 1, 0, 0, -1], [1, 0, 0, 0, -1], [0, 0, 0, 0, -1]]
computed_with[14] = ["PLD_sym", "PLD_num"]

################################
# Component 15
################################

D[15] = s
χ[15] = 19
weights[15] = [[-1, 0, -1, -1, 0], [0, 1, 0, 1, 2], [-1, -1, -1, 0, 1], [-1, 0, -1, 0, -1], [0, 1, 0, 2, 1], [-1, -1, -1, 1, 0], [-1, 0, -1, -1, -1], [0, 1, 0, 1, 1], [-1, -1, -1, 0, 0], [0, 1, 0, 0, 1], [-1, -1, -1, -1, 0], [0, 0, 0, 1, 2], [0, 1, 0, 1, 0], [-1, -1, -1, 0, -1], [0, 0, 0, 2, 1], [0, 1, 0, 0, 0], [0, 0, 0, 1, 1], [0, 0, 0, 0, 1], [0, 0, 0, 1, 0]]
computed_with[15] = ["PLD_sym", "PLD_num"]

################################
# Component 16
################################

D[16] = s + t
χ[16] = 28
weights[16] = [[-1, -1, -1, -1, -1], [0, 0, 0, 0, 0]]
computed_with[16] = ["PLD_num"]

################################
# Component 17
################################

D[17] = t
χ[17] = 19
weights[17] = [[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[17] = ["PLD_num"]