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

name = "pltrb"

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

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

χ_generic = 81
f_vector = [33, 108, 152, 115, 49, 11]

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

D[1] = M[3]
χ[1] = 64
weights[1] = [[0, -1, -1, 0, 1, 1], [1, -1, -1, 0, 0, 1], [1, -1, -1, 0, 1, 0], [-1, 0, -1, 0, 0, -1], [0, -1, -1, 0, 0, 1], [0, -1, -1, 0, 1, 0], [1, -1, -1, 0, 0, 0], [-1, -1, -1, 0, 0, -1], [0, 0, 0, 1, 1, 1], [-1, -1, -1, 0, -1, 0], [0, 1, 0, 1, 1, 0], [0, -1, -1, 0, 0, 0], [1, 0, 0, 1, 0, 1], [1, 0, 0, 1, 1, 0], [0, 0, 0, 1, 1, 0], [0, 0, 0, 1, 0, 1]]
computed_with[1] = ["PLD_sym", "PLD_num"]

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

D[2] = M[3] - 4*m2
χ[2] = 74
weights[2] = [[0, -1, -1, 0, 1, 1], [1, -1, -1, 0, 0, 1], [1, -1, -1, 0, 1, 0], [0, -1, -1, 0, 0, 1], [0, -1, -1, 0, 1, 0], [1, -1, -1, 0, 0, 0], [0, -1, -1, 0, 0, 0]]
computed_with[2] = ["PLD_sym", "PLD_num", "HyperInt"]

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

D[3] = M[3] - 9*m2
χ[3] = 80
weights[3] = [[-1, 0, -1, 0, 0, -1]]
computed_with[3] = ["PLD_sym", "HyperInt"]

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

D[4] = M[3] - m2
χ[4] = 78
weights[4] = [[-1, 0, -1, 0, 0, -1], [-1, 0, -1, 0, 0, -1], [-1, 0, -1, 0, 0, -1]]
computed_with[4] = ["PLD_sym", "HyperInt"]

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

D[5] = M[3]*M[4]*m2 - M[3]*M[4]*s - M[3]*m2^2 + M[3]*m2*s - M[4]^2*m2 + 2*M[4]*m2*s - m2*s^2
χ[5] = 79
weights[5] = []
computed_with[5] = ["HyperInt"]

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

D[6] = M[3]^2 - 2*M[3]*M[4] - 2*M[3]*s + M[4]^2 - 2*M[4]*s + s^2
χ[6] = 57
weights[6] = [[-1, -1, -1, -1, 0, 0], [-1, 0, -1, -1, 0, -1], [-1, 0, -1, 0, -1, -1], [1, 0, 0, 0, 2, 2], [0, -1, -1, 0, -1, -1], [0, -1, -1, -1, -1, 0], [2, 0, 0, 0, 1, 2], [0, -1, -1, -1, 0, -1], [2, 0, 0, 0, 2, 1], [-1, -1, -1, -1, 0, -1], [0, 0, 0, 0, 1, 1], [-1, -1, -1, -1, -1, 0], [-1, -1, -1, 0, -1, -1], [0, 1, 0, 0, 1, 0], [-1, 0, -1, -1, -1, -1], [0, 1, 0, 1, 0, 0], [1, 0, 0, 0, 1, 2], [1, 0, 0, 0, 2, 1], [1, 0, 0, 1, 0, 0], [0, -1, -1, -1, -1, -1], [2, 0, 0, 0, 1, 1], [1, 0, 0, 0, 0, 1], [1, 0, 0, 0, 1, 0], [-1, -1, -1, -1, -1, -1], [0, 0, 0, 0, 1, 0], [0, 0, 0, 0, 0, 1], [0, 0, 0, 1, 0, 0], [0, 1, 0, 0, 0, 0], [1, 0, 0, 0, 1, 1], [1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0]]
computed_with[6] = ["PLD_sym", "PLD_num", "HyperInt"]

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

D[7] = M[3]^2*m2 - 2*M[3]*M[4]*m2 + M[3]*M[4]*s - 2*M[3]*m2*s + M[4]^2*m2 - 2*M[4]*m2*s + m2*s^2
χ[7] = 74
weights[7] = [[0, -1, -1, -1, 1, 1], [1, -1, -1, -1, 0, 1], [1, -1, -1, -1, 1, 0], [0, -1, -1, -1, 0, 1], [0, -1, -1, -1, 1, 0], [1, -1, -1, -1, 0, 0], [0, -1, -1, -1, 0, 0]]
computed_with[7] = ["PLD_sym", "PLD_num", "HyperInt"]

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

D[8] = M[3]^2*m2 - 2*M[3]*M[4]*m2 + M[3]*M[4]*s - 5*M[3]*m2*s + M[4]^2*m2 - 5*M[4]*m2*s + 9*m2^2*s + 4*m2*s^2
χ[8] = 80
weights[8] = [[-1, 0, -1, 0, -1, -1]]
computed_with[8] = ["PLD_num", "HyperInt"]

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

D[9] = M[3]^2*m2 - 2*M[3]*M[4]*m2 + M[3]*M[4]*s - M[3]*m2*s + M[4]^2*m2 - M[4]*m2*s + m2^2*s
χ[9] = 79
weights[9] = []
computed_with[9] = ["HyperInt"]

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

D[10] = M[3]^2*m2 - 5*M[3]*M[4]*m2 + M[3]*M[4]*s - 2*M[3]*m2*s + 4*M[4]^2*m2 + 9*M[4]*m2^2 - 5*M[4]*m2*s + m2*s^2
χ[10] = 80
weights[10] = [[-1, 0, -1, -1, 0, -1]]
computed_with[10] = ["PLD_num", "HyperInt"]

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

D[11] = M[3]^2*m2 - 5//4*M[3]*M[4]*m2 + 1//4*M[3]*M[4]*s + 9//4*M[3]*m2^2 - 5//4*M[3]*m2*s + 1//4*M[4]^2*m2 - 1//2*M[4]*m2*s + 1//4*m2*s^2
χ[11] = 80
weights[11] = [[0, -1, -1, 0, -1, -1]]
computed_with[11] = ["PLD_num", "HyperInt"]

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

D[12] = M[3]^2*m2 - M[3]*M[4]*m2 + M[3]*M[4]*s - 2*M[3]*m2*s + M[4]*m2^2 - M[4]*m2*s + m2*s^2
χ[12] = 79
weights[12] = []
computed_with[12] = ["HyperInt"]

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

D[13] = M[3]^4*m2^2 - 2*M[3]^4*m2*s + M[3]^4*s^2 - 4*M[3]^3*M[4]*m2^2 + 5*M[3]^3*M[4]*m2*s - M[3]^3*M[4]*s^2 + 2*M[3]^3*m2^2*s - M[3]^3*m2*s^2 - M[3]^3*s^3 + 6*M[3]^2*M[4]^2*m2^2 - 4*M[3]^2*M[4]^2*m2*s + M[3]^2*M[4]^2*s^2 - 8*M[3]^2*M[4]*m2^2*s + 3*M[3]^2*M[4]*m2*s^2 - M[3]^2*M[4]*s^3 + 3*M[3]^2*m2^2*s^2 - M[3]^2*m2*s^3 + M[3]^2*s^4 - 4*M[3]*M[4]^3*m2^2 + M[3]*M[4]^3*m2*s + 10*M[3]*M[4]^2*m2^2*s - 4*M[3]*M[4]^2*m2*s^2 - 8*M[3]*M[4]*m2^2*s^2 + 5*M[3]*M[4]*m2*s^3 + 2*M[3]*m2^2*s^3 - 2*M[3]*m2*s^4 + M[4]^4*m2^2 - 4*M[4]^3*m2^2*s + 6*M[4]^2*m2^2*s^2 - 4*M[4]*m2^2*s^3 + m2^2*s^4
χ[13] = 80
weights[13] = [[-1, -1, -1, 0, -1, -1]]
computed_with[13] = ["PLD_num", "HyperInt"]

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

D[14] = M[3]^4*m2^2 - 4*M[3]^3*M[4]*m2^2 + M[3]^3*M[4]*m2*s - 4*M[3]^3*m2^2*s + 6*M[3]^2*M[4]^2*m2^2 - 4*M[3]^2*M[4]^2*m2*s + M[3]^2*M[4]^2*s^2 + 10*M[3]^2*M[4]*m2^2*s - 4*M[3]^2*M[4]*m2*s^2 + 6*M[3]^2*m2^2*s^2 - 4*M[3]*M[4]^3*m2^2 + 5*M[3]*M[4]^3*m2*s - M[3]*M[4]^3*s^2 - 8*M[3]*M[4]^2*m2^2*s + 3*M[3]*M[4]^2*m2*s^2 - M[3]*M[4]^2*s^3 - 8*M[3]*M[4]*m2^2*s^2 + 5*M[3]*M[4]*m2*s^3 - 4*M[3]*m2^2*s^3 + M[4]^4*m2^2 - 2*M[4]^4*m2*s + M[4]^4*s^2 + 2*M[4]^3*m2^2*s - M[4]^3*m2*s^2 - M[4]^3*s^3 + 3*M[4]^2*m2^2*s^2 - M[4]^2*m2*s^3 + M[4]^2*s^4 + 2*M[4]*m2^2*s^3 - 2*M[4]*m2*s^4 + m2^2*s^4
χ[14] = 80
weights[14] = [[-1, 0, -1, -1, -1, -1]]
computed_with[14] = ["PLD_num", "HyperInt"]

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

D[15] = M[4]
χ[15] = 64
weights[15] = [[0, 0, -1, -1, 1, 1], [1, 0, -1, -1, 0, 1], [1, 0, -1, -1, 1, 0], [0, 0, -1, -1, 0, 1], [0, 0, -1, -1, 1, 0], [0, 0, -1, 0, -1, -1], [1, 0, -1, -1, 0, 0], [0, 1, 0, 0, 1, 1], [-1, 0, -1, -1, -1, 0], [0, 0, -1, -1, 0, 0], [1, 1, 0, 1, 0, 0], [0, 0, -1, -1, -1, -1], [1, 1, 0, 0, 0, 1], [1, 1, 0, 0, 1, 0], [0, 1, 0, 0, 0, 1], [1, 1, 0, 0, 0, 0]]
computed_with[15] = ["PLD_sym", "PLD_num"]

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

D[16] = M[4] - 4*m2
χ[16] = 74
weights[16] = [[0, 0, -1, -1, 1, 1], [1, 0, -1, -1, 0, 1], [1, 0, -1, -1, 1, 0], [0, 0, -1, -1, 0, 1], [0, 0, -1, -1, 1, 0], [1, 0, -1, -1, 0, 0], [0, 0, -1, -1, 0, 0]]
computed_with[16] = ["PLD_sym", "PLD_num", "HyperInt"]

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

D[17] = M[4] - 9*m2
χ[17] = 80
weights[17] = [[0, 0, -1, 0, -1, -1]]
computed_with[17] = ["PLD_sym", "HyperInt"]

################################
# Component 18
################################

D[18] = M[4] - m2
χ[18] = 78
weights[18] = [[0, 0, -1, 0, -1, -1], [0, 0, -1, 0, -1, -1], [0, 0, -1, 0, -1, -1]]
computed_with[18] = ["PLD_sym", "HyperInt"]

################################
# Component 19
################################

D[19] = m2
χ[19] = 16
weights[19] = [[-1, 0, 1, 1, 0, 1], [-1, 1, 0, 1, 0, 1], [-1, 1, 1, 0, 0, 1], [-1, 1, 1, 1, 0, 0], [0, -1, 0, 0, 1, 1], [0, 0, -1, 0, 1, 1], [0, 0, 0, -1, 1, 1], [0, 1, 1, 1, 1, -1], [1, -1, 0, 0, 0, 1], [1, -1, 0, 0, 1, 0], [0, 0, 1, 1, -1, 1], [1, 0, 1, 1, 1, -1], [1, 0, -1, 0, 0, 1], [1, 0, -1, 0, 1, 0], [0, 1, 0, 1, -1, 1], [1, 1, 0, 1, 1, -1], [1, 0, 0, -1, 0, 1], [1, 0, 0, -1, 1, 0], [0, 1, 1, 0, -1, 1], [1, 1, 1, 0, 1, -1], [0, 1, 1, 1, -1, 0], [1, 1, 1, 1, 0, -1], [-1, 0, 0, 1, 0, 1], [-1, 0, 1, 0, 0, 1], [-1, 0, 1, 1, 0, 0], [-1, -1, 0, 0, 0, 0], [-1, 1, 0, 0, 0, 1], [-1, 1, 0, 1, 0, 0], [-1, 0, -1, 0, 0, 0], [-1, 1, 1, 0, 0, 0], [-1, 0, 0, -1, 0, 0], [-1, 0, 0, 0, 0, -1], [0, -1, 0, 0, 0, 1], [0, -1, 0, 0, 1, 0], [0, 0, -1, 0, 0, 1], [0, 0, -1, 0, 1, 0], [0, 0, 0, -1, 0, 1], [0, 0, 0, -1, 1, 0], [0, 0, 1, 1, 1, -1], [0, 1, 0, 1, 1, -1], [0, 1, 1, 0, 1, -1], [0, 1, 1, 1, 0, -1], [1, -1, 0, 0, 0, 0], [0, -1, 0, 0, -1, 0], [0, -1, 0, 0, 0, -1], [0, 0, 0, 1, -1, 1], [0, 0, 1, 0, -1, 1], [0, 0, 1, 1, -1, 0], [1, 0, 0, 1, 1, -1], [1, 0, 1, 0, 1, -1], [1, 0, 1, 1, 0, -1], [1, 0, -1, 0, 0, 0], [0, 0, -1, 0, -1, 0], [0, 0, -1, 0, 0, -1], [0, 1, 0, 0, -1, 1], [0, 1, 0, 1, -1, 0], [1, 1, 0, 0, 1, -1], [1, 1, 0, 1, 0, -1], [1, 0, 0, -1, 0, 0], [0, 0, 0, -1, -1, 0], [0, 0, 0, -1, 0, -1], [0, 1, 1, 0, -1, 0], [1, 1, 1, 0, 0, -1], [0, 0, 0, 0, -1, -1], [-1, 0, 0, 0, 0, 1], [-1, 0, 0, 1, 0, 0], [-1, 0, 1, 0, 0, 0], [-1, -1, 0, 0, 0, -1], [-1, 1, 0, 0, 0, 0], [0, -1, 0, 0, 0, 0], [0, 0, -1, 0, 0, 0], [0, 0, 0, -1, 0, 0], [0, 0, 0, 1, 1, -1], [0, 0, 1, 0, 1, -1], [0, 0, 1, 1, 0, -1], [0, 1, 0, 0, 1, -1], [0, 1, 0, 1, 0, -1], [0, 1, 1, 0, 0, -1], [0, 0, 0, 0, -1, 1], [0, 0, 0, 1, -1, 0], [0, 0, 1, 0, -1, 0], [1, 0, 0, 0, 1, -1], [1, 0, 0, 1, 0, -1], [1, 0, 1, 0, 0, -1], [0, 1, 0, 0, -1, 0], [1, 1, 0, 0, 0, -1], [0, 0, 0, -1, -1, -1], [-1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 1, -1], [0, 0, 0, 1, 0, -1], [0, 0, 1, 0, 0, -1], [0, 1, 0, 0, 0, -1], [0, 0, 0, 0, -1, 0], [1, 0, 0, 0, 0, -1], [0, 0, 0, 0, 0, -1]]
computed_with[19] = ["PLD_sym", "PLD_num"]

################################
# Component 20
################################

D[20] = m2 - 1//3*s
χ[20] = 79
weights[20] = []
computed_with[20] = ["HyperInt"]

################################
# Component 21
################################

D[21] = m2 - 1//4*s
χ[21] = 62
weights[21] = [[-1, 0, 1, 1, -1, 1], [-1, 1, 0, 1, -1, 1], [-1, 1, 1, 0, -1, 1], [-1, 1, 1, 1, -1, 0], [0, -1, 0, -1, 1, 1], [1, -1, 0, -1, 0, 1], [1, -1, 0, -1, 1, 0], [-1, 0, 0, 1, -1, 1], [-1, 0, 1, 0, -1, 1], [-1, 0, 1, 1, -1, 0], [-1, 1, 0, 0, -1, 1], [-1, 1, 0, 1, -1, 0], [-1, 1, 1, 0, -1, 0], [0, -1, 0, -1, 0, 1], [0, -1, 0, -1, 1, 0], [1, -1, 0, -1, 0, 0], [-1, 0, 0, 0, -1, 1], [-1, 0, 0, 1, -1, 0], [-1, 0, 1, 0, -1, 0], [-1, -1, 0, -1, -1, 0], [-1, 1, 0, 0, -1, 0], [0, -1, 0, -1, 0, 0], [-1, 0, 0, 0, -1, 0]]
computed_with[21] = ["PLD_sym", "PLD_num", "HyperInt"]

################################
# Component 22
################################

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

################################
# Component 23
################################

D[23] = m2 - s
χ[23] = 75
weights[23] = [[-1, 0, 0, -1, 0, -1], [-1, 0, 0, -1, 0, -1], [-1, 0, 0, -1, 0, -1], [0, -1, 0, 0, -1, -1], [0, -1, 0, 0, -1, -1], [0, -1, 0, 0, -1, -1]]
computed_with[23] = ["PLD_sym", "HyperInt"]

################################
# Component 24
################################

D[24] = s
χ[24] = 39
weights[24] = [[-1, 0, 1, 1, -1, 1], [-1, 1, 0, 1, -1, 1], [-1, 1, 1, 0, -1, 1], [-1, 1, 1, 1, -1, 0], [0, -1, 0, -1, 1, 1], [1, -1, 0, -1, 0, 1], [1, -1, 0, -1, 1, 0], [-1, 0, 0, 1, -1, 1], [-1, 0, 1, 0, -1, 1], [-1, 0, 1, 1, -1, 0], [-1, 1, 0, 0, -1, 1], [-1, 1, 0, 1, -1, 0], [-1, 0, 0, -1, 0, -1], [-1, 1, 1, 0, -1, 0], [0, -1, 0, -1, 0, 1], [0, -1, 0, -1, 1, 0], [0, -1, 0, 0, -1, -1], [1, -1, 0, -1, 0, 0], [-1, 0, 0, 0, -1, 1], [-1, 0, 0, 1, -1, 0], [-1, -1, -1, 0, -1, 0], [-1, -1, 0, -1, 0, -1], [-1, 0, 1, 0, -1, 0], [0, 0, 1, 0, 1, 1], [-1, -1, 0, -1, -1, 0], [-1, -1, 0, 0, -1, -1], [0, 0, 1, 1, 0, 1], [-1, 1, 0, 0, -1, 0], [-1, 0, -1, -1, -1, 0], [0, 1, 0, 1, 0, 1], [0, 1, 1, 0, 1, 0], [-1, 0, 0, -1, -1, -1], [0, 1, 1, 0, 0, 1], [0, 1, 1, 1, 0, 0], [0, -1, 0, -1, 0, 0], [1, 0, 1, 1, 0, 0], [0, -1, 0, -1, -1, -1], [1, 0, 1, 0, 0, 1], [1, 0, 1, 0, 1, 0], [-1, 0, 0, 0, -1, 0], [-1, -1, -1, -1, -1, 0], [0, 0, 0, 1, 0, 1], [0, 0, 1, 0, 1, 0], [-1, -1, 0, -1, -1, -1], [0, 0, 1, 0, 0, 1], [0, 0, 1, 1, 0, 0], [0, 1, 0, 0, 0, 1], [0, 1, 1, 0, 0, 0], [1, 0, 1, 0, 0, 0], [0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0]]
computed_with[24] = ["PLD_sym", "PLD_num"]