import mpmath
import numpy as np
[docs]
def q_b(m0, m1, m2, n0, n1, n2):
"""Stretch"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
return mpmath.sqrt((m0 - n0)**2 + (m1 - n1)**2 + (m2 - n2)**2)
[docs]
def dq_b(m0, m1, m2, n0, n1, n2):
"""Stretch, first derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = m0 - n0
x1 = m1 - n1
x2 = m2 - n2
x3 = 1/mpmath.sqrt(x0**2 + x1**2 + x2**2)
return np.array([x0*x3, x1*x3, x2*x3, -x0*x3, -x1*x3, -x2*x3], dtype=np.float64)
[docs]
def d2q_b(m0, m1, m2, n0, n1, n2):
"""Stretch, 2nd derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = m0 - n0
x1 = x0**2
x2 = m1 - n1
x3 = x2**2
x4 = m2 - n2
x5 = x4**2
x6 = x1 + x3 + x5
x7 = 1/mpmath.sqrt(x6)
x8 = x6**(mpmath.mpf(-3)/mpmath.mpf(2))
x9 = -x0
x10 = x8*x9
x11 = x0*x10 + x7
x12 = x10*x2
x13 = x10*x4
x14 = -x7
x15 = -x2
x16 = x10*x15
x17 = -x4
x18 = x10*x17
x19 = x15*x8
x20 = x0*x19
x21 = x19*x2 + x7
x22 = x19*x4
x23 = x17*x19
x24 = x17*x8
x25 = x0*x24
x26 = x2*x24
x27 = x24*x4 + x7
x28 = x0*x8
x29 = x2*x28
x30 = x28*x4
x31 = x2*x4*x8
return np.array([x11, x12, x13, x14 + x8*x9**2, x16, x18, x20, x21, x22, x16, x14 + x15**2*x8, x23, x25, x26, x27, x18, x23, x14 + x17**2*x8, x1*x8 + x14, x29, x30, x11, x20, x25, x29, x14 + x3*x8, x31, x12, x21, x26, x30, x31, x14 + x5*x8, x13, x22, x27], dtype=np.float64)
[docs]
def q_a(m0, m1, m2, o0, o1, o2, n0, n1, n2):
"""Bend"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -o0
x1 = m0 + x0
x2 = -o1
x3 = m1 + x2
x4 = -o2
x5 = m2 + x4
x6 = n0 + x0
x7 = n1 + x2
x8 = n2 + x4
return mpmath.acos((x1*x6 + x3*x7 + x5*x8)/(mpmath.sqrt(x1**2 + x3**2 + x5**2)*mpmath.sqrt(x6**2 + x7**2 + x8**2)))
[docs]
def dq_a(m0, m1, m2, o0, o1, o2, n0, n1, n2):
"""Bend, first derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -o0
x1 = m0 + x0
x2 = -o1
x3 = m1 + x2
x4 = -o2
x5 = m2 + x4
x6 = x1**2 + x3**2 + x5**2
x7 = n0 + x0
x8 = n1 + x2
x9 = n2 + x4
x10 = x7**2 + x8**2 + x9**2
x11 = x1*x7 + x3*x8 + x5*x9
x12 = 1/mpmath.sqrt(1 - x11**2/(x10*x6))
x13 = 1/mpmath.sqrt(x6)
x14 = 1/mpmath.sqrt(x10)
x15 = x13*x14
x16 = x11*x14/x6**(mpmath.mpf(3)/mpmath.mpf(2))
x17 = x11*x13/x10**(mpmath.mpf(3)/mpmath.mpf(2))
return np.array([-x12*(-x1*x16 + x15*x7), -x12*(x15*x8 - x16*x3), -x12*(x15*x9 - x16*x5), -x12*(x1*x16 + x15*(-m0 - n0 + 2*o0) + x17*x7), -x12*(x15*(-m1 - n1 + 2*o1) + x16*x3 + x17*x8), -x12*(x15*(-m2 - n2 + 2*o2) + x16*x5 + x17*x9), -x12*(x1*x15 - x17*x7), -x12*(x15*x3 - x17*x8), -x12*(x15*x5 - x17*x9)], dtype=np.float64)
[docs]
def d2q_a(m0, m1, m2, o0, o1, o2, n0, n1, n2):
"""Bend, 2nd derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -o0
x1 = m0 + x0
x2 = n0 + x0
x3 = x1*x2
x4 = -o1
x5 = m1 + x4
x6 = n1 + x4
x7 = x5*x6
x8 = -o2
x9 = m2 + x8
x10 = n2 + x8
x11 = x10*x9
x12 = x11 + x3 + x7
x13 = x12**2
x14 = x1**2
x15 = x5**2
x16 = x9**2
x17 = x14 + x15 + x16
x18 = x17**(-1)
x19 = x2**2
x20 = x6**2
x21 = x10**2
x22 = x19 + x20 + x21
x23 = x22**(-1)
x24 = x18*x23
x25 = -x13*x24 + 1
x26 = 1/mpmath.sqrt(x25)
x27 = -x1
x28 = 1/mpmath.sqrt(x22)
x29 = x17**(mpmath.mpf(-3)/mpmath.mpf(2))
x30 = x28*x29
x31 = x2*x30
x32 = x12*x30
x33 = -x32
x34 = x12*x27
x35 = -3*o0
x36 = 3*m0 + x35
x37 = x28/x17**(mpmath.mpf(5)/mpmath.mpf(2))
x38 = -x36*x37
x39 = 1/mpmath.sqrt(x17)
x40 = x28*x39
x41 = x2*x40 + x27*x32
x42 = x25**(mpmath.mpf(-3)/mpmath.mpf(2))
x43 = 2*n0
x44 = -2*o0
x45 = x43 + x44
x46 = (mpmath.mpf(1)/mpmath.mpf(2))*x45
x47 = x12*x24
x48 = 2*m0
x49 = x44 + x48
x50 = (mpmath.mpf(1)/mpmath.mpf(2))*x13
x51 = x23*x50/x17**2
x52 = x42*(x46*x47 - x49*x51)
x53 = -x5
x54 = x12*x38
x55 = x27*x30
x56 = x31*x53 + x55*x6
x57 = x32*x53 + x40*x6
x58 = -x9
x59 = x10*x55 + x31*x58
x60 = x10*x40 + x32*x58
x61 = -m0 - n0 - x44
x62 = x22**(mpmath.mpf(-3)/mpmath.mpf(2))
x63 = x39*x62
x64 = x12*x63
x65 = x1*x32 + x2*x64 + x40*x61
x66 = x29*x62
x67 = x34*x66
x68 = -x40
x69 = x32 + x68
x70 = x19*x63 + x2*x67 + x3*x30 + x55*x61 + x69
x71 = -2*o1
x72 = -m1 - n1 - x71
x73 = x2*x63
x74 = x6*x73
x75 = x31*x5 + x55*x72 + x6*x67 + x74
x76 = x32*x5 + x40*x72 + x6*x64
x77 = -2*o2
x78 = -m2 - n2 - x77
x79 = x10*x73
x80 = x10*x67 + x31*x9 + x55*x78 + x79
x81 = x10*x64 + x32*x9 + x40*x78
x82 = -x2
x83 = x26*(x1*x55 + x40 + x67*x82 + x73*x82)
x84 = x1*x40 + x64*x82
x85 = -x6
x86 = x26*(x5*x55 + x67*x85 + x73*x85)
x87 = x40*x5 + x64*x85
x88 = -x10
x89 = x26*(x55*x9 + x67*x88 + x73*x88)
x90 = x40*x9 + x64*x88
x91 = -3*o1
x92 = 3*m1 + x91
x93 = -x37*x92
x94 = 2*n1
x95 = x71 + x94
x96 = (mpmath.mpf(1)/mpmath.mpf(2))*x47
x97 = 2*m1
x98 = x71 + x97
x99 = x42*(-x51*x98 + x95*x96)
x100 = x30*x6
x101 = x12*x93
x102 = x30*x53
x103 = x10*x102 + x100*x58
x104 = x12*x53
x105 = x104*x66
x106 = x1*x100 + x102*x61 + x105*x2 + x74
x107 = x102*x72 + x105*x6 + x20*x63 + x30*x7 + x69
x108 = x6*x63
x109 = x10*x108
x110 = x10*x105 + x100*x9 + x102*x78 + x109
x111 = x26*(x1*x102 + x105*x82 + x108*x82)
x112 = x26*(x102*x5 + x105*x85 + x108*x85 + x40)
x113 = x26*(x102*x9 + x105*x88 + x108*x88)
x114 = -3*o2
x115 = 3*m2 + x114
x116 = -x115*x37
x117 = 2*n2
x118 = x117 + x77
x119 = 2*m2
x120 = x119 + x77
x121 = x42*(x118*x96 - x120*x51)
x122 = x10*x30
x123 = x116*x12
x124 = x30*x58
x125 = x12*x58
x126 = x125*x66
x127 = x1*x122 + x124*x61 + x126*x2 + x79
x128 = x109 + x122*x5 + x124*x72 + x126*x6
x129 = x10*x126 + x11*x30 + x124*x78 + x21*x63 + x69
x130 = x10*x63
x131 = x26*(x1*x124 + x126*x82 + x130*x82)
x132 = x26*(x124*x5 + x126*x85 + x130*x85)
x133 = x26*(x124*x9 + x126*x88 + x130*x88 + x40)
x134 = x18/x22**2
x135 = x42*(x13*x134*x46 + x49*x51 + x96*(4*o0 - x43 - x48))
x136 = x36*x37
x137 = x1*x30
x138 = 2*x61
x139 = x12*x136
x140 = x12*x2
x141 = 3*n0 + x35
x142 = x39/x22**(mpmath.mpf(5)/mpmath.mpf(2))
x143 = x141*x142
x144 = x12*x66
x145 = 2*x144
x146 = -x64
x147 = x146 + x33 + 2*x40
x148 = x12*x6
x149 = x30*x61
x150 = x1*x66
x151 = x140*x66
x152 = x108*x61 + x137*x72 + x148*x150 + x149*x5 + x151*x5 + x72*x73
x153 = x10*x12
x154 = x130*x61 + x137*x78 + x149*x9 + x150*x153 + x151*x9 + x73*x78
x155 = x12*x82
x156 = x61*x63
x157 = x64 + x68
x158 = x14*x30 + x150*x155 + x156*x82 + x157 + x3*x63
x159 = x12*x85
x160 = x137*x5
x161 = x150*x159 + x156*x85 + x160 + x5*x73
x162 = x12*x88
x163 = x137*x9
x164 = x150*x162 + x156*x88 + x163 + x73*x9
x165 = x37*x92
x166 = x134*x50
x167 = x42*(x166*x95 + x51*x98 + x96*(4*o1 - x94 - x97))
x168 = x12*x165
x169 = 3*n1 + x91
x170 = x142*x169
x171 = x30*x5
x172 = 2*x72
x173 = x30*x9
x174 = x5*x66
x175 = x66*x9
x176 = x108*x78 + x130*x72 + x148*x175 + x153*x174 + x171*x78 + x173*x72
x177 = x63*x72
x178 = x1*x108 + x155*x174 + x160 + x177*x82
x179 = x15*x30 + x157 + x159*x174 + x177*x85 + x63*x7
x180 = x171*x9
x181 = x5*x88
x182 = x108*x9 + x144*x181 + x177*x88 + x180
x183 = x115*x37
x184 = x42*(x118*x166 + x120*x51 + x96*(4*o2 - x117 - x119))
x185 = x12*x183
x186 = 3*n2 + x114
x187 = x142*x186
x188 = 2*x78
x189 = x63*x78
x190 = x1*x130 + x155*x175 + x163 + x189*x82
x191 = x85*x9
x192 = x130*x5 + x144*x191 + x180 + x189*x85
x193 = x11*x63 + x157 + x16*x30 + x162*x175 + x189*x88
x194 = x42*(-x166*x45 + x49*x96)
x195 = -x141*x142
x196 = x1*x63
x197 = x63*x82
x198 = x196*x85 + x197*x5
x199 = x196*x88 + x197*x9
x200 = x42*(-x166*x95 + x96*x98)
x201 = -x142*x169
x202 = 2*x63
x203 = x181*x63 + x191*x63
x204 = x42*(-x118*x166 + x120*x96)
x205 = -x142*x186
return np.array([-x26*(2*x27*x31 + x33 + x34*x38) - x41*x52, -x26*(x53*x54 + x56) - x52*x57, -x26*(x54*x58 + x59) - x52*x60, -x26*(x1*x54 + x70) - x52*x65, -x26*(x5*x54 + x75) - x52*x76, -x26*(x54*x9 + x80) - x52*x81, -x52*x84 - x83, -x52*x87 - x86, -x52*x90 - x89, -x26*(x34*x93 + x56) - x41*x99, -x26*(2*x100*x53 + x101*x53 + x33) - x57*x99, -x26*(x101*x58 + x103) - x60*x99, -x26*(x1*x101 + x106) - x65*x99, -x26*(x101*x5 + x107) - x76*x99, -x26*(x101*x9 + x110) - x81*x99, -x111 - x84*x99, -x112 - x87*x99, -x113 - x90*x99, -x121*x41 - x26*(x116*x34 + x59), -x121*x57 - x26*(x103 + x104*x116), -x121*x60 - x26*(2*x122*x58 + x123*x58 + x33), -x121*x65 - x26*(x1*x123 + x127), -x121*x76 - x26*(x123*x5 + x128), -x121*x81 - x26*(x123*x9 + x129), -x121*x84 - x131, -x121*x87 - x132, -x121*x90 - x133, -x135*x41 - x26*(x136*x34 + x70), -x135*x57 - x26*(x104*x136 + x106), -x135*x60 - x26*(x125*x136 + x127), -x135*x65 - x26*(x1*x139 + x137*x138 + x138*x73 + x140*x143 + x145*x3 + x147), -x135*x76 - x26*(x139*x5 + x143*x148 + x152), -x135*x81 - x26*(x139*x9 + x143*x153 + x154), -x135*x84 - x26*(x143*x155 + x158), -x135*x87 - x26*(x143*x159 + x161), -x135*x90 - x26*(x143*x162 + x164), -x167*x41 - x26*(x165*x34 + x75), -x167*x57 - x26*(x104*x165 + x107), -x167*x60 - x26*(x125*x165 + x128), -x167*x65 - x26*(x1*x168 + x140*x170 + x152), -x167*x76 - x26*(x108*x172 + x145*x7 + x147 + x148*x170 + x168*x5 + x171*x172), -x167*x81 - x26*(x153*x170 + x168*x9 + x176), -x167*x84 - x26*(x155*x170 + x178), -x167*x87 - x26*(x159*x170 + x179), -x167*x90 - x26*(x162*x170 + x182), -x184*x41 - x26*(x183*x34 + x80), -x184*x57 - x26*(x104*x183 + x110), -x184*x60 - x26*(x125*x183 + x129), -x184*x65 - x26*(x1*x185 + x140*x187 + x154), -x184*x76 - x26*(x148*x187 + x176 + x185*x5), -x184*x81 - x26*(x11*x145 + x130*x188 + x147 + x153*x187 + x173*x188 + x185*x9), -x184*x84 - x26*(x155*x187 + x190), -x184*x87 - x26*(x159*x187 + x192), -x184*x90 - x26*(x162*x187 + x193), -x194*x41 - x83, -x111 - x194*x57, -x131 - x194*x60, -x194*x65 - x26*(x140*x195 + x158), -x194*x76 - x26*(x148*x195 + x178), -x194*x81 - x26*(x153*x195 + x190), -x194*x84 - x26*(x146 + x155*x195 + 2*x196*x82), -x194*x87 - x26*(x159*x195 + x198), -x194*x90 - x26*(x162*x195 + x199), -x200*x41 - x86, -x112 - x200*x57, -x132 - x200*x60, -x200*x65 - x26*(x140*x201 + x161), -x200*x76 - x26*(x148*x201 + x179), -x200*x81 - x26*(x153*x201 + x192), -x200*x84 - x26*(x155*x201 + x198), -x200*x87 - x26*(x146 + x159*x201 + x202*x5*x85), -x200*x90 - x26*(x162*x201 + x203), -x204*x41 - x89, -x113 - x204*x57, -x133 - x204*x60, -x204*x65 - x26*(x140*x205 + x164), -x204*x76 - x26*(x148*x205 + x182), -x204*x81 - x26*(x153*x205 + x193), -x204*x84 - x26*(x155*x205 + x199), -x204*x87 - x26*(x159*x205 + x203), -x204*x90 - x26*(x146 + x162*x205 + x202*x88*x9)], dtype=np.float64)
[docs]
def q_a2(m0, m1, m2, o0, o1, o2, n0, n1, n2):
"""Bend2"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -o0
x1 = m0 + x0
x2 = -o1
x3 = n1 + x2
x4 = m1 + x2
x5 = n0 + x0
x6 = -o2
x7 = n2 + x6
x8 = m2 + x6
return mpmath.atan2(mpmath.sqrt((x1*x3 - x4*x5)**2 + (-x1*x7 + x5*x8)**2 + (-x3*x8 + x4*x7)**2), x1*x5 + x3*x4 + x7*x8)
[docs]
def dq_a2(m0, m1, m2, o0, o1, o2, n0, n1, n2):
"""Bend2, first derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -o0
x1 = n0 + x0
x2 = m0 + x0
x3 = -o1
x4 = n1 + x3
x5 = m1 + x3
x6 = -x1*x5 + x2*x4
x7 = -o2
x8 = n2 + x7
x9 = m2 + x7
x10 = x1*x9 - x2*x8
x11 = -x4*x9 + x5*x8
x12 = x10**2 + x11**2 + x6**2
x13 = mpmath.sqrt(x12)
x14 = x1*x2 + x4*x5 + x8*x9
x15 = (x12 + x14**2)**(-1)
x16 = x13*x15
x17 = 2*n1
x18 = -2*o1
x19 = x17 + x18
x20 = (mpmath.mpf(1)/mpmath.mpf(2))*x6
x21 = 2*n2
x22 = -2*o2
x23 = x21 + x22
x24 = (mpmath.mpf(1)/mpmath.mpf(2))*x10
x25 = x13**(-1)
x26 = 2*n0
x27 = -2*o0
x28 = x26 + x27
x29 = (mpmath.mpf(1)/mpmath.mpf(2))*x11
x30 = 2*m1
x31 = -x17 + x30
x32 = 2*m2
x33 = -x21 + x32
x34 = x14*x15*x25
x35 = 2*m0
x36 = -x26 + x35
x37 = x18 + x30
x38 = x22 + x32
x39 = x27 + x35
return np.array([-x1*x16 + x14*x15*x25*(x19*x20 - x23*x24), x14*x15*x25*(-x20*x28 + x23*x29) - x16*x4, x14*x15*x25*(-x19*x29 + x24*x28) - x16*x8, -x16*(-m0 - n0 - x27) + x34*(x20*x31 - x24*x33), -x16*(-m1 - n1 - x18) + x34*(-x20*x36 + x29*x33), -x16*(-m2 - n2 - x22) + x34*(x24*x36 - x29*x31), x14*x15*x25*(-x20*x37 + x24*x38) - x16*x2, x14*x15*x25*(x20*x39 - x29*x38) - x16*x5, x14*x15*x25*(-x24*x39 + x29*x37) - x16*x9], dtype=np.float64)
[docs]
def d2q_a2(m0, m1, m2, o0, o1, o2, n0, n1, n2):
"""Bend2, 2nd derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -o1
x1 = n1 + x0
x2 = -o2
x3 = n2 + x2
x4 = -x3
x5 = -o0
x6 = m0 + x5
x7 = n0 + x5
x8 = m1 + x0
x9 = m2 + x2
x10 = x1*x8 + x3*x9 + x6*x7
x11 = x1*x6
x12 = -x7
x13 = x12*x8
x14 = x11 + x13
x15 = x4*x6
x16 = x7*x9
x17 = x15 + x16
x18 = x3*x8
x19 = -x1
x20 = x19*x9
x21 = x18 + x20
x22 = x14**2 + x17**2 + x21**2
x23 = x10**2 + x22
x24 = x23**(-1)
x25 = mpmath.sqrt(x22)
x26 = x25**(-1)
x27 = x24*x26
x28 = x10*x27
x29 = 2*n1
x30 = -2*o1
x31 = x29 + x30
x32 = x14*x31
x33 = 2*n2
x34 = -2*o2
x35 = x33 + x34
x36 = -x17*x35
x37 = 2*n0
x38 = -2*o0
x39 = x37 + x38
x40 = -x10*x39 - x32 - x36
x41 = x23**(-2)
x42 = x25*x41
x43 = x40*x42
x44 = (mpmath.mpf(1)/mpmath.mpf(2))*x32 + (mpmath.mpf(1)/mpmath.mpf(2))*x36
x45 = x10*x44
x46 = x24/x22**(mpmath.mpf(3)/mpmath.mpf(2))
x47 = -x44*x46
x48 = x26*x41
x49 = x40*x48
x50 = x1*x12*x28
x51 = -x14*x39
x52 = x21*x35
x53 = (mpmath.mpf(1)/mpmath.mpf(2))*x51 + (mpmath.mpf(1)/mpmath.mpf(2))*x52
x54 = x27*x7
x55 = x53*x54
x56 = x27*x44
x57 = x1*x56
x58 = x10*x53
x59 = x28*x4*x7
x60 = x17*x39
x61 = -x21*x31
x62 = (mpmath.mpf(1)/mpmath.mpf(2))*x60 + (mpmath.mpf(1)/mpmath.mpf(2))*x61
x63 = x54*x62
x64 = x3*x56
x65 = x10*x62
x66 = 2*m1
x67 = -x29 + x66
x68 = x14*x67
x69 = 2*m2
x70 = -x33 + x69
x71 = -x17*x70
x72 = (mpmath.mpf(1)/mpmath.mpf(2))*x68 + (mpmath.mpf(1)/mpmath.mpf(2))*x71
x73 = x54*x72
x74 = -m0 - n0 - x38
x75 = x56*x74
x76 = x10*x72
x77 = x24*x25
x78 = m1 - n1
x79 = m2 - n2
x80 = -x79
x81 = x28*(x1*x78 + x4*x80) + x77
x82 = 2*m0
x83 = -x37 + x82
x84 = -x14*x83
x85 = x21*x70
x86 = (mpmath.mpf(1)/mpmath.mpf(2))*x84 + (mpmath.mpf(1)/mpmath.mpf(2))*x85
x87 = x54*x86
x88 = m0 - n0
x89 = -x88
x90 = x28*(x1*x89 - x14)
x91 = -m1 - n1 - x30
x92 = x56*x91
x93 = x10*x86
x94 = x28*(x17 + x4*x88)
x95 = x17*x83
x96 = -x21*x67
x97 = (mpmath.mpf(1)/mpmath.mpf(2))*x95 + (mpmath.mpf(1)/mpmath.mpf(2))*x96
x98 = x54*x97
x99 = -m2 - n2 - x34
x100 = x56*x99
x101 = x10*x97
x102 = x30 + x66
x103 = -x102*x14
x104 = x34 + x69
x105 = x104*x17
x106 = (mpmath.mpf(1)/mpmath.mpf(2))*x103 + (mpmath.mpf(1)/mpmath.mpf(2))*x105
x107 = x106*x54
x108 = x56*x6
x109 = x10*x106
x110 = -x77
x111 = -x8
x112 = x110 + x28*(x1*x111 + x4*x9)
x113 = x28*(2*x11 + x13)
x114 = x38 + x82
x115 = x114*x14
x116 = -x104*x21
x117 = (mpmath.mpf(1)/mpmath.mpf(2))*x115 + (mpmath.mpf(1)/mpmath.mpf(2))*x116
x118 = x117*x54
x119 = x56*x8
x120 = x10*x117
x121 = -x114*x17
x122 = x102*x21
x123 = (mpmath.mpf(1)/mpmath.mpf(2))*x121 + (mpmath.mpf(1)/mpmath.mpf(2))*x122
x124 = x123*x54
x125 = -x6
x126 = x28*(x125*x4 - x17)
x127 = x56*x9
x128 = x10*x123
x129 = -x10*x31 - x51 - x52
x130 = x129*x42
x131 = -x46*x53
x132 = x129*x48
x133 = x19*x28*x3
x134 = x1*x27
x135 = x134*x62
x136 = x27*x53
x137 = x136*x3
x138 = x28*(x12*x78 + x14)
x139 = x134*x72
x140 = x136*x74
x141 = x134*x86
x142 = x136*x91
x143 = x28*(x12*x89 + x3*x79) + x77
x144 = x134*x97
x145 = -x78
x146 = x28*(x145*x3 - x21)
x147 = x136*x99
x148 = x106*x134
x149 = x28*(x111*x12 - x14)
x150 = x136*x6
x151 = x117*x134
x152 = x136*x8
x153 = -x9
x154 = x110 + x28*(x12*x6 + x153*x3)
x155 = x28*(2*x18 + x20)
x156 = x123*x134
x157 = x136*x9
x158 = -x10*x35 - x60 - x61
x159 = x158*x42
x160 = -x46*x62
x161 = x158*x48
x162 = x27*x3
x163 = x162*x72
x164 = x28*(-x17 + x7*x80)
x165 = x27*x62
x166 = x165*x74
x167 = x28*(x19*x79 + x21)
x168 = x162*x86
x169 = x165*x91
x170 = x162*x97
x171 = x165*x99
x172 = x28*(x145*x19 + x7*x88) + x77
x173 = x28*(x15 + 2*x16)
x174 = x106*x162
x175 = x165*x6
x176 = x117*x162
x177 = x28*(x153*x19 - x21)
x178 = x165*x8
x179 = x123*x162
x180 = x165*x9
x181 = x110 + x28*(x125*x7 + x19*x8)
x182 = -x10*(4*o0 - x37 - x82) - x68 - x71
x183 = x182*x42
x184 = -x46*x72
x185 = x182*x48
x186 = -2*x77
x187 = x28*x78*x89
x188 = x27*x74
x189 = x188*x86
x190 = x27*x72
x191 = x190*x91
x192 = x28*x80*x88
x193 = x188*x97
x194 = x190*x99
x195 = x106*x188
x196 = x190*x6
x197 = x28*(x111*x78 + x80*x9) + x77
x198 = x28*(-x14 + x6*x78)
x199 = x190*x8
x200 = x117*x188
x201 = x28*(x125*x80 + x17)
x202 = x190*x9
x203 = x123*x188
x204 = -x10*(4*o1 - x29 - x66) - x84 - x85
x205 = x204*x42
x206 = -x46*x86
x207 = x204*x48
x208 = x145*x28*x79
x209 = x27*x91
x210 = x209*x97
x211 = x27*x86
x212 = x211*x99
x213 = x28*(x111*x89 + x14)
x214 = x211*x6
x215 = x106*x209
x216 = x117*x209
x217 = x211*x8
x218 = x28*(x153*x79 + x6*x89) + x77
x219 = x28*(-x21 + x79*x8)
x220 = x211*x9
x221 = x123*x209
x222 = -x10*(4*o2 - x33 - x69) - x95 - x96
x223 = x222*x42
x224 = -x46*x97
x225 = x222*x48
x226 = x28*(-x17 + x88*x9)
x227 = x27*x97
x228 = x227*x6
x229 = x27*x99
x230 = x106*x229
x231 = x28*(x145*x153 + x21)
x232 = x227*x8
x233 = x117*x229
x234 = x123*x229
x235 = x227*x9
x236 = x28*(x125*x88 + x145*x8) + x77
x237 = -x10*x114 - x103 - x105
x238 = x237*x42
x239 = -x106*x46
x240 = x237*x48
x241 = x111*x28*x6
x242 = x27*x6
x243 = x117*x242
x244 = x106*x27
x245 = x244*x8
x246 = x125*x28*x9
x247 = x123*x242
x248 = x244*x9
x249 = -x10*x102 - x115 - x116
x250 = x249*x42
x251 = -x117*x46
x252 = x249*x48
x253 = x153*x28*x8
x254 = x123*x27*x8
x255 = x117*x27*x9
x256 = -x10*x104 - x121 - x122
x257 = x256*x42
x258 = -x123*x46
x259 = x256*x48
return np.array([x28*(x1**2 + x4**2) - x43*x7 + x45*x47 + x45*x49, -x1*x43 + x47*x58 + x49*x58 + x50 + x55 - x57, -x3*x43 + x47*x65 + x49*x65 + x59 + x63 - x64, -x43*x74 + x47*x76 + x49*x76 + x73 - x75 + x81, -x43*x91 + x47*x93 + x49*x93 + x87 + x90 - x92, -x100 + x101*x47 + x101*x49 - x43*x99 + x94 + x98, x107 - x108 + x109*x47 + x109*x49 + x112 - x43*x6, x113 + x118 - x119 + x120*x47 + x120*x49 - x43*x8, x124 + x126 - x127 + x128*x47 + x128*x49 - x43*x9, -x130*x7 + x131*x45 + x132*x45 + x50 - x55 + x57, -x1*x130 + x131*x58 + x132*x58 + x28*(x12**2 + x3**2), -x130*x3 + x131*x65 + x132*x65 + x133 + x135 - x137, -x130*x74 + x131*x76 + x132*x76 + x138 + x139 - x140, -x130*x91 + x131*x93 + x132*x93 + x141 - x142 + x143, x101*x131 + x101*x132 - x130*x99 + x144 + x146 - x147, x109*x131 + x109*x132 - x130*x6 + x148 + x149 - x150, x120*x131 + x120*x132 - x130*x8 + x151 - x152 + x154, x128*x131 + x128*x132 - x130*x9 + x155 + x156 - x157, -x159*x7 + x160*x45 + x161*x45 + x59 - x63 + x64, -x1*x159 + x133 - x135 + x137 + x160*x58 + x161*x58, -x159*x3 + x160*x65 + x161*x65 + x28*(x19**2 + x7**2), -x159*x74 + x160*x76 + x161*x76 + x163 + x164 - x166, -x159*x91 + x160*x93 + x161*x93 + x167 + x168 - x169, x101*x160 + x101*x161 - x159*x99 + x170 - x171 + x172, x109*x160 + x109*x161 - x159*x6 + x173 + x174 - x175, x120*x160 + x120*x161 - x159*x8 + x176 + x177 - x178, x128*x160 + x128*x161 - x159*x9 + x179 - x180 + x181, -x183*x7 + x184*x45 + x185*x45 - x73 + x75 + x81, -x1*x183 + x138 - x139 + x140 + x184*x58 + x185*x58, -x163 + x164 + x166 - x183*x3 + x184*x65 + x185*x65, -x183*x74 + x184*x76 + x185*x76 + x186 + x28*(x78**2 + x80**2), -x183*x91 + x184*x93 + x185*x93 + x187 + x189 - x191, x101*x184 + x101*x185 - x183*x99 + x192 + x193 - x194, x109*x184 + x109*x185 - x183*x6 + x195 - x196 + x197, x120*x184 + x120*x185 - x183*x8 + x198 - x199 + x200, x128*x184 + x128*x185 - x183*x9 + x201 - x202 + x203, -x205*x7 + x206*x45 + x207*x45 - x87 + x90 + x92, -x1*x205 - x141 + x142 + x143 + x206*x58 + x207*x58, x167 - x168 + x169 - x205*x3 + x206*x65 + x207*x65, x187 - x189 + x191 - x205*x74 + x206*x76 + x207*x76, x186 - x205*x91 + x206*x93 + x207*x93 + x28*(x79**2 + x89**2), x101*x206 + x101*x207 - x205*x99 + x208 + x210 - x212, x109*x206 + x109*x207 - x205*x6 + x213 - x214 + x215, x120*x206 + x120*x207 - x205*x8 + x216 - x217 + x218, x128*x206 + x128*x207 - x205*x9 + x219 - x220 + x221, x100 - x223*x7 + x224*x45 + x225*x45 + x94 - x98, -x1*x223 - x144 + x146 + x147 + x224*x58 + x225*x58, -x170 + x171 + x172 - x223*x3 + x224*x65 + x225*x65, x192 - x193 + x194 - x223*x74 + x224*x76 + x225*x76, x208 - x210 + x212 - x223*x91 + x224*x93 + x225*x93, x101*x224 + x101*x225 + x186 - x223*x99 + x28*(x145**2 + x88**2), x109*x224 + x109*x225 - x223*x6 + x226 - x228 + x230, x120*x224 + x120*x225 - x223*x8 + x231 - x232 + x233, x128*x224 + x128*x225 - x223*x9 + x234 - x235 + x236, -x107 + x108 + x112 - x238*x7 + x239*x45 + x240*x45, -x1*x238 - x148 + x149 + x150 + x239*x58 + x240*x58, x173 - x174 + x175 - x238*x3 + x239*x65 + x240*x65, -x195 + x196 + x197 - x238*x74 + x239*x76 + x240*x76, x213 + x214 - x215 - x238*x91 + x239*x93 + x240*x93, x101*x239 + x101*x240 + x226 + x228 - x230 - x238*x99, x109*x239 + x109*x240 - x238*x6 + x28*(x111**2 + x9**2), x120*x239 + x120*x240 - x238*x8 + x241 + x243 - x245, x128*x239 + x128*x240 - x238*x9 + x246 + x247 - x248, x113 - x118 + x119 - x250*x7 + x251*x45 + x252*x45, -x1*x250 - x151 + x152 + x154 + x251*x58 + x252*x58, -x176 + x177 + x178 - x250*x3 + x251*x65 + x252*x65, x198 + x199 - x200 - x250*x74 + x251*x76 + x252*x76, -x216 + x217 + x218 - x250*x91 + x251*x93 + x252*x93, x101*x251 + x101*x252 + x231 + x232 - x233 - x250*x99, x109*x251 + x109*x252 + x241 - x243 + x245 - x250*x6, x120*x251 + x120*x252 - x250*x8 + x28*(x153**2 + x6**2), x128*x251 + x128*x252 - x250*x9 + x253 + x254 - x255, -x124 + x126 + x127 - x257*x7 + x258*x45 + x259*x45, -x1*x257 + x155 - x156 + x157 + x258*x58 + x259*x58, -x179 + x180 + x181 - x257*x3 + x258*x65 + x259*x65, x201 + x202 - x203 - x257*x74 + x258*x76 + x259*x76, x219 + x220 - x221 - x257*x91 + x258*x93 + x259*x93, x101*x258 + x101*x259 - x234 + x235 + x236 - x257*x99, x109*x258 + x109*x259 + x246 - x247 + x248 - x257*x6, x120*x258 + x120*x259 + x253 - x254 + x255 - x257*x8, x128*x258 + x128*x259 - x257*x9 + x28*(x125**2 + x8**2)], dtype=np.float64)
[docs]
def q_d(m0, m1, m2, o0, o1, o2, p0, p1, p2, n0, n1, n2):
"""Torsion"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
p0 = mpmath.mpf(p0)
p1 = mpmath.mpf(p1)
p2 = mpmath.mpf(p2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -p0
x1 = -o0 - x0
x2 = m0 - o0
x3 = m1 - o1
x4 = m2 - o2
x5 = 1/mpmath.sqrt(x2**2 + x3**2 + x4**2)
x6 = -p1
x7 = -o1 - x6
x8 = -p2
x9 = -o2 - x8
x10 = 1/mpmath.sqrt(x1**2 + x7**2 + x9**2)
x11 = x10*x5
x12 = x11*x2
x13 = x11*x3
x14 = x11*x4
x15 = n0 + x0
x16 = n1 + x6
x17 = n2 + x8
x18 = 1/mpmath.sqrt(x15**2 + x16**2 + x17**2)
x19 = x10*x18
x20 = x15*x19
x21 = x16*x19
x22 = x17*x19
return mpmath.acos(((-x1*x13 + x12*x7)*(-x1*x21 + x20*x7) + (x13*x9 - x14*x7)*(x21*x9 - x22*x7) + (x1*x10*x17*x18 - x20*x9)*(x1*x10*x4*x5 - x12*x9))/(mpmath.sqrt(1 - (x1*x12 + x13*x7 + x14*x9)**2)*mpmath.sqrt(1 - (-x1*x20 - x21*x7 - x22*x9)**2)))
[docs]
def dq_d(m0, m1, m2, o0, o1, o2, p0, p1, p2, n0, n1, n2):
"""Torsion, first derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
p0 = mpmath.mpf(p0)
p1 = mpmath.mpf(p1)
p2 = mpmath.mpf(p2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = m0 - o0
x1 = -p0
x2 = o0 + x1
x3 = -x2
x4 = x0**2
x5 = m1 - o1
x6 = x5**2
x7 = m2 - o2
x8 = x7**2
x9 = x4 + x6 + x8
x10 = 1/mpmath.sqrt(x9)
x11 = x3**2
x12 = -p1
x13 = o1 + x12
x14 = -x13
x15 = x14**2
x16 = -p2
x17 = o2 + x16
x18 = -x17
x19 = x18**2
x20 = x11 + x15 + x19
x21 = 1/mpmath.sqrt(x20)
x22 = x10*x21
x23 = x22*x3
x24 = x14*x22
x25 = x18*x22
x26 = x0*x23 + x24*x5 + x25*x7
x27 = 1 - x26**2
x28 = n0 + x1
x29 = x28**2
x30 = n1 + x12
x31 = x30**2
x32 = n2 + x16
x33 = x32**2
x34 = x29 + x31 + x33
x35 = 1/mpmath.sqrt(x34)
x36 = x21*x35
x37 = x28*x36
x38 = x30*x36
x39 = x32*x36
x40 = -x14*x38 - x18*x39 - x3*x37
x41 = 1 - x40**2
x42 = x0*x24 - x23*x5
x43 = x14*x37 - x3*x38
x44 = -x0*x25 + x10*x21*x3*x7
x45 = -x18*x37 + x21*x3*x32*x35
x46 = -x24*x7 + x25*x5
x47 = -x14*x39 + x18*x38
x48 = x42*x43 + x44*x45 + x46*x47
x49 = 1/mpmath.sqrt(1 - x48**2/(x27*x41))
x50 = -x0
x51 = x9**(mpmath.mpf(-3)/mpmath.mpf(2))
x52 = x21*x51
x53 = x50*x52
x54 = x5*x53
x55 = x53*x7
x56 = x0*x14
x57 = x0*x18
x58 = 1/mpmath.sqrt(x27)
x59 = 1/mpmath.sqrt(x41)
x60 = x58*x59
x61 = 2*x23
x62 = x0*x3
x63 = 2*x62
x64 = 2*x14
x65 = 2*x18
x66 = (mpmath.mpf(1)/mpmath.mpf(2))*x48
x67 = x26*x59*x66/x27**(mpmath.mpf(3)/mpmath.mpf(2))
x68 = -x5
x69 = x52*x68
x70 = x5*x69
x71 = x69*x7
x72 = 2*x24
x73 = -x7
x74 = x52*x73
x75 = x5*x74
x76 = x7*x74
x77 = 2*x25
x78 = 2*x37
x79 = x20**(mpmath.mpf(-3)/mpmath.mpf(2))
x80 = x11*x79
x81 = x35*x80
x82 = 2*x28
x83 = x3*x79
x84 = x35*x83
x85 = x14*x84
x86 = x30*x85
x87 = x18*x84
x88 = x32*x87
x89 = x40*x58*x66/x41**(mpmath.mpf(3)/mpmath.mpf(2))
x90 = x0*x22
x91 = 2*x90
x92 = x10*x80
x93 = 2*x0
x94 = x4*x52
x95 = 2*x3
x96 = x5*x52
x97 = x56*x96
x98 = x52*x7
x99 = x57*x98
x100 = x10*x83
x101 = x100*x5
x102 = x101*x14
x103 = x100*x7
x104 = x30*x87
x105 = -x32*x85
x106 = x14*x28
x107 = x106*x84
x108 = x28*x87
x109 = x101*x18
x110 = x56*x98
x111 = x14*x7
x112 = -x100*x111 + x57*x96
x113 = x22*x5
x114 = x100*x56
x115 = x62*x96
x116 = x100*x57
x117 = x22*x7
x118 = x117 - x25
x119 = 2*x38
x120 = x15*x79
x121 = x120*x35
x122 = 2*x30
x123 = x18*x79
x124 = x123*x35
x125 = 2*x113
x126 = x10*x120
x127 = x52*x6
x128 = x7*x96
x129 = x128*x18
x130 = x10*x123
x131 = x111*x130
x132 = x106*x124
x133 = x14*x30
x134 = x124*x133
x135 = -x128*x3 + x130*x56
x136 = x23 - x90
x137 = x130*x5
x138 = x137*x14
x139 = x111*x96
x140 = 2*x39
x141 = x19*x79
x142 = x141*x35
x143 = 2*x32
x144 = 2*x117
x145 = x10*x141
x146 = 2*x7
x147 = x52*x8
x148 = x14*x32
x149 = x100*x93
x150 = x2*x79
x151 = x10*x150
x152 = x5*x64
x153 = x130*x146
x154 = x3*x36
x155 = 2*x154
x156 = x34**(mpmath.mpf(-3)/mpmath.mpf(2))
x157 = x156*x21
x158 = x157*x29
x159 = x157*x32
x160 = x159*x18
x161 = x82*x84
x162 = x150*x35
x163 = 2*x133
x164 = x124*x143
x165 = x157*x30
x166 = x165*x28
x167 = x166*x18
x168 = x2*x30
x169 = x106*x159
x170 = x159*x28
x171 = x124*x28
x172 = x18*x36
x173 = x172 + x39
x174 = x14*x36
x175 = x174 + x38
x176 = x13*x79
x177 = x10*x176
x178 = 2*x174
x179 = x157*x31
x180 = x176*x35
x181 = x165*x32
x182 = -x181*x3
x183 = x154 + x37
x184 = x17*x79
x185 = x10*x184
x186 = 2*x172
x187 = x157*x33
x188 = x184*x35
x189 = x17*x30
x190 = -x28
x191 = x165*x190
x192 = x157*x190
x193 = x18*x28
x194 = x28*x95
x195 = x159*x65
x196 = -x30
x197 = x157*x196
x198 = x165*x196
x199 = -x32
x200 = x157*x199
return np.array([-x49*(x60*(x43*(x24 - x3*x54 + x53*x56) + x45*(x21*x3*x50*x51*x7 - x25 - x53*x57) + x47*(-x14*x55 + x18*x54)) + x67*(x53*x63 + x54*x64 + x55*x65 + x61)), -x49*(x60*(x43*(x0*x14*x21*x51*x68 - x23 - x3*x70) + x45*(x21*x3*x51*x68*x7 - x57*x69) + x47*(-x14*x71 + x18*x70 + x25)) + x67*(x63*x69 + x64*x70 + x65*x71 + x72)), -x49*(x60*(x43*(-x3*x75 + x56*x74) + x45*(x23 + x3*x76 - x57*x74) + x47*(-x14*x76 + x18*x21*x5*x51*x73 - x24)) + x67*(x63*x74 + x64*x75 + x65*x76 + x77)), -x49*(x60*(x42*(x107 - x30*x81 + x38) + x43*(x113 + x114 - x115 + x14*x94 - x24 - x5*x92) + x44*(-x108 + x11*x32*x35*x79 - x39) + x45*(x0*x21*x3*x51*x7 + x10*x11*x7*x79 - x116 - x118 - x18*x94) + x46*(x104 + x105) + x47*(x109 - x110 + x112)) + x67*(2*x102 + 2*x103*x18 - x61 - x91 + x92*x93 + x94*x95 + 2*x97 + 2*x99) + x89*(x78 - x81*x82 - 2*x86 - 2*x88)), -x49*(x60*(x42*(x15*x28*x35*x79 - x37 - x86) + x43*(x0*x126 - x102 - x127*x3 + x136 + x97) + x44*(-x105 - x132) + x45*(-x112 - x135) + x46*(-x121*x32 + x134 + x39) + x47*(x118 - x126*x7 + x127*x18 + x138 - x139)) + x67*(2*x114 + 2*x115 - x125 + 2*x126*x5 + x127*x64 + 2*x129 + 2*x131 - x72) + x89*(-2*x107 + x119 - x121*x122 - x124*x32*x64)), -x49*(x60*(x42*(-x104 + x132) + x43*(-x109 + x110 + x135) + x44*(-x142*x28 + x37 + x88) + x45*(-x0*x145 + x10*x18*x3*x7*x79 - x136 + x21*x3*x51*x8 - x99) + x46*(-x124*x148 + x19*x30*x35*x79 - x38) + x47*(-x113 + x129 - x131 - x14*x147 + x145*x5 + x24)) + x67*(2*x116 + 2*x138 + 2*x139 - x144 + x145*x146 + x147*x65 + 2*x62*x98 - x77) + x89*(-2*x108 - 2*x134 + x140 - x142*x143)), -x49*(x60*(x42*(x14*x156*x21*x29 + x14*x2*x28*x35*x79 - x166*x3 - x168*x84 - x175) + x43*(x0*x10*x14*x2*x79 - x101*x2 - x113) + x44*(-x158*x18 + x170*x3 - x171*x2 + x173 + x2*x32*x84) + x45*(x103*x2 + x117 - x151*x57) + x46*(x124*x168 - x148*x162 + x167 - x169) + x47*(-x111*x151 + x137*x2)) + x67*(x149*x2 + x151*x152 + x153*x2 + x91) + x89*(-x106*x122*x157 + x155 - x158*x95 - x160*x82 - x161*x2 - x162*x163 - x164*x2 - x78)), -x49*(x60*(x42*(x106*x165 + x106*x180 - x13*x30*x84 - x179*x3 + x183) + x43*(-x101*x13 + x177*x56 + x90) + x44*(-x13*x171 + x13*x3*x32*x35*x79 - x167 - x182) + x45*(x10*x13*x3*x7*x79 - x177*x57) + x46*(x13*x18*x30*x35*x79 - x133*x159 - x148*x180 + x156*x18*x21*x31 - x173) + x47*(x10*x13*x18*x5*x79 - x111*x177 - x117)) + x67*(x125 + x13*x149 + x13*x153 + x152*x177) + x89*(-x119 - x122*x160 - x13*x161 - x13*x164 - x163*x180 - x166*x95 + x178 - x179*x64)), -x49*(x60*(x42*(x106*x188 + x169 + x182 - x189*x84) + x43*(-x101*x17 + x185*x56) + x44*(x156*x21*x3*x33 - x17*x171 + x17*x3*x32*x35*x79 - x170*x18 - x183) + x45*(x10*x17*x3*x7*x79 - x185*x57 - x90) + x46*(x124*x189 - x14*x187 - x148*x188 + x175 + x18*x181) + x47*(-x111*x185 + x113 + x137*x17)) + x67*(x144 + x149*x17 + x152*x185 + x153*x17) + x89*(-x133*x143*x157 - x140 - x161*x17 - x163*x188 - x164*x17 - x170*x95 + x186 - x187*x65)), -x49*(x60*(x42*(x106*x192 + x174 - x191*x3) + x44*(x156*x190*x21*x3*x32 - x172 - x192*x193) + x46*(-x148*x192 + x18*x191)) + x89*(-x155 - x163*x192 - x190*x195 - x192*x194)), -x49*(x60*(x42*(x14*x156*x196*x21*x28 - x154 - x198*x3) + x44*(x156*x196*x21*x3*x32 - x193*x197) + x46*(-x148*x197 + x172 + x18*x198)) + x89*(-x163*x197 - x178 - x194*x197 - x195*x196)), -x49*(x60*(x42*(x106*x200 - x165*x199*x3) + x44*(x154 + x159*x199*x3 - x193*x200) + x46*(-x148*x200 + x156*x18*x199*x21*x30 - x174)) + x89*(-x163*x200 - x186 - x194*x200 - x195*x199))], dtype=np.float64)
[docs]
def d2q_d(m0, m1, m2, o0, o1, o2, p0, p1, p2, n0, n1, n2):
"""Torsion, 2nd derivative wrt. Cartesians"""
m0 = mpmath.mpf(m0)
m1 = mpmath.mpf(m1)
m2 = mpmath.mpf(m2)
o0 = mpmath.mpf(o0)
o1 = mpmath.mpf(o1)
o2 = mpmath.mpf(o2)
p0 = mpmath.mpf(p0)
p1 = mpmath.mpf(p1)
p2 = mpmath.mpf(p2)
n0 = mpmath.mpf(n0)
n1 = mpmath.mpf(n1)
n2 = mpmath.mpf(n2)
x0 = -p2
x1 = o2 + x0
x2 = -x1
x3 = -p1
x4 = n1 + x3
x5 = -p0
x6 = n0 + x5
x7 = x6**2
x8 = x4**2
x9 = n2 + x0
x10 = x9**2
x11 = x10 + x7 + x8
x12 = 1/mpmath.sqrt(x11)
x13 = o0 + x5
x14 = -x13
x15 = x14**2
x16 = o1 + x3
x17 = -x16
x18 = x17**2
x19 = x2**2
x20 = x15 + x18 + x19
x21 = 1/mpmath.sqrt(x20)
x22 = x12*x21
x23 = x22*x4
x24 = x22*x9
x25 = -x17*x24 + x2*x23
x26 = m0 - o0
x27 = -x26
x28 = m1 - o1
x29 = x26**2
x30 = x28**2
x31 = m2 - o2
x32 = x31**2
x33 = x29 + x30 + x32
x34 = x33**(mpmath.mpf(-3)/mpmath.mpf(2))
x35 = x21*x34
x36 = x28*x35
x37 = x2*x36
x38 = x17*x35
x39 = x31*x38
x40 = x27*x37 - x27*x39
x41 = x25*x40
x42 = x22*x6
x43 = -x14*x23 + x17*x42
x44 = 1/mpmath.sqrt(x33)
x45 = x21*x44
x46 = x17*x45
x47 = x26*x38
x48 = x14*x36
x49 = x27*x47 - x27*x48 + x46
x50 = x43*x49
x51 = x12*x14*x21*x9 - x2*x42
x52 = x2*x45
x53 = x2*x35
x54 = x26*x53
x55 = x14*x21*x27*x31*x34 - x27*x54 - x52
x56 = x51*x55
x57 = x41 + x50 + x56
x58 = x14*x45
x59 = x26*x58 + x28*x46 + x31*x52
x60 = x59**2
x61 = 1 - x60
x62 = 1/mpmath.sqrt(x61)
x63 = -x14*x42 - x17*x23 - x2*x24
x64 = x63**2
x65 = 1 - x64
x66 = 1/mpmath.sqrt(x65)
x67 = x62*x66
x68 = 2*x58
x69 = x14*x35
x70 = x27*x69
x71 = x26*x70
x72 = x27*x36
x73 = x17*x72
x74 = x27*x53
x75 = 2*x74
x76 = x31*x75 + x68 + 2*x71 + 2*x73
x77 = x59*x76
x78 = x61**(mpmath.mpf(-3)/mpmath.mpf(2))
x79 = x66*x78
x80 = x77*x79
x81 = x26*x46 - x28*x58
x82 = x14*x21*x31*x44 - x26*x52
x83 = x28*x52 - x31*x46
x84 = x25*x83 + x43*x81 + x51*x82
x85 = (mpmath.mpf(1)/mpmath.mpf(2))*x84
x86 = x57*x67 + x80*x85
x87 = x84**2
x88 = x61**(-1)
x89 = x65**(-1)
x90 = x88*x89
x91 = -x87*x90 + 1
x92 = x91**(mpmath.mpf(-3)/mpmath.mpf(2))
x93 = (mpmath.mpf(1)/mpmath.mpf(2))*x87
x94 = x89*x93/x61**2
x95 = x85*x90
x96 = x92*(x77*x94 + x95*(2*x41 + 2*x50 + 2*x56))
x97 = 1/mpmath.sqrt(x91)
x98 = (mpmath.mpf(3)/mpmath.mpf(4))*x84
x99 = x60*x66*x98/x61**(mpmath.mpf(5)/mpmath.mpf(2))
x100 = x79*x85
x101 = x100*(x31*x74 + x58 + x71 + x73)
x102 = 3*o0
x103 = 3*m0 - x102
x104 = -x103
x105 = x33**(mpmath.mpf(-5)/mpmath.mpf(2))
x106 = x105*x21
x107 = x104*x106
x108 = x107*x27
x109 = x108*x28
x110 = x108*x31
x111 = -x37
x112 = x111 + x39
x113 = x27*x38
x114 = 2*x113
x115 = x17*x26
x116 = -x47
x117 = x116 + x48
x118 = x2*x26
x119 = x31*x69
x120 = -x54
x121 = x119 + x120
x122 = 2*x14
x123 = x122*x26
x124 = 2*x17
x125 = 2*x2
x126 = x26*x69
x127 = 2*x126
x128 = x17*x36
x129 = 2*x128
x130 = x31*x53
x131 = 2*x130
x132 = -x127 - x129 - x131
x133 = x59*x79
x134 = x133*x85
x135 = -x28
x136 = x135*x14*x21*x31*x34 - x135*x54
x137 = x136*x51
x138 = x135*x37 - x135*x39 + x52
x139 = x138*x25
x140 = x135*x17*x21*x26*x34 - x135*x48 - x58
x141 = x140*x43
x142 = x137 + x139 + x141
x143 = 2*x46
x144 = x135*x69
x145 = 2*x144
x146 = x135*x53
x147 = x146*x31
x148 = x129*x135 + x143 + x145*x26 + 2*x147
x149 = x133*x148
x150 = x142*x67 + x149*x85
x151 = x107*x135
x152 = x151*x28
x153 = x151*x31
x154 = x135*x38
x155 = -x70
x156 = x154 + x155
x157 = x114 + x145
x158 = (mpmath.mpf(1)/mpmath.mpf(2))*x57
x159 = (mpmath.mpf(1)/mpmath.mpf(2))*x80
x160 = x76*x99
x161 = x142*x159 + x148*x160 + x149*x158
x162 = -x31
x163 = x162*x47 - x162*x48
x164 = x163*x43
x165 = x119*x162 - x162*x54 + x58
x166 = x165*x51
x167 = x162*x2*x21*x28*x34 - x162*x39 - x46
x168 = x167*x25
x169 = x164 + x166 + x168
x170 = 2*x52
x171 = x162*x69
x172 = 2*x171
x173 = x162*x53
x174 = x173*x31
x175 = x129*x162 + x170 + x172*x26 + 2*x174
x176 = x133*x175
x177 = x169*x67 + x176*x85
x178 = x162*x38
x179 = x107*x162
x180 = x179*x28
x181 = x179*x31
x182 = x155 + x173
x183 = x172 + x75
x184 = x158*x176 + x159*x169 + x160*x175
x185 = 2*x42
x186 = x20**(mpmath.mpf(-3)/mpmath.mpf(2))
x187 = x12*x186
x188 = x15*x187
x189 = x188*x6
x190 = x187*x4
x191 = x14*x190
x192 = x17*x191
x193 = x187*x9
x194 = x14*x193
x195 = x194*x2
x196 = x185 - 2*x189 - 2*x192 - 2*x195
x197 = x196*x63
x198 = x65**(mpmath.mpf(-3)/mpmath.mpf(2))
x199 = x198*x62
x200 = x197*x199
x201 = 2*x45
x202 = x201*x26
x203 = x186*x44
x204 = x15*x203
x205 = 2*x204
x206 = x29*x35
x207 = x14*x206
x208 = x31*x54
x209 = x14*x203
x210 = x17*x209
x211 = 2*x210
x212 = x2*x209
x213 = 2*x212
x214 = x129*x26 - x202 + x205*x26 + 2*x207 + 2*x208 + x211*x28 + x213*x31 - x68
x215 = x133*x214
x216 = x191*x2
x217 = x17*x187
x218 = x217*x9
x219 = -x14*x218
x220 = x216 + x219
x221 = x220*x83
x222 = x217*x6
x223 = x14*x222
x224 = -x188*x4 + x223 + x23
x225 = x224*x81
x226 = x187*x6
x227 = x2*x226
x228 = x14*x227
x229 = x12*x15*x186*x9 - x228 - x24
x230 = x229*x82
x231 = x212*x28
x232 = x31*x47
x233 = -x210*x31 + x26*x37
x234 = x231 - x232 + x233
x235 = x234*x25
x236 = x28*x45
x237 = x26*x48
x238 = x210*x26
x239 = x238 - x46
x240 = x17*x206 - x204*x28 + x236 - x237 + x239
x241 = x240*x43
x242 = x212*x26
x243 = x31*x45
x244 = -x52
x245 = x243 + x244
x246 = x14*x21*x26*x31*x34 + x15*x186*x31*x44 - x2*x206 - x242 - x245
x247 = x246*x51
x248 = x221 + x225 + x230 + x235 + x241 + x247
x249 = x200*x85 + x215*x85 + x248*x67
x250 = 2*o0
x251 = 2*m0 - x250
x252 = x107*x29
x253 = x26*x28
x254 = x107*x253
x255 = x14*x254
x256 = -x48
x257 = x186*x34
x258 = x15*x257
x259 = x258*x27
x260 = x14*x257
x261 = x260*x27
x262 = x115*x261
x263 = -x113 + x210 + x256 - x259*x28 + x262 + x72
x264 = x26*x31
x265 = x107*x264
x266 = x17*x265
x267 = x17*x31
x268 = -x261*x267
x269 = x2*x254 + x268
x270 = x261*x28
x271 = x2*x270
x272 = -x39
x273 = x272 + x37
x274 = x271 + x273
x275 = x27*x35
x276 = x275*x31
x277 = -x276
x278 = -x212
x279 = -x251
x280 = x118*x261
x281 = x14*x265
x282 = x220*x40 + x224*x49 + x229*x55
x283 = 4*o0
x284 = 4*m0 - x283
x285 = x17*x254
x286 = x2*x265
x287 = x26*x275
x288 = 2*x287
x289 = 2*x26
x290 = x17*x270
x291 = -x201
x292 = x131 + x291
x293 = x129 + 2*x2*x261*x31 + x205 + x259*x289 - x288 + 2*x290 + x292 - 2*x70
x294 = x198*x78
x295 = x197*x294
x296 = (mpmath.mpf(1)/mpmath.mpf(4))*x84
x297 = x296*x77
x298 = x158*x200 + x158*x215 + x159*x248 + x160*x214 + x295*x297
x299 = 2*x23
x300 = x18*x187
x301 = x300*x4
x302 = x2*x218
x303 = -2*x223 + x299 - 2*x301 - 2*x302
x304 = x199*x63
x305 = x303*x304
x306 = x201*x28
x307 = x18*x203
x308 = 2*x307
x309 = x30*x35
x310 = x17*x309
x311 = x31*x37
x312 = x17*x203
x313 = x2*x312
x314 = x31*x313
x315 = -x143 + 2*x237 + 2*x238 + x28*x308 - x306 + 2*x310 + 2*x311 + 2*x314
x316 = x133*x315
x317 = x2*x222
x318 = -x219 - x317
x319 = x318*x82
x320 = x190*x2
x321 = x17*x320
x322 = x24 - x300*x9 + x321
x323 = x322*x83
x324 = x12*x18*x186*x6 - x192 - x42
x325 = x324*x81
x326 = x26*x313 - x31*x48
x327 = -x233 - x326
x328 = x327*x51
x329 = x128*x26
x330 = x210*x28
x331 = x26*x45
x332 = -x331
x333 = x332 + x58
x334 = -x14*x309 + x26*x307 + x329 - x330 + x333
x335 = x334*x43
x336 = x28*x313
x337 = x128*x31
x338 = x2*x309 + x245 - x307*x31 + x336 - x337
x339 = x25*x338
x340 = x319 + x323 + x325 + x328 + x335 + x339
x341 = x305*x85 + x316*x85 + x340*x67
x342 = x107*x30
x343 = x18*x257
x344 = x27*x343
x345 = x128 + x307
x346 = -x45
x347 = -x287 + x346 + x70
x348 = x26*x344 - x290 + x345 + x347
x349 = x28*x31
x350 = x107*x349
x351 = x17*x350
x352 = x257*x27
x353 = x2*x28
x354 = x17*x353
x355 = x352*x354
x356 = x276 - x74
x357 = -x31*x344 + x355 + x356
x358 = x2*x352
x359 = x115*x358 + x313
x360 = -x14*x350 + x359
x361 = x318*x55 + x322*x40 + x324*x49
x362 = x2*x350
x363 = 2*x48
x364 = 2*x72
x365 = 2*x28
x366 = x267*x358
x367 = -x114 + x211 + 2*x262 + x344*x365 + x363 - x364 + 2*x366
x368 = x294*x63
x369 = x297*x368
x370 = x158*x305 + x158*x316 + x159*x340 + x160*x315 + x303*x369
x371 = 2*x24
x372 = x187*x19
x373 = x372*x9
x374 = -2*x228 - 2*x321 + x371 - 2*x373
x375 = x304*x374
x376 = x201*x31
x377 = x19*x203
x378 = 2*x377
x379 = x32*x35
x380 = x2*x379
x381 = x119*x26
x382 = x129*x31 - x170 + 2*x242 + x31*x378 + 2*x336 - x376 + 2*x380 + 2*x381
x383 = x133*x382
x384 = -x216 + x317
x385 = x384*x81
x386 = x195 - x372*x6 + x42
x387 = x386*x82
x388 = x12*x186*x19*x4 - x23 - x302
x389 = x388*x83
x390 = -x231 + x232 + x326
x391 = x390*x43
x392 = x14*x186*x2*x31*x44 + x14*x21*x32*x34 - x208 - x26*x377 - x333
x393 = x392*x51
x394 = -x236 + x311
x395 = -x17*x379 + x28*x377 - x314 + x394 + x46
x396 = x25*x395
x397 = x385 + x387 + x389 + x391 + x393 + x396
x398 = x375*x85 + x383*x85 + x397*x67
x399 = -x271 + x39
x400 = x107*x32
x401 = -x72
x402 = x19*x257
x403 = x27*x402
x404 = x113 + x28*x403 - x366 + x401
x405 = x130 + x377
x406 = -x14*x186*x2*x27*x31*x34 + x26*x403 + x347 + x405
x407 = x384*x49 + x386*x55 + x388*x40
x408 = 2*x119
x409 = 2*x276
x410 = 2*x31
x411 = x213 + 2*x280 + 2*x355 + x403*x410 + x408 - x409 - x75
x412 = x158*x375 + x158*x383 + x159*x397 + x160*x382 + x369*x374
x413 = x13*x209
x414 = 2*x413
x415 = x13*x312
x416 = 2*x415
x417 = x2*x203
x418 = x13*x417
x419 = 2*x418
x420 = x202 + x26*x414 + x28*x416 + x31*x419
x421 = x133*x420
x422 = x14*x22
x423 = 2*x422
x424 = x11**(mpmath.mpf(-3)/mpmath.mpf(2))
x425 = x21*x424
x426 = x425*x7
x427 = x425*x6
x428 = x4*x427
x429 = x17*x428
x430 = x427*x9
x431 = x2*x430
x432 = x14*x226
x433 = 2*x432
x434 = x17*x190
x435 = 2*x434
x436 = x193*x2
x437 = 2*x436
x438 = -x13*x433 - x13*x435 - x13*x437 - 2*x14*x426 - x185 + x423 - 2*x429 - 2*x431
x439 = x304*x438
x440 = x28*x418 - x31*x415
x441 = x25*x440
x442 = x243 - x26*x418 + x31*x413
x443 = x442*x51
x444 = x13*x17*x186*x26*x44 - x236 - x28*x413
x445 = x43*x444
x446 = x2*x428
x447 = x17*x430
x448 = -x13*x218 + x13*x320 + x446 - x447
x449 = x448*x83
x450 = x14*x430
x451 = x2*x22
x452 = x24 + x451
x453 = x13*x194 - x13*x227 - x2*x426 + x450 + x452
x454 = x453*x82
x455 = x14*x428
x456 = x17*x22
x457 = x23 + x456
x458 = x12*x13*x17*x186*x6 - x13*x191 + x17*x21*x424*x7 - x455 - x457
x459 = x458*x81
x460 = x441 + x443 + x445 + x449 + x454 + x459
x461 = x421*x85 + x439*x85 + x460*x67
x462 = x13*x352
x463 = x261*x289
x464 = x124*x28
x465 = x358*x410
x466 = x134*(x13*x463 + x13*x465 + x201 + x288 + x414 + x462*x464) + x158*x421 + x158*x439 + x159*x460 + x160*x420 + x369*x438 + x67*(x25*(-x267*x462 + x353*x462) + x40*x448 + x43*(x115*x462 - x13*x270 + x401 + x415) + x453*x55 + x458*x49 + x51*(-x118*x462 + x13*x14*x186*x27*x31*x34 - x277 - x418))
x467 = x16*x209
x468 = 2*x467
x469 = x16*x312
x470 = 2*x469
x471 = x16*x417
x472 = 2*x471
x473 = x26*x468 + x28*x470 + x306 + x31*x472
x474 = x133*x473
x475 = 2*x456
x476 = x425*x8
x477 = x4*x425
x478 = x477*x9
x479 = x2*x478
x480 = -x16*x433 - x16*x435 - x16*x437 - 2*x17*x476 - x299 - 2*x455 + x475 - 2*x479
x481 = x304*x480
x482 = x14*x16*x186*x31*x44 - x26*x471
x483 = x482*x51
x484 = x26*x469 - x28*x467 + x331
x485 = x43*x484
x486 = x16*x186*x2*x28*x44 - x243 - x31*x469
x487 = x25*x486
x488 = -x14*x478
x489 = x12*x14*x16*x186*x9 - x16*x227 - x446 - x488
x490 = x489*x82
x491 = x42 + x422
x492 = -x14*x476 - x16*x191 + x16*x222 + x429 + x491
x493 = x492*x81
x494 = x17*x478
x495 = x12*x16*x186*x2*x4 - x16*x218 + x2*x21*x424*x8 - x452 - x494
x496 = x495*x83
x497 = x483 + x485 + x487 + x490 + x493 + x496
x498 = x474*x85 + x481*x85 + x497*x67
x499 = x16*x352
x500 = x287 + x45
x501 = x134*(x16*x463 + x16*x465 + x364 + x464*x499 + x468) + x158*x474 + x158*x481 + x159*x497 + x160*x473 + x369*x480 + x67*(x25*(x16*x186*x2*x27*x28*x34 - x267*x499 - x276) + x40*x495 + x43*(x115*x499 - x16*x270 + x469 + x500) + x489*x55 + x49*x492 + x51*(-x118*x499 + x14*x16*x186*x27*x31*x34 - x471))
x502 = x1*x209
x503 = 2*x502
x504 = x1*x312
x505 = 2*x504
x506 = x1*x417
x507 = 2*x506
x508 = x26*x503 + x28*x505 + x31*x507 + x376
x509 = x133*x508
x510 = 2*x451
x511 = x10*x425
x512 = -x1*x433 - x1*x435 - x1*x437 - 2*x2*x511 - x371 - 2*x450 - 2*x494 + x510
x513 = x304*x512
x514 = x26*x504 - x28*x502
x515 = x43*x514
x516 = x236 + x28*x506 - x31*x504
x517 = x25*x516
x518 = x1*x14*x186*x31*x44 - x26*x506 - x331
x519 = x51*x518
x520 = -x1*x191 + x1*x222 + x447 + x488
x521 = x520*x81
x522 = -x1*x218 + x1*x320 - x17*x511 + x457 + x479
x523 = x522*x83
x524 = x1*x12*x14*x186*x9 - x1*x227 + x10*x14*x21*x424 - x431 - x491
x525 = x524*x82
x526 = x515 + x517 + x519 + x521 + x523 + x525
x527 = x509*x85 + x513*x85 + x526*x67
x528 = x1*x352
x529 = x134*(x1*x463 + x1*x465 + x409 + x464*x528 + x503) + x158*x509 + x158*x513 + x159*x526 + x160*x508 + x369*x512 + x67*(x25*(-x267*x528 + x353*x528 + x72) + x40*x522 + x43*(-x1*x270 + x115*x528 + x504) + x49*x520 + x51*(x1*x14*x186*x27*x31*x34 - x118*x528 - x500 - x506) + x524*x55)
x530 = -x6
x531 = x477*x530
x532 = x425*x9
x533 = x530*x532
x534 = -x17*x533 + x2*x531
x535 = x534*x83
x536 = x427*x530
x537 = -x14*x531 + x17*x536 + x456
x538 = x537*x81
x539 = x14*x21*x424*x530*x9 - x2*x536 - x451
x540 = x539*x82
x541 = x535 + x538 + x540
x542 = 2*x536
x543 = 2*x531
x544 = 2*x533
x545 = -x14*x542 - x17*x543 - x2*x544 - x423
x546 = x304*x545
x547 = x541*x67 + x546*x85
x548 = x97*(x158*x546 + x159*x541 + x369*x545 + x67*(x40*x534 + x49*x537 + x539*x55))
x549 = -x4
x550 = x427*x549
x551 = x14*x21*x424*x549*x9 - x2*x550
x552 = x551*x82
x553 = x477*x549
x554 = x532*x549
x555 = -x17*x554 + x2*x553 + x451
x556 = x555*x83
x557 = -x14*x553 + x17*x21*x424*x549*x6 - x422
x558 = x557*x81
x559 = x552 + x556 + x558
x560 = 2*x550
x561 = 2*x553
x562 = 2*x554
x563 = -x14*x560 - x17*x561 - x2*x562 - x475
x564 = x304*x563
x565 = x559*x67 + x564*x85
x566 = x97*(x158*x564 + x159*x559 + x369*x563 + x67*(x40*x555 + x49*x557 + x55*x551))
x567 = -x9
x568 = x427*x567
x569 = x477*x567
x570 = -x14*x569 + x17*x568
x571 = x570*x81
x572 = x532*x567
x573 = x14*x572 - x2*x568 + x422
x574 = x573*x82
x575 = -x17*x572 + x2*x21*x4*x424*x567 - x456
x576 = x575*x83
x577 = x571 + x574 + x576
x578 = 2*x568
x579 = 2*x569
x580 = 2*x572
x581 = -x14*x578 - x17*x579 - x2*x580 - x510
x582 = x304*x581
x583 = x577*x67 + x582*x85
x584 = x97*(x158*x582 + x159*x577 + x369*x581 + x67*(x40*x575 + x49*x570 + x55*x573))
x585 = x59*x94
x586 = x92*(x148*x585 + x95*(2*x137 + 2*x139 + 2*x141))
x587 = 3*o1
x588 = 3*m1 - x587
x589 = -x588
x590 = x106*x589
x591 = x27*x590
x592 = x125*x31
x593 = x100*(x128*x135 + x144*x26 + x147 + x46)
x594 = x135*x590
x595 = x256 + x47
x596 = 2*x146
x597 = x162*x590
x598 = x14*x597
x599 = -x154 + x173
x600 = 2*x178
x601 = x596 + x600
x602 = (mpmath.mpf(1)/mpmath.mpf(2))*x149
x603 = (mpmath.mpf(1)/mpmath.mpf(2))*x142
x604 = x148*x99
x605 = x169*x602 + x175*x604 + x176*x603
x606 = x253*x590
x607 = x135*x258
x608 = x126 + x204
x609 = x135*x36
x610 = x154 + x346 - x609
x611 = -x135*x14*x17*x186*x26*x34 + x28*x607 + x608 + x610
x612 = x135*x260
x613 = -x267*x612
x614 = x2*x606 + x613
x615 = x264*x590
x616 = x212 + x353*x612
x617 = -x17*x615 + x616
x618 = x29*x590
x619 = x118*x612
x620 = x135*x35
x621 = x31*x620
x622 = -x146 + x621
x623 = -x135*x15*x186*x31*x34 + x619 + x622
x624 = x136*x229 + x138*x220 + x140*x224
x625 = 2*x47
x626 = x26*x620
x627 = 2*x626
x628 = x28*x612
x629 = x31*x612
x630 = x124*x628 + x125*x629 - x145 + x211 + x289*x607 + x625 - x627
x631 = x296*x59
x632 = x148*x631
x633 = x200*x603 + x214*x604 + x215*x603 + x248*x602 + x295*x632
x634 = 2*o1
x635 = 2*m1 - x634
x636 = -x635
x637 = x135*x343
x638 = x30*x590
x639 = x17*x628
x640 = x144 - x626
x641 = x349*x590
x642 = x135*x257
x643 = x354*x642
x644 = x272 - x31*x637 + x313 + x622 + x643
x645 = x14*x641
x646 = x2*x642
x647 = x115*x646
x648 = -x119
x649 = x54 + x648
x650 = x647 + x649
x651 = x136*x318 + x138*x322 + x140*x324
x652 = 4*o1
x653 = 4*m1 - x652
x654 = 2*x609
x655 = x125*x267
x656 = 2*x115*x612 + x127 - 2*x154 + x292 + x308 + x365*x637 + x642*x655 - x654
x657 = x368*x632
x658 = x303*x657 + x305*x603 + x315*x604 + x316*x603 + x340*x602
x659 = x119 - x647
x660 = x32*x590
x661 = x135*x402
x662 = -x267*x646 + x28*x661 + x405 + x610
x663 = -x135*x14*x186*x2*x31*x34 + x26*x661 + x640
x664 = x136*x386 + x138*x388 + x140*x384
x665 = 2*x39
x666 = 2*x621
x667 = 2*x313
x668 = x410*x661 - x596 + 2*x619 + 2*x643 + x665 - x666 + x667
x669 = x374*x657 + x375*x603 + x382*x604 + x383*x603 + x397*x602
x670 = x13*x642
x671 = x45 + x609
x672 = x289*x612
x673 = x134*(x13*x672 + x416 + x464*x670 + x592*x670 + x627) + x420*x604 + x421*x603 + x438*x657 + x439*x603 + x460*x602 + x67*(x136*x453 + x138*x448 + x140*x458 + x25*(-x267*x670 + x353*x670 + x418) + x43*(x13*x135*x17*x186*x26*x34 - x13*x628 - x413 - x671) + x51*(-x118*x670 + x13*x629 + x621))
x674 = x16*x642
x675 = x134*(x16*x672 + x201 + x464*x674 + x470 + x592*x674 + x654) + x473*x604 + x474*x603 + x480*x657 + x481*x603 + x497*x602 + x67*(x136*x489 + x138*x495 + x140*x492 + x25*(-x267*x674 + x353*x674 + x471 - x621) + x43*(x115*x674 - x16*x628 - x467 + x626) + x51*(-x118*x674 + x135*x14*x16*x186*x31*x34))
x676 = x1*x642
x677 = x134*(x1*x672 + x464*x676 + x505 + x592*x676 + x666) + x508*x604 + x509*x603 + x512*x657 + x513*x603 + x526*x602 + x67*(x136*x524 + x138*x522 + x140*x520 + x25*(-x267*x676 + x353*x676 + x506 + x671) + x43*(x1*x135*x17*x186*x26*x34 - x1*x628 - x502) + x51*(x1*x135*x14*x186*x31*x34 - x118*x676 - x626))
x678 = x97*(x541*x602 + x545*x657 + x546*x603 + x67*(x136*x539 + x138*x534 + x140*x537))
x679 = x97*(x559*x602 + x563*x657 + x564*x603 + x67*(x136*x551 + x138*x555 + x140*x557))
x680 = x97*(x577*x602 + x581*x657 + x582*x603 + x67*(x136*x573 + x138*x575 + x140*x570))
x681 = x92*(x175*x585 + x95*(2*x164 + 2*x166 + 2*x168))
x682 = 3*o2
x683 = 3*m2 - x682
x684 = -x683
x685 = x106*x684
x686 = x27*x685
x687 = x100*(x128*x162 + x171*x26 + x174 + x52)
x688 = x135*x685
x689 = x14*x688
x690 = x162*x685
x691 = x14*x690
x692 = x29*x685
x693 = x253*x685
x694 = x162*x258
x695 = x162*x260
x696 = x115*x695
x697 = x162*x36
x698 = -x178 + x697
x699 = -x28*x694 + x696 + x698
x700 = x264*x685
x701 = -x353*x695
x702 = x17*x700 + x701
x703 = x210 + x267*x695
x704 = -x2*x693 + x703
x705 = x118*x695
x706 = x162*x35
x707 = x31*x706
x708 = x173 + x346 - x707
x709 = x31*x694 + x608 - x705 + x708
x710 = x163*x224 + x165*x229 + x167*x220
x711 = 2*x54
x712 = x26*x706
x713 = 2*x712
x714 = x28*x695
x715 = x31*x695
x716 = x124*x714 + x125*x715 - x172 + x213 + x289*x694 + x711 - x713
x717 = (mpmath.mpf(1)/mpmath.mpf(2))*x169
x718 = (mpmath.mpf(1)/mpmath.mpf(2))*x176
x719 = x175*x99
x720 = x175*x631
x721 = x200*x717 + x214*x719 + x215*x717 + x248*x718 + x295*x720
x722 = x30*x685
x723 = x162*x343
x724 = -x17*x714 + x171 + x26*x723 - x712
x725 = x349*x685
x726 = -x162*x17*x186*x2*x28*x34 + x31*x723 + x345 + x708
x727 = x14*x725
x728 = x162*x257
x729 = x2*x728
x730 = x115*x729
x731 = x48 - x730
x732 = x163*x324 + x165*x318 + x167*x322
x733 = 2*x37
x734 = 2*x697
x735 = x365*x723 - x600 + x655*x728 + x667 + 2*x696 + x733 - x734
x736 = x368*x720
x737 = x303*x736 + x305*x717 + x315*x719 + x316*x717 + x340*x718
x738 = x595 + x730
x739 = -x313
x740 = 2*o2
x741 = 2*m2 - x740
x742 = -x741
x743 = x162*x402
x744 = x32*x685
x745 = x267*x729
x746 = x120 - x171 + x2*x715 + x212 - x26*x743 + x712
x747 = x163*x384 + x165*x386 + x167*x388
x748 = 4*o2
x749 = 4*m2 - x748
x750 = 2*x707
x751 = x127 + x129 - 2*x173 + x291 + 2*x354*x728 + x378 + x410*x743 + 2*x705 - x750
x752 = x374*x736 + x375*x717 + x382*x719 + x383*x717 + x397*x718
x753 = x13*x728
x754 = x45 + x707
x755 = x289*x695
x756 = x134*(x13*x755 + x419 + x464*x753 + x592*x753 + x713) + x420*x719 + x421*x717 + x438*x736 + x439*x717 + x460*x718 + x67*(x163*x458 + x165*x453 + x167*x448 + x25*(x13*x162*x186*x2*x28*x34 - x267*x753 - x415) + x43*(x13*x162*x17*x186*x26*x34 - x13*x714 - x697) + x51*(-x118*x753 + x13*x715 + x413 + x754))
x757 = x16*x728
x758 = x134*(x16*x755 + x464*x757 + x472 + x592*x757 + x734) + x473*x719 + x474*x717 + x480*x736 + x481*x717 + x497*x718 + x67*(x163*x492 + x165*x489 + x167*x495 + x25*(x16*x162*x186*x2*x28*x34 - x267*x757 - x469 - x754) + x43*(x115*x757 - x16*x714 + x712) + x51*(-x118*x757 + x16*x715 + x467))
x759 = x1*x728
x760 = -x502
x761 = -x504
x762 = x134*(x1*x755 + x201 + x464*x759 + x507 + x592*x759 + x750) + x508*x719 + x509*x717 + x512*x736 + x513*x717 + x526*x718 + x67*(x163*x520 + x165*x524 + x167*x522 + x25*(-x267*x759 + x353*x759 + x697 + x761) + x43*(-x1*x714 + x115*x759) + x51*(x1*x14*x162*x186*x31*x34 - x118*x759 - x712 - x760))
x763 = x97*(x541*x718 + x545*x736 + x546*x717 + x67*(x163*x537 + x165*x539 + x167*x534))
x764 = x97*(x559*x718 + x563*x736 + x564*x717 + x67*(x163*x557 + x165*x551 + x167*x555))
x765 = x97*(x577*x718 + x581*x736 + x582*x717 + x67*(x163*x570 + x165*x573 + x167*x575))
x766 = x88*x93/x65**2
x767 = x92*(x197*x766 + x214*x585 + x95*(2*x221 + 2*x225 + 2*x230 + 2*x235 + 2*x241 + 2*x247))
x768 = x103*x106
x769 = x27*x768
x770 = x100*(x204*x26 + x207 + x208 + x212*x31 + x329 + x330 + x332 - x58)
x771 = x135*x768
x772 = x162*x768
x773 = x14*x772
x774 = x62*x64*x98/x65**(mpmath.mpf(5)/mpmath.mpf(2))
x775 = x199*x85
x776 = x775*(-x189 - x192 - x195 + x42)
x777 = -4*p0
x778 = x283 + x777
x779 = -x778
x780 = x20**(mpmath.mpf(-5)/mpmath.mpf(2))
x781 = -3*p0
x782 = x102 + x781
x783 = -x782
x784 = x780*x783
x785 = x15*x784
x786 = x12*x6
x787 = 2*x786
x788 = x12*x4
x789 = x17*x788
x790 = x122*x784
x791 = x12*x9
x792 = x2*x784
x793 = x791*x792
x794 = x435 + x437
x795 = x433 + x794
x796 = x304*x85
x797 = x59*x85
x798 = x203*x26
x799 = x257*x29
x800 = x29*x768
x801 = x44*x785
x802 = x253*x768
x803 = x264*x768
x804 = x28*x44
x805 = x17*x804
x806 = x44*x792
x807 = x260*x264
x808 = 4*x2*x807
x809 = 4*x45
x810 = x253*x260
x811 = x209*x26
x812 = 2*x811
x813 = x28*x312
x814 = 2*x813
x815 = x31*x417
x816 = 2*x815
x817 = -x812 - x814 - x816
x818 = x132 + 4*x17*x810 + x809 + x817
x819 = -x320
x820 = x14*x784
x821 = x2*x820
x822 = x788*x821
x823 = x17*x820
x824 = x218 - x791*x823
x825 = -2*p0
x826 = x250 + x825
x827 = -x826
x828 = x17*x786
x829 = -x222
x830 = x191 + x829
x831 = -x194 + x227
x832 = x44*x820
x833 = x353*x832
x834 = x17*x803
x835 = x124*x807
x836 = x125*x810
x837 = x2*x802 - x267*x832
x838 = x31*x312
x839 = x28*x417
x840 = -x839
x841 = x838 + x840
x842 = x112 + x841
x843 = x26*x36
x844 = 2*x843
x845 = x203*x28
x846 = 2*x258
x847 = x253*x846
x848 = x17*x799
x849 = x122*x848
x850 = x209*x28
x851 = x26*x312
x852 = -x851
x853 = x850 + x852
x854 = x117 + x853
x855 = x203*x31
x856 = x264*x35
x857 = 2*x856
x858 = x2*x799
x859 = x122*x858
x860 = x264*x846
x861 = x26*x417
x862 = x209*x31
x863 = -x862
x864 = x861 + x863
x865 = x649 + x864
x866 = x18*x784
x867 = -x432
x868 = x434 + x867
x869 = x792*x828
x870 = x44*x866
x871 = x30*x768
x872 = x309 + x813
x873 = -x206 + x608 - x811
x874 = -x128 + x18*x799 - x258*x30 - x307 + x872 + x873
x875 = x44*x784
x876 = x354*x875
x877 = x17*x349
x878 = x2*x30
x879 = x260*x878
x880 = x264*x343
x881 = x2*x253
x882 = x17*x257
x883 = -x260*x877 + x881*x882
x884 = x120 + x278 + x856 + x862 + x879 - x880 + x883
x885 = x349*x768
x886 = x31*x36
x887 = x2*x848 - x258*x349 + x739 + x886
x888 = x115*x806 - x14*x885 + x887
x889 = x260*x881
x890 = x111 - x14*x17*x186*x26*x31*x34 + x838 + x889
x891 = x220*x338 + x224*x334 + x229*x327 + x234*x322 + x240*x324 + x246*x318
x892 = x122*x875
x893 = x2*x264
x894 = x124*x257*x893
x895 = 2*x343
x896 = x253*x895
x897 = 2*x851
x898 = 2*x850
x899 = x260*x30
x900 = x124*x899
x901 = x260*x349
x902 = x125*x901
x903 = -4*x210 - x363 - x625 - 4*x843 + x847 + x849 + x894 + x896 - x897 - x898 + x900 + x902
x904 = 2*x788
x905 = 2*x222
x906 = -x905
x907 = 2*x191
x908 = -x907
x909 = x906 + x908
x910 = (mpmath.mpf(1)/mpmath.mpf(2))*x248
x911 = (mpmath.mpf(1)/mpmath.mpf(2))*x340
x912 = x214*x99
x913 = x196*x774
x914 = x295*x631
x915 = x303*x368
x916 = x214*x631
x917 = x200*x911 + x215*x911 + x303*x913 + x305*x910 + x315*x912 + x315*x914 + x316*x910 + x915*x916
x918 = x17*x791
x919 = x19*x784
x920 = x436 + x867
x921 = x17*x807 + x272 + x839 - x889
x922 = x44*x919
x923 = x17*x32
x924 = -x850
x925 = x2*x349
x926 = x264*x882
x927 = x2*x926
x928 = -x843
x929 = x210 + x253*x402 - x260*x923 + x928
x930 = x260*x925 + x47 + x924 - x927 + x929
x931 = x379 + x815
x932 = -x130 - x377 + x931
x933 = -x15*x186*x32*x34 + x19*x799 + x873 + x932
x934 = x220*x395 + x224*x390 + x229*x392 + x234*x388 + x240*x384 + x246*x386
x935 = x125*x32
x936 = x257*x881
x937 = x124*x936
x938 = x260*x32
x939 = x125*x938
x940 = 2*x861
x941 = 2*x862
x942 = 2*x402
x943 = x264*x942
x944 = x124*x901 - 4*x212 - x408 - x711 - 4*x856 + x859 + x860 + x937 + x939 - x940 - x941 + x943
x945 = 2*x791
x946 = x786*x792
x947 = x124*x788
x948 = 2*x227
x949 = -x948
x950 = 2*x194
x951 = -x950
x952 = x949 + x951
x953 = (mpmath.mpf(1)/mpmath.mpf(2))*x397
x954 = x368*x916
x955 = x200*x953 + x215*x953 + x374*x913 + x374*x954 + x375*x910 + x382*x912 + x382*x914 + x383*x910
x956 = x13*x784
x957 = x13*x788
x958 = x13*x190
x959 = x186*x424
x960 = x14*x959
x961 = x17*x7
x962 = x960*x961
x963 = x15*x959
x964 = x6*x963
x965 = x4*x964
x966 = x14*x217
x967 = x428 - x966
x968 = x222 + x958 + x962 - x965 + x967
x969 = -x218
x970 = x17*x960
x971 = x9*x970
x972 = x6*x971
x973 = x2*x960
x974 = x4*x6
x975 = x973*x974
x976 = x320 + x975
x977 = x969 - x972 + x976
x978 = x13*x193
x979 = x7*x973
x980 = x14*x187
x981 = x2*x980
x982 = x430 - x981
x983 = -x15*x186*x424*x6*x9 + x227 + x978 + x979 + x982
x984 = x13*x810
x985 = x13*x832
x986 = x13*x845
x987 = x415 - x986
x988 = x13*x875
x989 = x13*x936
x990 = -x13*x926
x991 = -x838
x992 = x839 + x991
x993 = x989 + x990 + x992
x994 = -x861
x995 = x13*x855
x996 = x13*x807
x997 = -x13*x858 + x418 + x856 + x994 - x995 + x996
x998 = x220*x440 + x224*x444 + x229*x442 + x234*x448 + x240*x458 + x246*x453
x999 = x124*x804
x1000 = x13*x798
x1001 = x122*x799
x1002 = x13*x257
x1003 = x124*x253
x1004 = x125*x264
x1005 = x291 + x816
x1006 = -2*x1000 + x1001*x13 + x1002*x1003 + x1002*x1004 + x1005 + 2*x206 - x414 + x814
x1007 = x122*x786
x1008 = x125*x791
x1009 = 2*x188
x1010 = 2*x7
x1011 = x9*x973
x1012 = x1011*x6
x1013 = 2*x1012
x1014 = 2*x22
x1015 = x970*x974
x1016 = x1014 + 2*x1015
x1017 = -x1009 + x1010*x963 + x1013 + x1016 - 2*x12*x13*x186*x6 - 2*x21*x424*x7 + x794
x1018 = (mpmath.mpf(1)/mpmath.mpf(2))*x460
x1019 = x1018*x200 + x1018*x215 + x420*x912 + x420*x914 + x421*x910 + x438*x913 + x438*x954 + x439*x910
x1020 = x16*x784
x1021 = x16*x788
x1022 = -x22
x1023 = x1015 + x1022
x1024 = x16*x190
x1025 = x1024 + x476
x1026 = x188 + x432
x1027 = x1023 + x1025 + x1026 - x8*x963
x1028 = x4*x971
x1029 = -x8*x973 + x981
x1030 = x1028 + x1029 + x194
x1031 = x16*x193
x1032 = x4*x963
x1033 = -x1032*x9 + x478
x1034 = x1031 + x1033 + x975
x1035 = x16*x875
x1036 = x16*x810
x1037 = x206 + x346 + x811
x1038 = x16*x845
x1039 = x1038 - x469
x1040 = -x1036 + x1037 + x1039 + x16*x848
x1041 = x16*x936
x1042 = x16*x926 + x856
x1043 = -x1041 + x1042 + x862
x1044 = x16*x855
x1045 = x1044 - x471
x1046 = x1045 - x14*x16*x186*x26*x31*x34 + x16*x858
x1047 = x220*x486 + x224*x484 + x229*x482 + x234*x495 + x240*x492 + x246*x489
x1048 = x16*x798
x1049 = 2*x1048
x1050 = x16*x257
x1051 = x1001*x16 + x1003*x1050 + x1004*x1050 - x1049 - x468 + x844 + x898
x1052 = x16*x226
x1053 = 2*x1052
x1054 = x8*x970
x1055 = x1011*x4
x1056 = 2*x428
x1057 = 2*x966
x1058 = -x1056 - x1057
x1059 = -x1053 + 2*x1054 + 2*x1055 + x1058 + x907 + 2*x965
x1060 = (mpmath.mpf(1)/mpmath.mpf(2))*x497
x1061 = x1060*x200 + x1060*x215 + x473*x912 + x473*x914 + x474*x910 + x480*x913 + x480*x954 + x481*x910
x1062 = x1*x784
x1063 = x1*x788
x1064 = x1*x190
x1065 = x1033 + x1064 + x972
x1066 = -x10*x970 + x966
x1067 = x1055 + x1066 + x191
x1068 = x1012 + x1022
x1069 = x1*x193
x1070 = x1069 + x511
x1071 = -x10*x15*x186*x424 + x1026 + x1068 + x1070
x1072 = x1*x875
x1073 = x1*x845
x1074 = x1*x810
x1075 = x1*x848 + x1073 - x1074 + x761
x1076 = x1*x926
x1077 = x1*x936 + x843
x1078 = -x1076 + x1077 + x850
x1079 = x1*x855
x1080 = x1079 - x506
x1081 = -x1*x14*x186*x26*x31*x34 + x1*x858 + x1037 + x1080
x1082 = x220*x516 + x224*x514 + x229*x518 + x234*x522 + x240*x520 + x246*x524
x1083 = x1*x798
x1084 = 2*x1083
x1085 = x1*x257
x1086 = x1*x1001 + x1003*x1085 + x1004*x1085 - x1084 - x503 + x857 + x941
x1087 = x1*x226
x1088 = 2*x1087
x1089 = x10*x973
x1090 = 2*x430
x1091 = 2*x981
x1092 = -x1090 - x1091
x1093 = 2*x1028 - x1088 + 2*x1089 + x1092 + 2*x9*x964 + x950
x1094 = (mpmath.mpf(1)/mpmath.mpf(2))*x526
x1095 = x1094*x200 + x1094*x215 + x508*x912 + x508*x914 + x509*x910 + x512*x913 + x512*x954 + x513*x910
x1096 = x4*x973
x1097 = x1096*x530
x1098 = -x530*x971
x1099 = x530*x6
x1100 = x1099*x970
x1101 = x1099*x973
x1102 = (mpmath.mpf(1)/mpmath.mpf(2))*x541
x1103 = -x1014
x1104 = 2*x964
x1105 = x4*x530
x1106 = x1105*x970
x1107 = x1011*x530
x1108 = x1102*x200 + x1102*x215 + x545*x913 + x545*x954 + x546*x910 + x67*(x234*x534 + x240*x537 + x246*x539 + x81*(-x1032*x530 + x1100 + x531 + x966) + x82*(-x1101 + x15*x186*x424*x530*x9 - x533 - x981) + x83*(x1097 + x1098)) + x796*(-x1009 - x1103 - x1104*x530 - 2*x1106 - 2*x1107 + x542)
x1109 = x549*x6
x1110 = x1109*x973
x1111 = -x549*x971
x1112 = x1096*x549 + x981
x1113 = x1109*x970
x1114 = -x188 + x22
x1115 = (mpmath.mpf(1)/mpmath.mpf(2))*x559
x1116 = x4*x549
x1117 = x1116*x970
x1118 = x1011*x549
x1119 = x1115*x200 + x1115*x215 + x563*x913 + x563*x954 + x564*x910 + x67*(x234*x555 + x240*x557 + x246*x551 + x81*(-x1032*x549 + x1113 + x1114 + x553) + x82*(-x1110 + x15*x186*x424*x549*x9 - x554) + x83*(x1111 + x1112)) + x796*(-x1057 - x1104*x549 - 2*x1117 - 2*x1118 + x560)
x1120 = x567*x6
x1121 = x1120*x970
x1122 = -x1096*x567
x1123 = x567*x971 + x966
x1124 = x1120*x973
x1125 = (mpmath.mpf(1)/mpmath.mpf(2))*x577
x1126 = x4*x567
x1127 = x1126*x970
x1128 = x1011*x567
x1129 = x1125*x200 + x1125*x215 + x581*x913 + x581*x954 + x582*x910 + x67*(x234*x575 + x240*x570 + x246*x573 + x81*(-x1032*x567 + x1121 + x569) + x82*(-x1114 - x1124 + x15*x186*x424*x567*x9 - x572) + x83*(-x1122 - x1123)) + x796*(-x1091 - x1104*x567 - 2*x1127 - 2*x1128 + x578)
x1130 = x63*x766
x1131 = x92*(x1130*x303 + x315*x585 + x95*(2*x319 + 2*x323 + 2*x325 + 2*x328 + 2*x335 + 2*x339))
x1132 = x106*x588
x1133 = x1132*x27
x1134 = x100*(x237 + x239 + x28*x307 + x310 + x314 + x394)
x1135 = x1132*x135
x1136 = x1132*x162
x1137 = x1136*x14
x1138 = -3*p1
x1139 = x1138 + x587
x1140 = -x1139
x1141 = x1140*x780
x1142 = x1141*x15
x1143 = x1141*x14
x1144 = x1143*x2
x1145 = x1144*x788
x1146 = -x1143*x918 + x194
x1147 = x1132*x29
x1148 = x1142*x44
x1149 = x1143*x44
x1150 = x1132*x14
x1151 = x1149*x353
x1152 = x1132*x17
x1153 = x1152*x264
x1154 = x1132*x881 - x1149*x267
x1155 = x1141*x122
x1156 = x31*x44
x1157 = x1141*x2
x1158 = x1157*x122
x1159 = x775*(-x223 + x23 - x301 - x302)
x1160 = -4*p1
x1161 = x1160 + x652
x1162 = -x1161
x1163 = x1141*x18
x1164 = x124*x791
x1165 = x124*x30
x1166 = x1163*x44
x1167 = x1132*x122
x1168 = x1132*x349
x1169 = x1141*x44
x1170 = x1169*x122
x1171 = 4*x882*x925
x1172 = -x227
x1173 = x1157*x828
x1174 = -2*p1
x1175 = x1174 + x634
x1176 = -x1175
x1177 = x218 + x819
x1178 = x862 + x994
x1179 = x1157*x44
x1180 = x115*x1179 - x1150*x349
x1181 = x14*x30
x1182 = x1143*x804
x1183 = 2*x886
x1184 = x349*x895
x1185 = x257*x878
x1186 = x1185*x124
x1187 = x1169*x354
x1188 = x1141*x19
x1189 = -x434 + x436
x1190 = -x856
x1191 = x119 + x1190 + x212 - x879 + x880 + x883 + x994
x1192 = x1188*x44
x1193 = x30*x402 - x309 - x32*x343 + x345 - x813 + x932
x1194 = -x14*x186*x2*x28*x31*x34 + x48 + x852 + x927 + x929
x1195 = x318*x392 + x322*x395 + x324*x390 + x327*x386 + x334*x384 + x338*x388
x1196 = x349*x942
x1197 = 2*x839
x1198 = 2*x838
x1199 = x257*x923
x1200 = x1199*x125
x1201 = x1184 + x1186 + x1196 - x1197 - x1198 + x1200 - 4*x313 - x665 - x733 + x835 + x836 - 4*x886
x1202 = x1157*x786
x1203 = 2*x320
x1204 = -x1203
x1205 = 2*x218
x1206 = -x1205
x1207 = x1204 + x1206
x1208 = x315*x99
x1209 = x303*x774
x1210 = x631*x915
x1211 = x368*x374
x1212 = x315*x631
x1213 = x1208*x382 + x1209*x374 + x1210*x382 + x1211*x1212 + x305*x953 + x316*x953 + x375*x911 + x383*x911
x1214 = x300 + x434
x1215 = x13*x226
x1216 = x1215 + x426
x1217 = x1023 + x1214 + x1216 - x18*x186*x424*x7
x1218 = x1141*x13
x1219 = x2*x959
x1220 = x1219*x17
x1221 = x1220*x974
x1222 = x18*x959
x1223 = x1222*x6
x1224 = -x1223*x9 + x430
x1225 = x1221 + x1224 + x978
x1226 = x1143*x13
x1227 = x2*x217
x1228 = -x1219*x961 + x1227
x1229 = x1228 + x218 + x972
x1230 = x346 + x872
x1231 = x1000 - x413
x1232 = x1230 + x1231 - x13*x17*x186*x26*x28*x34 + x13*x899
x1233 = x1169*x13
x1234 = x1185*x13 - x13*x349*x882 - x418 + x995
x1235 = x13*x901 + x886
x1236 = x1235 + x838 - x989
x1237 = x318*x442 + x322*x440 + x324*x444 + x327*x453 + x334*x458 + x338*x448
x1238 = x1156*x125
x1239 = 2*x986
x1240 = x1165*x257
x1241 = x125*x349
x1242 = x1002*x1241 - x1239 + x1240*x13 - x416 + x844 + x897 + 2*x984
x1243 = x1141*x124
x1244 = 2*x958
x1245 = x1220*x9
x1246 = x1245*x6
x1247 = x1058 + 2*x1223*x4 - x1244 + 2*x1246 + x905 + 2*x962
x1248 = x1212*x368
x1249 = x1018*x305 + x1018*x316 + x1208*x420 + x1209*x438 + x1210*x420 + x1248*x438 + x421*x911 + x439*x911
x1250 = x1052 + x1054 - x18*x186*x4*x424*x6 + x191 + x967
x1251 = x1141*x16
x1252 = x1220*x8
x1253 = x1222*x9
x1254 = x1253*x4
x1255 = -x1227 + x478
x1256 = x1031 + x1252 - x1254 + x1255 + x320
x1257 = -x1028 + x1221 + x831
x1258 = x1169*x16
x1259 = x16*x882
x1260 = -x1048 + x467
x1261 = x1259*x253 + x1260 - x16*x899 + x843 + x924
x1262 = x1259*x349
x1263 = -x16*x901
x1264 = x1041 + x1263 + x864
x1265 = x318*x482 + x322*x486 + x324*x484 + x327*x489 + x334*x492 + x338*x495
x1266 = x1005 + 2*x1036 - 2*x1038 + x1050*x1241 + x1240*x16 + 2*x309 - x470 + x812
x1267 = 2*x300
x1268 = 2*x1222
x1269 = x1245*x4
x1270 = 2*x1269 + x433
x1271 = x1016 - 2*x12*x16*x186*x4 - x1267 + x1268*x8 + x1270 - 2*x21*x424*x8 + x437
x1272 = x1060*x305 + x1060*x316 + x1208*x473 + x1209*x480 + x1210*x473 + x1248*x480 + x474*x911 + x481*x911
x1273 = x1028 + x1087 + x1224
x1274 = x1*x1141
x1275 = x1022 + x1269
x1276 = -x10*x1222 + x1070 + x1214 + x1275
x1277 = x1246 + x222
x1278 = x1066 + x1277
x1279 = x1*x1169
x1280 = x1*x882
x1281 = -x1083 + x502
x1282 = -x1*x899 + x1280*x253 + x1281
x1283 = x1*x1185 + x1080 + x1230 - x1280*x349
x1284 = -x1*x901
x1285 = x1077 + x1284 + x851
x1286 = x318*x518 + x322*x516 + x324*x514 + x327*x524 + x334*x520 + x338*x522
x1287 = 2*x1073
x1288 = x1*x1240 + 2*x1074 + x1085*x1241 + x1183 + x1198 - x1287 - x505
x1289 = 2*x1064
x1290 = x10*x1220
x1291 = 2*x478
x1292 = 2*x1227
x1293 = -x1291 - x1292
x1294 = x1205 + 2*x1254 - x1289 + 2*x1290 + x1293 + 2*x972
x1295 = x1094*x305 + x1094*x316 + x1208*x508 + x1209*x512 + x1210*x508 + x1248*x512 + x509*x911 + x513*x911
x1296 = x1105*x1220
x1297 = x1099*x1220 + x1227
x1298 = x22 - x300
x1299 = x1245*x530
x1300 = x1102*x305 + x1102*x316 + x1209*x545 + x1248*x545 + x546*x911 + x67*(x327*x539 + x334*x537 + x338*x534 + x81*(-x1106 - x1298 + x18*x186*x424*x530*x6 - x536) + x82*(-x1098 - x1297) + x83*(-x1253*x530 + x1296 + x533)) + x796*(-x1057 - 2*x1100 - x1105*x1268 - 2*x1299 + x543)
x1301 = x1109*x1220
x1302 = x1116*x1220
x1303 = x1245*x549
x1304 = x1115*x305 + x1115*x316 + x1209*x563 + x1248*x563 + x564*x911 + x67*(x327*x551 + x334*x557 + x338*x555 + x81*(-x1117 + x18*x186*x424*x549*x6 - x550 - x966) + x82*(-x1111 - x1301) + x83*(x1227 - x1253*x549 + x1302 + x554)) + x796*(-x1103 - 2*x1113 - x1116*x1268 - x1267 - 2*x1303 + x561)
x1305 = x1120*x1220
x1306 = x1126*x1220
x1307 = x1245*x567
x1308 = x1125*x305 + x1125*x316 + x1209*x581 + x1248*x581 + x582*x911 + x67*(x327*x573 + x334*x570 + x338*x575 + x81*(-x1127 + x18*x186*x424*x567*x6 - x568) + x82*(x1123 - x1305) + x83*(-x1253*x567 + x1298 + x1306 + x572)) + x796*(-2*x1121 - x1126*x1268 - x1292 - 2*x1307 + x579)
x1309 = x92*(x1130*x374 + x382*x585 + x95*(2*x385 + 2*x387 + 2*x389 + 2*x391 + 2*x393 + 2*x396))
x1310 = x106*x683
x1311 = x1310*x27
x1312 = x100*(x242 - x243 + x244 + x31*x377 + x336 + x337 + x380 + x381)
x1313 = x1310*x135
x1314 = x1310*x162
x1315 = x1314*x14
x1316 = -3*p2
x1317 = x1316 + x682
x1318 = -x1317
x1319 = x1318*x780
x1320 = x1319*x15
x1321 = x1319*x14
x1322 = x1321*x918
x1323 = x1321*x2
x1324 = -x1323*x788 + x191
x1325 = x1310*x29
x1326 = x1320*x44
x1327 = x1321*x44
x1328 = x1310*x14
x1329 = x1327*x353
x1330 = x1310*x17
x1331 = x1330*x264
x1332 = x1310*x881 - x1327*x267
x1333 = x122*x1319
x1334 = x1319*x2
x1335 = x122*x1334
x1336 = x775*(-x228 + x24 - x321 - x373)
x1337 = x1319*x18
x1338 = x1334*x828
x1339 = x1337*x44
x1340 = x1321*x804
x1341 = x1319*x44
x1342 = x1341*x354
x1343 = x1334*x44
x1344 = x115*x1343 - x1328*x349
x1345 = x122*x1310
x1346 = x122*x1341
x1347 = -4*p2
x1348 = x1347 + x748
x1349 = -x1348
x1350 = x1319*x19
x1351 = x1334*x786
x1352 = x1350*x44
x1353 = x124*x349
x1354 = -2*p2
x1355 = x1354 + x740
x1356 = -x1355
x1357 = x1228 + x976
x1358 = x13*x1319
x1359 = x19*x959
x1360 = x1359*x4
x1361 = -x1360*x6 + x428
x1362 = x1246 + x1361 + x958
x1363 = x13*x1321
x1364 = x372 + x436
x1365 = x1068 + x1216 - x1359*x7 + x1364
x1366 = x1235 + x839 + x990
x1367 = x13*x1341
x1368 = x1002*x925 - x1199*x13 + x987
x1369 = x346 + x931
x1370 = -x1002*x893 + x1231 + x13*x938 + x1369
x1371 = x384*x444 + x386*x442 + x388*x440 + x390*x458 + x392*x453 + x395*x448
x1372 = 2*x995
x1373 = x1002*x1353 + x1002*x935 - x1372 - x419 + x857 + x940 + 2*x996
x1374 = x124*x1319
x1375 = 2*x978
x1376 = x1359*x6
x1377 = x1376*x9
x1378 = x1092 + 2*x1221 - x1375 + 2*x1377 + x948 + 2*x979
x1379 = x382*x99
x1380 = x374*x774
x1381 = x1211*x631
x1382 = x368*x438
x1383 = x382*x631
x1384 = x1018*x375 + x1018*x383 + x1379*x420 + x1380*x438 + x1381*x420 + x1382*x1383 + x421*x953 + x439*x953
x1385 = x1319*x16
x1386 = x1029 + x1221 + x227
x1387 = x1025 + x1275 + x1364 - x186*x19*x424*x8
x1388 = x1052 + x1055 + x1361
x1389 = x1341*x16
x1390 = x1042 + x1263 + x861
x1391 = x1039 + x1199*x16 + x1369 - x16*x186*x2*x28*x31*x34
x1392 = x1050*x893 + x1260 - x14*x16*x186*x32*x34
x1393 = x384*x484 + x386*x482 + x388*x486 + x390*x492 + x392*x489 + x395*x495
x1394 = 2*x1044
x1395 = x1050*x1353 + x1050*x935 + x1183 + x1197 - x1394 + 2*x16*x807 - x472
x1396 = 2*x1031
x1397 = x1203 + 2*x1252 + x1293 + 2*x1360*x9 - x1396 + 2*x975
x1398 = x1383*x368
x1399 = x1060*x375 + x1060*x383 + x1379*x473 + x1380*x480 + x1381*x473 + x1398*x480 + x474*x953 + x481*x953
x1400 = x1*x1319
x1401 = -x191
x1402 = -x1055 + x1277 + x1401
x1403 = x1064 + x1255 + x1290 - x186*x19*x4*x424*x9 + x218
x1404 = x1087 + x1089 - x1377 + x194 + x982
x1405 = x1*x1341
x1406 = x851 + x924
x1407 = x1076 + x1284 + x1406
x1408 = -x1*x1199 - x1073 + x1085*x925 + x504 + x886 + x991
x1409 = x1085*x893
x1410 = x384*x514 + x386*x518 + x388*x516 + x390*x520 + x392*x524 + x395*x522
x1411 = 2*x1*x807 - 2*x1079 + x1085*x1353 + x1085*x935 + x291 + 2*x379 - x507 + x812 + x814
x1412 = 2*x372
x1413 = 2*x1359
x1414 = -2*x1*x12*x186*x9 + x10*x1413 - 2*x10*x21*x424 + x1013 + x1014 + x1270 - x1412 + x435
x1415 = x1094*x375 + x1094*x383 + x1379*x508 + x1380*x512 + x1381*x508 + x1398*x512 + x509*x953 + x513*x953
x1416 = x22 - x372
x1417 = x1413*x9
x1418 = x1102*x375 + x1102*x383 + x1380*x545 + x1398*x545 + x546*x953 + x67*(x390*x537 + x392*x539 + x395*x534 + x81*(-x1097 + x1297) + x82*(x1107 - x1376*x530 + x1416 + x536) + x83*(-x1299 + x186*x19*x4*x424*x530 - x531)) + x796*(-x1091 - 2*x1101 - 2*x1296 - x1417*x530 + x544)
x1419 = x1115*x375 + x1115*x383 + x1380*x563 + x1398*x563 + x564*x953 + x67*(x390*x557 + x392*x551 + x395*x555 + x81*(-x1112 + x1301) + x82*(x1118 - x1376*x549 + x550) + x83*(-x1303 - x1416 + x186*x19*x4*x424*x549 - x553)) + x796*(-2*x1110 - x1292 - 2*x1302 - x1417*x549 + x562)
x1420 = x1125*x375 + x1125*x383 + x1380*x581 + x1398*x581 + x582*x953 + x67*(x390*x570 + x392*x573 + x395*x575 + x81*(x1122 + x1305) + x82*(x1128 - x1376*x567 + x568 + x981) + x83*(-x1227 - x1307 + x186*x19*x4*x424*x567 - x569)) + x796*(-x1103 - 2*x1124 - 2*x1306 - x1412 - x1417*x567 + x580)
x1421 = x92*(x1130*x438 + x420*x585 + x95*(2*x441 + 2*x443 + 2*x445 + 2*x449 + 2*x454 + 2*x459))
x1422 = x100*(x26*x413 + x28*x415 + x31*x418 + x331)
x1423 = x780*x782
x1424 = x1423*x15
x1425 = x14*x1423
x1426 = x1425*x2
x1427 = x1426*x788
x1428 = x1425*x918
x1429 = x1424*x44
x1430 = x1425*x44
x1431 = x1430*x353
x1432 = x1430*x267
x1433 = x122*x1423
x1434 = x1423*x2
x1435 = x122*x1434
x1436 = x775*(-x13*x432 - x13*x434 - x13*x436 - x14*x426 - x42 + x422 - x429 - x431)
x1437 = x1423*x18
x1438 = -x1434*x828
x1439 = x1425*x804
x1440 = x1437*x44
x1441 = x1423*x44
x1442 = x1441*x354
x1443 = x1434*x44
x1444 = -x115*x1443
x1445 = x122*x1441
x1446 = x1423*x19
x1447 = x1446*x44
x1448 = x1434*x786
x1449 = x13*x1441
x1450 = x13*x1423
x1451 = 4*n0 + x777
x1452 = 3*n0 + x781
x1453 = x11**(mpmath.mpf(-5)/mpmath.mpf(2))
x1454 = x1453*x21
x1455 = x1452*x1454
x1456 = x14*x1455
x1457 = x13*x960
x1458 = x124*x1455
x1459 = x1455*x2
x1460 = x6*x9
x1461 = x1459*x1460
x1462 = x124*x1423
x1463 = x17*x959
x1464 = x13*x1463
x1465 = x1464*x974
x1466 = x1219*x13
x1467 = x1460*x1466
x1468 = x14*x427
x1469 = 2*x1468
x1470 = x17*x477
x1471 = 2*x1470
x1472 = x2*x532
x1473 = 2*x1472
x1474 = -x1469 - x1471 - x1473
x1475 = -4*x12*x21 + x1474 - x433 - x435 - x437
x1476 = x13*x1425
x1477 = x1459*x974
x1478 = x1455*x17
x1479 = -x1460*x1478
x1480 = x124*x959
x1481 = x13*x1480
x1482 = 2*x13
x1483 = x1482*x974
x1484 = x17*x532
x1485 = x2*x477
x1486 = x1484 - x1485
x1487 = x1177 + x1486
x1488 = x13*x217
x1489 = 2*x1488
x1490 = 2*n0 + x825
x1491 = -x1490
x1492 = x1456*x974
x1493 = x17*x427
x1494 = x14*x477
x1495 = x1493 - x1494
x1496 = x1244 + x1483*x960
x1497 = x187*x2
x1498 = x13*x1497
x1499 = 2*x1498
x1500 = x2*x425
x1501 = x1456*x1460
x1502 = x2*x427
x1503 = x14*x532
x1504 = x1502 - x1503
x1505 = x1482*x960
x1506 = x1375 + x1460*x1505
x1507 = x1423*x16
x1508 = x16*x217
x1509 = x13*x8
x1510 = x16*x974
x1511 = x1510*x960
x1512 = x13*x980
x1513 = x1468 + x1512
x1514 = x1216 + x1513
x1515 = -x1024 + x1465 - x1470 - x1508 - x1509*x960 - x1511 + x1514 + x16*x959*x961 - x476
x1516 = x4*x9
x1517 = x1498 + x1502
x1518 = x1460*x16
x1519 = x1463*x1518 + x430
x1520 = -x13*x186*x2*x424*x8 + x1464*x1516 + x1517 + x1519 - x16*x186*x2*x4*x424*x6 + x978
x1521 = x1425*x16
x1522 = x1457*x1516 + x478
x1523 = x1456*x1516 + x1522
x1524 = x1219*x7
x1525 = x1518*x960
x1526 = x1497*x16
x1527 = x1485 + x1526
x1528 = x1031 - x1466*x974 - x1524*x16 + x1525 + x1527
x1529 = x1441*x16
x1530 = x1000 - x1038
x1531 = x440*x495 + x442*x489 + x444*x492 + x448*x486 + x453*x482 + x458*x484
x1532 = x1049 + x1239
x1533 = x124*x8
x1534 = x1459*x1516
x1535 = 2*x1493
x1536 = -x1535
x1537 = 2*x1494
x1538 = -x1537
x1539 = x16*x980
x1540 = 2*x1539
x1541 = x1010*x960
x1542 = 2*x1219
x1543 = x1219*x1482
x1544 = x1053 + x1480*x1510
x1545 = x1480*x1509 - x1489 + x1496 + x1516*x1543 + x1518*x1542 + x1536 + x1538 - x1540 + x1541*x16 + x1544 + 4*x428
x1546 = x420*x99
x1547 = x438*x774
x1548 = x368*x480
x1549 = x420*x631
x1550 = x1382*x631
x1551 = x1018*x474 + x1018*x481 + x1060*x421 + x1060*x439 + x1546*x473 + x1547*x480 + x1548*x1549 + x1550*x473
x1552 = x1*x974
x1553 = x1552*x960
x1554 = x1*x217
x1555 = x1484 + x1554
x1556 = -x1*x17*x186*x424*x7 + x1064 - x13*x17*x186*x424*x6*x9 + x1553 + x1555
x1557 = x10*x17
x1558 = x1*x1423
x1559 = x1557*x959
x1560 = x1*x1460
x1561 = x1488 + x1493
x1562 = x1219*x1552 + x428
x1563 = -x13*x1559 - x1463*x1560 + x1466*x1516 + x1561 + x1562 + x958
x1564 = x1*x1497
x1565 = -x1069 - x1472 - x1564 - x511
x1566 = -x1*x14*x186*x424*x6*x9 + x1*x1524 - x10*x13*x14*x186*x424 + x1467 + x1514 + x1565
x1567 = x1*x1441
x1568 = -x1079
x1569 = x1000 + x1568
x1570 = x440*x522 + x442*x524 + x444*x520 + x448*x516 + x453*x518 + x458*x514
x1571 = x1084 + x1372
x1572 = 2*x10
x1573 = 2*x1502
x1574 = -x1573
x1575 = 2*x1503
x1576 = -x1575
x1577 = x1*x980
x1578 = 2*x1577
x1579 = x10*x1219
x1580 = x1088 + x1542*x1560
x1581 = x1*x1541 + x1480*x1552 + x1481*x1516 + x1482*x1579 - x1499 + x1506 + x1574 + x1576 - x1578 + x1580 + 4*x430
x1582 = x1549*x368
x1583 = x1018*x509 + x1018*x513 + x1094*x421 + x1094*x439 + x1546*x508 + x1547*x512 + x1550*x508 + x1582*x512
x1584 = x17*x425
x1585 = x1584*x530
x1586 = x1585 + x531
x1587 = x1105*x1457 - x13*x17*x186*x424*x530*x6 - x1488 + x1494 + x1586
x1588 = x530*x9
x1589 = -x1484 + x1485
x1590 = x1105*x1466 - x1464*x1588 + x1589
x1591 = x1500*x530
x1592 = x1591 + x533
x1593 = -x1099*x1466 + x1457*x1588 - x1498 + x1503 + x1592
x1594 = x440*x534 + x442*x539 + x444*x537
x1595 = 2*x1456
x1596 = 2*x1459
x1597 = x1014 + x1473
x1598 = x1099*x1505 + x1105*x1481 - 2*x14*x21*x424*x530 + x1471 + 2*x1512 + x1543*x1588 + x1597 + x542
x1599 = x1018*x546 + x1102*x421 + x1102*x439 + x1547*x545 + x1582*x545
x1600 = x1513 + x22
x1601 = x1584*x549
x1602 = x1601 + x553
x1603 = x1116*x1457 - x13*x17*x186*x424*x549*x6 + x1600 + x1602
x1604 = x549*x9
x1605 = x1116*x1466 - x1464*x1604 + x1517
x1606 = x1500*x549
x1607 = x1606 + x554
x1608 = -x1109*x1466 + x1457*x1604 + x1607
x1609 = x440*x555 + x442*x551 + x444*x557
x1610 = x14*x425
x1611 = x1610*x549
x1612 = 2*x1611
x1613 = x1109*x1505 + x1116*x1481 + x1489 + x1535 + x1543*x1604 - x1612 + x560
x1614 = x1018*x564 + x1115*x421 + x1115*x439 + x1547*x563 + x1582*x563
x1615 = x1584*x567
x1616 = x1615 + x569
x1617 = x1126*x1457 - x13*x17*x186*x424*x567*x6 + x1616
x1618 = x567*x9
x1619 = -x13*x186*x2*x4*x424*x567 + x1464*x1618 + x1561
x1620 = x1500*x567
x1621 = x1620 + x572
x1622 = -x1120*x1466 + x1457*x1618 + x1600 + x1621
x1623 = x440*x575 + x442*x573 + x444*x570
x1624 = x1610*x567
x1625 = 2*x1624
x1626 = x1120*x1505 + x1126*x1481 + x1499 + x1543*x1618 + x1573 - x1625 + x578
x1627 = x1018*x582 + x1125*x421 + x1125*x439 + x1547*x581 + x1582*x581
x1628 = x92*(x1130*x480 + x473*x585 + x95*(2*x483 + 2*x485 + 2*x487 + 2*x490 + 2*x493 + 2*x496))
x1629 = x100*(x236 + x26*x467 + x28*x469 + x31*x471)
x1630 = x1139*x780
x1631 = x15*x1630
x1632 = x14*x1630
x1633 = x1632*x918
x1634 = x1632*x2
x1635 = -x1634*x788
x1636 = x1631*x44
x1637 = x1632*x44
x1638 = x1637*x267
x1639 = -x1637*x353
x1640 = x122*x1630
x1641 = x1630*x2
x1642 = x122*x1641
x1643 = x775*(-x16*x432 - x16*x434 - x16*x436 - x17*x476 - x23 - x455 + x456 - x479)
x1644 = x1630*x18
x1645 = x1641*x828
x1646 = x1644*x44
x1647 = x1632*x804
x1648 = x1630*x44
x1649 = x1648*x354
x1650 = x1641*x44
x1651 = x115*x1650
x1652 = x122*x1648
x1653 = x1630*x19
x1654 = x1653*x44
x1655 = x1641*x786
x1656 = 3*n1 + x1138
x1657 = x1454*x1656
x1658 = x13*x1630
x1659 = x14*x1657
x1660 = x1659*x974
x1661 = x1657*x17
x1662 = x1460*x1661
x1663 = x1657*x2
x1664 = x1663*x974
x1665 = x1460*x1659
x1666 = x13*x1632
x1667 = x13*x1648
x1668 = x124*x1657
x1669 = x1460*x1663
x1670 = x124*x1630
x1671 = x16*x1648
x1672 = x16*x1630
x1673 = 4*n1 + x1160
x1674 = x16*x8
x1675 = 4*x1463
x1676 = x1516*x1663
x1677 = x1516*x16
x1678 = x1219*x1677
x1679 = -x1516*x1659
x1680 = -x1502 + x1503
x1681 = x1172 + x1680 + x194
x1682 = 2*n1 + x1174
x1683 = 2*x960
x1684 = -x1493 + x1494
x1685 = x1684 + x830
x1686 = 2*x1526
x1687 = -x1682
x1688 = x1396 + x1480*x1677
x1689 = x1*x1630
x1690 = x1*x8
x1691 = x1503 + x1577
x1692 = x1087 + x1463*x1552 + x1519 - x1677*x960 - x1690*x960 + x1691
x1693 = x1*x1516
x1694 = x1470 + x1508
x1695 = x1025 + x1219*x1690 - x1463*x1693 - x1559*x16 + x1565 + x1678 + x1694
x1696 = x1494 + x1539
x1697 = -x1*x14*x186*x4*x424*x9 - x10*x14*x16*x186*x424 + x1052 + x1219*x1518 + x1562 + x1696
x1698 = x1*x1648
x1699 = x1038 + x1568
x1700 = x482*x524 + x484*x520 + x486*x522 + x489*x518 + x492*x514 + x495*x516
x1701 = x1287 + x1394
x1702 = 2*x1485
x1703 = -x1702
x1704 = 2*x1484
x1705 = -x1704
x1706 = 2*x1554
x1707 = x1219*x16
x1708 = x1289 + x1542*x1693
x1709 = x1480*x1690 + 2*x1525 + 2*x1553 + x1572*x1707 - x1686 + x1688 + x1703 + x1705 - x1706 + x1708 + 4*x478
x1710 = x480*x774
x1711 = x368*x512
x1712 = x473*x631
x1713 = x508*x631
x1714 = x1060*x509 + x1060*x513 + x1094*x474 + x1094*x481 + x1548*x1713 + x1710*x512 + x1711*x1712 + x473*x508*x99
x1715 = x1463*x16
x1716 = x16*x960
x1717 = x1694 + x22
x1718 = x1610*x530
x1719 = x1718 + x536
x1720 = x1099*x1715 - x1105*x1716 + x1717 + x1719
x1721 = x1588*x1715 + x1592 - x16*x186*x2*x4*x424*x530
x1722 = x1099*x1707 - x14*x16*x186*x424*x530*x9 + x1527
x1723 = x482*x539 + x484*x537 + x486*x534
x1724 = 2*x1659
x1725 = 2*x1663
x1726 = 2*x1585
x1727 = x16*x1683
x1728 = x1480*x16
x1729 = x1542*x16
x1730 = x1099*x1727 + x1105*x1728 + x1537 + x1540 + x1588*x1729 - x1726 + x543
x1731 = x1712*x368
x1732 = x1060*x546 + x1102*x474 + x1102*x481 + x1710*x545 + x1731*x545
x1733 = x1611 + x550
x1734 = x1109*x1715 - x1116*x1716 + x1493 - x1539 + x1733
x1735 = x1484 - x1526 - x16*x186*x2*x4*x424*x549 + x1604*x1715 + x1607
x1736 = x1109*x1707 - x14*x16*x186*x424*x549*x9 + x1504
x1737 = x482*x551 + x484*x557 + x486*x555
x1738 = x1109*x1727 + x1116*x1728 + x1469 + 2*x1508 + x1597 + x1604*x1729 - 2*x17*x21*x424*x549 + x561
x1739 = x1060*x564 + x1115*x474 + x1115*x481 + x1710*x563 + x1731*x563
x1740 = x1624 + x568
x1741 = x1120*x1715 - x1126*x1716 + x1740
x1742 = -x16*x186*x2*x4*x424*x567 + x1618*x1715 + x1621 + x1717
x1743 = -x1120*x1707 + x1618*x1716 + x1696
x1744 = x482*x573 + x484*x570 + x486*x575
x1745 = 2*x1615
x1746 = x1120*x1727 + x1126*x1728 + x1618*x1729 + x1686 + x1702 - x1745 + x579
x1747 = x1060*x582 + x1125*x474 + x1125*x481 + x1710*x581 + x1731*x581
x1748 = x92*(x1130*x512 + x508*x585 + x95*(2*x515 + 2*x517 + 2*x519 + 2*x521 + 2*x523 + 2*x525))
x1749 = x100*(x243 + x26*x502 + x28*x504 + x31*x506)
x1750 = x1317*x780
x1751 = x15*x1750
x1752 = x14*x1750
x1753 = x1752*x2
x1754 = x1753*x788
x1755 = -x1752*x918
x1756 = x1751*x44
x1757 = x1752*x44
x1758 = x1757*x353
x1759 = -x1757*x267
x1760 = x122*x1750
x1761 = x1750*x2
x1762 = x122*x1761
x1763 = x775*(-x1*x432 - x1*x434 - x1*x436 - x2*x511 - x24 - x450 + x451 - x494)
x1764 = x1750*x18
x1765 = x1761*x828
x1766 = x1764*x44
x1767 = x1752*x804
x1768 = x1750*x44
x1769 = x1768*x354
x1770 = x1761*x44
x1771 = x115*x1770
x1772 = x122*x1768
x1773 = x1750*x19
x1774 = x1773*x44
x1775 = x1761*x786
x1776 = 3*n2 + x1316
x1777 = x1454*x1776
x1778 = x14*x1777
x1779 = x1778*x974
x1780 = x1777*x2
x1781 = x1780*x974
x1782 = x17*x1777
x1783 = x1460*x1782
x1784 = x13*x1750
x1785 = x13*x1768
x1786 = x124*x1777
x1787 = x1460*x1780
x1788 = x124*x1750
x1789 = x16*x1750
x1790 = -x1516*x1778
x1791 = x16*x1768
x1792 = x1516*x1780
x1793 = x1*x1768
x1794 = x1*x1750
x1795 = 4*n2 + x1347
x1796 = 2*n2 + x1354
x1797 = -x1796
x1798 = x1*x1463
x1799 = x1*x1105
x1800 = x1099*x1798 + x1555 - x1799*x960
x1801 = x1219*x1799 + x1586 - x1588*x1798
x1802 = x1*x1219
x1803 = x1472 + x1564 + x22
x1804 = -x1*x14*x186*x424*x530*x9 + x1099*x1802 + x1719 + x1803
x1805 = x514*x537 + x516*x534 + x518*x539
x1806 = 2*x1778
x1807 = 2*x1780
x1808 = 2*x1591
x1809 = x1*x1683
x1810 = x1*x1542
x1811 = x1099*x1809 + x1480*x1799 + x1575 + x1578 + x1588*x1810 - x1808 + x544
x1812 = x512*x774
x1813 = x1713*x368
x1814 = x1094*x546 + x1102*x509 + x1102*x513 + x1812*x545 + x1813*x545
x1815 = x1*x960
x1816 = -x1*x17*x186*x424*x549*x6 + x1116*x1815 + x1691
x1817 = x1116*x1802 + x1602 - x1604*x1798 + x1803
x1818 = -x1*x14*x186*x424*x549*x9 + x1109*x1802 + x1733
x1819 = x514*x557 + x516*x555 + x518*x551
x1820 = 2*x1606
x1821 = x1*x1480
x1822 = x1109*x1809 + x1116*x1821 + x1604*x1810 + x1704 + x1706 - x1820 + x562
x1823 = x1094*x564 + x1115*x509 + x1115*x513 + x1812*x563 + x1813*x563
x1824 = x1120*x1798 - x1126*x1815 + x1495
x1825 = x1126*x1802 + x1485 - x1554 + x1616 - x1618*x1798
x1826 = -x1*x14*x186*x424*x567*x9 + x1120*x1802 + x1502 - x1577 + x1740
x1827 = x514*x570 + x516*x575 + x518*x573
x1828 = x1014 + x1120*x1809 + x1126*x1821 + x1469 + x1471 + 2*x1564 + x1618*x1810 - 2*x2*x21*x424*x567 + x580
x1829 = x1094*x582 + x1125*x509 + x1125*x513 + x1812*x581 + x1813*x581
x1830 = x92*(x1130*x545 + x95*(2*x535 + 2*x538 + 2*x540))
x1831 = x775*(-x14*x536 - x17*x531 - x2*x533 - x422)
x1832 = -x1452
x1833 = x1454*x1832
x1834 = x14*x1833
x1835 = x1834*x974
x1836 = x1833*x2
x1837 = x1836*x974
x1838 = x17*x1833
x1839 = x1460*x1838
x1840 = x1460*x1834
x1841 = x124*x1833
x1842 = x1460*x1836
x1843 = -x1516*x1834
x1844 = x1516*x1836
x1845 = 2*x1834
x1846 = 2*x1836
x1847 = -x1718
x1848 = x1601 + x1847
x1849 = x1612 + x1726
x1850 = x545*x774
x1851 = x1102*x564 + x1115*x546 + x1850*x563
x1852 = x1620 + x1847
x1853 = x1625 + x1808
x1854 = x1102*x582 + x1125*x546 + x1850*x581
x1855 = x92*(x1130*x563 + x95*(2*x552 + 2*x556 + 2*x558))
x1856 = x775*(-x14*x550 - x17*x553 - x2*x554 - x456)
x1857 = -x1656
x1858 = x1454*x1857
x1859 = x14*x1858
x1860 = x1859*x974
x1861 = x1858*x2
x1862 = x1861*x974
x1863 = x17*x1858
x1864 = -x1460*x1863
x1865 = x1460*x1859
x1866 = x124*x1858
x1867 = x1460*x1861
x1868 = x1516*x1859
x1869 = x1516*x1861
x1870 = 2*x1859
x1871 = 2*x1861
x1872 = -x1601 + x1620
x1873 = x1745 + x1820
x1874 = x1115*x582 + x1125*x564 + x563*x581*x774
x1875 = x92*(x1130*x581 + x95*(2*x571 + 2*x574 + 2*x576))
x1876 = x775*(-x14*x568 - x17*x569 - x2*x572 - x451)
x1877 = -x1776
x1878 = x1454*x1877
x1879 = x14*x1878
x1880 = x1879*x974
x1881 = x17*x1878
x1882 = x1460*x1881
x1883 = x1878*x2
x1884 = -x1883*x974
x1885 = x1460*x1879
x1886 = x124*x1878
x1887 = x1460*x1883
x1888 = x1516*x1879
x1889 = x1516*x1883
x1890 = 2*x1879
x1891 = 2*x1883
return np.array([-x86*x96 - x97*(x101*x76 + x134*(x108*x123 + x109*x124 + x110*x125 + x132 + 4*x70) + x57*x80 + x67*(x25*(x109*x2 - x110*x17 + x112) + x43*(x108*x115 - x109*x14 + x114 + x117) + x51*(x104*x105*x14*x21*x27*x31 - x108*x118 - x121 - x75)) + x76**2*x99), -x150*x96 - x97*(x101*x148 + x134*(x123*x151 + x124*x152 + x125*x153 + x157) + x161 + x67*(x25*(x152*x2 - x153*x17 + x74) + x43*(x115*x151 - x14*x152 + x156) + x51*(x104*x105*x135*x14*x21*x31 - x118*x151 - x146))), -x177*x96 - x97*(x101*x175 + x134*(x123*x179 + x124*x180 + x125*x181 + x183) + x184 + x67*(x25*(x104*x105*x162*x2*x21*x28 - x113 - x17*x181) + x43*(x115*x179 - x14*x180 + x178) + x51*(x104*x105*x14*x162*x21*x31 - x118*x179 - x182))), -x249*x96 - x97*(x101*x214 + x134*(x122*x252 + x284*x69 + 2*x285 + 2*x286 + x293) + x298 + x67*(x25*(-x266 + x269 + x274) + x282 + x43*(x17*x252 + x251*x38 - x255 + x263) + x51*(x119 - x2*x252 + x259*x31 + x277 + x278 + x279*x53 - x280 + x281 + x74))), -x341*x96 - x97*(x101*x315 + x134*(x124*x342 + 2*x255 + 2*x362 + x367) + x370 + x67*(x25*(x2*x342 - x351 + x357) + x361 + x43*(-x14*x342 + x285 + x348) + x51*(-x269 - x360 - x37))), -x398*x96 - x97*(x101*x382 + x134*(x125*x400 + 2*x281 + 2*x351 + x411) + x412 + x67*(x25*(-x17*x400 + x362 + x404) + x407 + x43*(x266 + x360 + x399) + x51*(x104*x105*x14*x21*x32 - x286 - x406))), -x461*x96 - x97*(x101*x420 + x466), -x498*x96 - x97*(x101*x473 + x501), -x527*x96 - x97*(x101*x508 + x529), -x547*x96 - x548, -x565*x96 - x566, -x583*x96 - x584, -x586*x86 - x97*(x134*(x123*x591 + x157 + x464*x591 + x591*x592) + x161 + x593*x76 + x67*(x25*(-x267*x591 + x353*x591 + x74) + x43*(x115*x591 - x14*x28*x591 + x156) + x51*(x105*x14*x21*x27*x31*x589 - x118*x591 - x146))), -x150*x586 - x97*(x134*(x123*x594 + x132 + 4*x154 + x464*x594 + x592*x594) + x142*x149 + x148**2*x99 + x148*x593 + x67*(x25*(x112 - x267*x594 + x353*x594 + x596) + x43*(x105*x135*x17*x21*x26*x589 - x14*x28*x594 - x145 - x595) + x51*(x105*x135*x14*x21*x31*x589 - x118*x594 - x121))), -x177*x586 - x97*(x134*(x123*x597 + x464*x597 + x592*x597 + x601) + x175*x593 + x605 + x67*(x25*(-x267*x597 + x353*x597 + x599) + x43*(x105*x162*x17*x21*x26*x589 - x171 - x28*x598) + x51*(-x118*x597 + x144 + x31*x598))), -x249*x586 - x97*(x134*(x122*x618 + x124*x606 + x125*x615 + x630) + x214*x593 + x633 + x67*(x25*(x54 + x614 + x617) + x43*(x105*x17*x21*x29*x589 - x14*x606 - x611) + x51*(x105*x14*x21*x26*x31*x589 - x2*x618 - x623) + x624)), -x341*x586 - x97*(x134*(x122*x606 + x124*x638 + x125*x641 + x38*x653 + x656) + x315*x593 + x658 + x67*(x25*(-x17*x641 + x2*x638 + x53*x635 + x644) + x43*(-x14*x638 + x17*x606 - x210 + x26*x637 + x47 + x636*x69 - x639 + x640) + x51*(-x614 + x645 - x650) + x651)), -x398*x586 - x97*(x134*(x122*x615 + x124*x641 + x125*x660 + x668) + x382*x593 + x669 + x67*(x25*(-x17*x660 + x2*x641 + x662) + x43*(-x617 - x645 - x659) + x51*(x105*x14*x21*x32*x589 - x2*x615 - x663) + x664)), -x461*x586 - x97*(x420*x593 + x673), -x498*x586 - x97*(x473*x593 + x675), -x527*x586 - x97*(x508*x593 + x677), -x547*x586 - x678, -x565*x586 - x679, -x583*x586 - x680, -x681*x86 - x97*(x134*(x123*x686 + x183 + x464*x686 + x592*x686) + x184 + x67*(x25*(x105*x2*x21*x27*x28*x684 - x113 - x267*x686) + x43*(x115*x686 - x14*x28*x686 + x178) + x51*(x105*x14*x21*x27*x31*x684 - x118*x686 - x182)) + x687*x76), -x150*x681 - x97*(x134*(x123*x688 + x464*x688 + x592*x688 + x601) + x148*x687 + x605 + x67*(x25*(-x267*x688 + x353*x688 + x599) + x43*(x105*x135*x17*x21*x26*x684 - x171 - x28*x689) + x51*(-x118*x688 + x144 + x31*x689))), -x177*x681 - x97*(x134*(x123*x690 + x132 + 4*x173 + x464*x690 + x592*x690) + x169*x176 + x175**2*x99 + x175*x687 + x67*(x25*(x105*x162*x2*x21*x28*x684 - x267*x690 - x273 - x600) + x43*(x115*x690 + x117 - x28*x691) + x51*(-x118*x690 + x172 + x31*x691 + x649))), -x249*x681 - x97*(x134*(x122*x692 + x124*x693 + x125*x700 + x716) + x214*x687 + x67*(x25*(-x47 - x702 - x704) + x43*(-x14*x693 + x17*x692 + x699) + x51*(x14*x700 - x2*x692 + x709) + x710) + x721), -x341*x681 - x97*(x134*(x122*x693 + x124*x722 + x125*x725 + x735) + x315*x687 + x67*(x25*(x105*x2*x21*x30*x684 - x17*x725 - x726) + x43*(-x14*x722 + x17*x693 + x724) + x51*(x704 + x727 + x731) + x732) + x737), -x398*x681 - x97*(x134*(x122*x700 + x124*x725 + x125*x744 + x53*x749 + x751) + x382*x687 + x67*(x25*(-x17*x744 + x178 + x2*x725 + x28*x743 + x37 + x38*x742 - x697 + x739 - x745) + x43*(x702 - x727 + x738) + x51*(x14*x744 - x2*x700 + x69*x741 + x746) + x747) + x752), -x461*x681 - x97*(x420*x687 + x756), -x498*x681 - x97*(x473*x687 + x758), -x527*x681 - x97*(x508*x687 + x762), -x547*x681 - x763, -x565*x681 - x764, -x583*x681 - x765, -x767*x86 - x97*(x134*(x123*x769 + 4*x126 + x293 + x464*x769 + x592*x769) + x298 + x67*(x25*(-x267*x769 + x268 + x274 + x353*x769) + x282 + x43*(x115*x769 - x14*x28*x769 + x263 + x625) + x51*(x103*x105*x14*x21*x27*x31 - x118*x769 + x15*x186*x27*x31*x34 - x212 - x280 - x356 - x648 - x711)) + x76*x770), -x150*x767 - x97*(x134*(x123*x771 + x464*x771 + x592*x771 + x630) + x148*x770 + x633 + x67*(x25*(-x267*x771 + x353*x771 + x54 + x613 + x616) + x43*(x103*x105*x135*x17*x21*x26 - x14*x28*x771 - x611) + x51*(x103*x105*x135*x14*x21*x31 - x118*x771 - x623) + x624)), -x177*x767 - x97*(x134*(x123*x772 + x464*x772 + x592*x772 + x716) + x175*x770 + x67*(x25*(x103*x105*x162*x2*x21*x28 - x267*x772 - x47 - x701 - x703) + x43*(x115*x772 - x28*x773 + x699) + x51*(-x118*x772 + x31*x773 + x709) + x710) + x721), -x249*x767 - x97*(x134*(x122*x31*x806 + x122*x800 + x124*x802 + x125*x803 + 4*x15*x799 - 4*x204 - 4*x206 - x284*x69 + x289*x801 + x778*x798 + x790*x805 + x808 + x818) + x196**2*x774 + x196*x776 + x200*x248 + x214**2*x99 + x214*x295*x797 + x214*x770 + x215*x248 + x67*(2*x220*x234 + 2*x224*x240 + 2*x229*x246 + x25*(x833 - x834 - x835 + x836 + x837 + x842) + x43*(x115*x832 - x14*x802 + x17*x800 - x211 + x279*x38 - x28*x801 + x827*x845 + x844 - x847 + x849 + x854) + x51*(-x118*x832 + x14*x803 - x2*x800 + x213 + x251*x53 + x31*x801 + x826*x855 - x857 - x859 + x860 + x865) + x81*(x190*x827 - x785*x788 + x820*x828 + x830) + x82*(x193*x826 + x785*x791 - x786*x821 + x831) + x83*(x819 + x822 + x824)) + x796*(-x122*x793 + x226*x779 - x785*x787 - x789*x790 + x795)), -x341*x767 - x97*(x134*(x115*x892 + x122*x802 + x124*x871 + x125*x885 + x365*x870 + x655*x875 + x903) + x303*x776 + x315*x770 + x67*(x25*(x2*x871 - x31*x870 - x768*x877 + x876 + x884) + x43*(-x14*x871 + x17*x802 + x26*x870 - x804*x823 + x874) + x51*(-x837 - x888 - x890) + x81*(x786*x866 - x788*x823 + x868) + x82*(-x824 - x869) + x83*(x194 + x789*x792 - x791*x866) + x891) + x796*(-x124*x793 - x790*x828 - x866*x904 - x909) + x917), -x398*x767 - x97*(x134*(x118*x892 + x122*x803 + x124*x885 + x410*x922 + x768*x935 + 2*x876 + x944) + x374*x776 + x382*x770 + x67*(x25*(x2*x885 - x267*x806 + x28*x922 - x768*x923 + x930) + x43*(-x833 + x834 + x888 + x921) + x51*(x103*x105*x14*x21*x32 + x14*x2*x31*x44*x780*x783 - x2*x803 - x26*x922 - x933) + x81*(x320 - x822 + x869) + x82*(x12*x14*x2*x780*x783*x9 - x786*x919 - x920) + x83*(x12*x19*x4*x780*x783 - x191 - x792*x918) + x934) + x796*(-x122*x946 - x792*x947 - x919*x945 - x952) + x955), -x461*x767 - x97*(x1019 + x134*(x1006 + x123*x988 + x592*x988 + 4*x811 + x956*x999) + x420*x770 + x438*x776 + x67*(x25*(-x267*x988 + x353*x988 + x993) + x43*(x13*x17*x186*x29*x34 + x13*x17*x26*x44*x780*x783 - x28*x985 - x843 - x852 - x898 - x984 - x987) + x51*(-x118*x988 + x31*x985 + x941 + x997) + x81*(-x820*x957 + x828*x956 + x908 + x968) + x82*(x12*x13*x14*x780*x783*x9 - x13*x946 - x951 - x983) + x83*(x792*x957 - x918*x956 + x977) + x998) + x796*(-x1007*x956 - x1008*x956 - x1017 - 4*x432 - x947*x956)), -x498*x767 - x97*(x1061 + x134*(x1020*x999 + x1035*x123 + x1035*x592 + x1051) + x473*x770 + x480*x776 + x67*(x1047 + x25*(-x1035*x267 - x1043 + x16*x2*x28*x44*x780*x783) + x43*(x1035*x115 + x1040 - x16*x28*x832) + x51*(-x1035*x118 - x1046 + x14*x16*x31*x44*x780*x783) + x81*(x1020*x828 - x1021*x820 + x1027) + x82*(-x1034 + x12*x14*x16*x780*x783*x9 - x16*x946) + x83*(-x1020*x918 - x1030 + x12*x16*x2*x4*x780*x783)) + x796*(-x1007*x1020 - x1008*x1020 - x1020*x947 - x1059)), -x527*x767 - x97*(x1095 + x134*(x1062*x999 + x1072*x123 + x1072*x592 + x1086) + x508*x770 + x512*x776 + x67*(x1082 + x25*(-x1072*x267 + x1072*x353 + x1078) + x43*(-x1*x28*x832 + x1072*x115 + x1075) + x51*(x1*x14*x31*x44*x780*x783 - x1072*x118 - x1081) + x81*(x1062*x828 - x1063*x820 + x1065) + x82*(x1*x12*x14*x780*x783*x9 - x1*x946 - x1071) + x83*(-x1062*x918 + x1063*x792 + x1067)) + x796*(-x1007*x1062 - x1008*x1062 - x1062*x947 - x1093)), -x547*x767 - x97*(x1108 + x545*x776), -x565*x767 - x97*(x1119 + x563*x776), -x583*x767 - x97*(x1129 + x581*x776), -x1131*x86 - x97*(x1134*x76 + x134*(x1133*x123 + x1133*x464 + x1133*x592 + x367) + x370 + x67*(x25*(-x1133*x267 + x1133*x353 + x357) + x361 + x43*(x1133*x115 - x1133*x14*x28 + x348) + x51*(x105*x14*x21*x27*x31*x588 - x1133*x118 - x268 - x359 - x37))), -x1131*x150 - x97*(x1134*x148 + x134*(x1135*x123 + x1135*x464 + x1135*x592 + 4*x128 + x656) + x658 + x67*(x25*(-x1135*x267 + x1135*x353 + x644 + x733) + x43*(x105*x135*x17*x21*x26*x588 - x1135*x14*x28 - x116 + x135*x18*x186*x26*x34 + x144 - x210 - x363 - x626 - x639) + x51*(x105*x135*x14*x21*x31*x588 - x1135*x118 - x613 - x650) + x651)), -x1131*x177 - x97*(x1134*x175 + x134*(x1136*x123 + x1136*x464 + x1136*x592 + x735) + x67*(x25*(x105*x162*x2*x21*x28*x588 - x1136*x267 - x726) + x43*(x1136*x115 - x1137*x28 + x724) + x51*(-x1136*x118 + x1137*x31 + x703 + x731) + x732) + x737), -x1131*x249 - x97*(x1134*x214 + x1159*x196 + x134*(x1003*x1132 + x1004*x1132 + x1147*x122 + x1148*x289 + x1155*x805 + x1156*x1158 + x903) + x67*(x25*(x1151 - x1153 + x1154 + x884) + x43*(x1147*x17 - x1148*x28 + x1149*x115 - x1150*x253 + x874) + x51*(x105*x14*x21*x26*x31*x588 + x1140*x15*x31*x44*x780 - x1147*x2 - x1149*x118 - x887 - x890) + x81*(-x1142*x788 + x1143*x828 + x868) + x82*(x1140*x12*x15*x780*x9 - x1144*x786 - x218) + x83*(x1145 + x1146) + x891) + x796*(-x1142*x787 - x1155*x789 - x1158*x791 - x909) + x917), -x1131*x341 - x97*(x1134*x315 + x1159*x303 + x134*(x1132*x1165 + x115*x1170 + x1161*x845 + x1166*x365 + x1167*x253 + x1168*x125 + x1169*x655 + x1171 + 4*x30*x343 - 4*x307 - 4*x309 - x38*x653 + x818) + x303**2*x774 + x305*x340 + x315**2*x99 + x315*x797*x915 + x316*x340 + x67*(x25*(-x1132*x877 + x1132*x878 - x1166*x31 + x1176*x855 + x1183 - x1184 + x1186 + x1187 + x53*x636 - x667 + x842) + 2*x318*x327 + 2*x322*x338 + 2*x324*x334 + x43*(-x1132*x1181 + x1152*x253 + x1166*x26 + x1175*x798 - x1182*x17 + x211 + x635*x69 - x844 + x854 + x896 - x900) + x51*(-x1154 - x1178 - x1180 - x121 + 2*x14*x17*x186*x28*x31*x34 - x937) + x81*(-x1143*x789 + x1163*x786 + x1175*x226 + x830) + x82*(-x1146 - x1172 - x1173) + x83*(x1157*x789 - x1163*x791 + x1176*x193 + x1177)) + x796*(-x1155*x828 - x1157*x1164 + x1162*x190 - x1163*x904 + x795)), -x1131*x398 - x97*(x1134*x382 + x1159*x374 + x1213 + x134*(x1132*x935 + x1167*x264 + x1168*x124 + x1170*x118 + 2*x1187 + x1192*x410 + x1201) + x67*(x1195 + x25*(-x1132*x923 + x1132*x925 - x1179*x267 + x1192*x28 + x1193) + x43*(-x1151 + x1153 + x1180 + x1191) + x51*(x105*x14*x21*x32*x588 - x1132*x893 + x1140*x14*x2*x31*x44*x780 - x1192*x26 - x1194) + x81*(-x1145 + x1173 - x227) + x82*(x1144*x791 - x1188*x786 + x222) + x83*(-x1157*x918 + x1188*x788 + x1189)) + x796*(-x1157*x947 - x1188*x945 - x1202*x122 - x1207)), -x1131*x461 - x97*(x1134*x420 + x1159*x438 + x1249 + x134*(x1218*x1238 + x1218*x999 + x123*x1233 + x1242) + x67*(x1237 + x25*(-x1233*x267 + x1233*x353 + x1234) + x43*(x1140*x13*x17*x26*x44*x780 - x1182*x13 - x1232) + x51*(x1156*x1226 - x118*x1233 + x1236) + x81*(x1140*x12*x13*x17*x6*x780 - x1143*x957 - x1217) + x82*(-x1202*x13 + x1226*x791 + x1229) + x83*(x1157*x957 - x1218*x918 + x1225)) + x796*(-x1007*x1218 - x1008*x1218 - x1243*x957 - x1247)), -x1131*x498 - x97*(x1134*x473 + x1159*x480 + x1272 + x134*(x123*x1258 + x1238*x1251 + x1251*x999 + x1266 + 4*x813) + x67*(x1265 + x25*(x1044 + x1140*x16*x2*x28*x44*x780 - x1198 - x1258*x267 - x1262 + x16*x186*x2*x30*x34 - x471 - x840 - x886) + x43*(x115*x1258 - x1182*x16 + x1261 + x897) + x51*(x1140*x14*x16*x31*x44*x780 - x118*x1258 - x1264) + x81*(-x1021*x1143 + x1140*x12*x16*x17*x6*x780 - x1250 - x906) + x82*(x1140*x12*x14*x16*x780*x9 - x1202*x16 - x1257) + x83*(x1021*x1157 + x1206 - x1251*x918 + x1256)) + x796*(-x1007*x1251 - x1008*x1251 - x1021*x1243 - x1271 - 4*x434)), -x1131*x527 - x97*(x1134*x508 + x1159*x512 + x1295 + x134*(x123*x1279 + x1238*x1274 + x1274*x999 + x1288) + x67*(x1286 + x25*(-x1279*x267 + x1279*x353 + x1283) + x43*(-x1*x1182 + x115*x1279 + x1282) + x51*(x1*x1140*x14*x31*x44*x780 - x118*x1279 - x1285) + x81*(x1*x1140*x12*x17*x6*x780 - x1063*x1143 - x1273) + x82*(x1*x1140*x12*x14*x780*x9 - x1*x1202 - x1278) + x83*(x1063*x1157 - x1274*x918 + x1276)) + x796*(-x1007*x1274 - x1008*x1274 - x1063*x1243 - x1294)), -x1131*x547 - x97*(x1159*x545 + x1300), -x1131*x565 - x97*(x1159*x563 + x1304), -x1131*x583 - x97*(x1159*x581 + x1308), -x1309*x86 - x97*(x1312*x76 + x134*(x123*x1311 + x1311*x464 + x1311*x592 + x411) + x412 + x67*(x25*(-x1311*x267 + x1311*x353 + x404) + x407 + x43*(x115*x1311 - x1311*x14*x28 + x359 + x399) + x51*(x105*x14*x21*x27*x31*x683 - x118*x1311 - x406))), -x1309*x150 - x97*(x1312*x148 + x134*(x123*x1313 + x1313*x464 + x1313*x592 + x668) + x669 + x67*(x25*(-x1313*x267 + x1313*x353 + x662) + x43*(x105*x135*x17*x21*x26*x683 - x1313*x14*x28 - x616 - x659) + x51*(x105*x135*x14*x21*x31*x683 - x118*x1313 - x663) + x664)), -x1309*x177 - x97*(x1312*x175 + x134*(x123*x1314 + 4*x130 + x1314*x464 + x1314*x592 + x751) + x67*(x25*(x105*x162*x2*x21*x28*x683 - x111 - x1314*x267 + x162*x186*x19*x28*x34 - x313 - x665 - x698 - x745) + x43*(x115*x1314 - x1315*x28 + x701 + x738) + x51*(-x118*x1314 + x1315*x31 + x408 + x746) + x747) + x752), -x1309*x249 - x97*(x1312*x214 + x1336*x196 + x134*(x1003*x1310 + x1004*x1310 + x1156*x1335 + x122*x1325 + x1326*x289 + x1333*x805 + x944) + x67*(x25*(x1329 - x1331 + x1332 + x930) + x43*(x115*x1327 + x1325*x17 - x1326*x28 - x1328*x253 + x887 + x921) + x51*(x105*x14*x21*x26*x31*x683 - x118*x1327 + x1318*x15*x31*x44*x780 - x1325*x2 - x933) + x81*(-x1320*x788 + x1321*x828 + x320) + x82*(x12*x1318*x15*x780*x9 - x1323*x786 - x920) + x83*(-x1322 - x1324) + x934) + x796*(-x1320*x787 - x1333*x789 - x1335*x791 - x952) + x955), -x1309*x341 - x97*(x1213 + x1312*x315 + x1336*x303 + x134*(x115*x1346 + x1165*x1310 + x1201 + x1241*x1310 + x1339*x365 + x1341*x655 + x1345*x253) + x67*(x1195 + x25*(x1193 - x1310*x877 + x1310*x878 - x1339*x31 + x1342) + x43*(-x1181*x1310 + x1191 + x1330*x253 + x1339*x26 - x1340*x17) + x51*(-x1194 - x1332 - x1344) + x81*(x12*x1318*x18*x6*x780 - x1321*x789 - x227) + x82*(x1322 - x1338 + x222) + x83*(x1189 + x1334*x789 - x1337*x791)) + x796*(-x1164*x1334 - x1207 - x1333*x828 - x1337*x904)), -x1309*x398 - x97*(x1211*x382*x797 + x1312*x382 + x1336*x374 + x134*(x1171 + x118*x1346 + x1310*x1353 + x1310*x935 + x132 + 2*x1342 + x1345*x264 + x1348*x855 + x1352*x410 + 4*x32*x402 - 4*x377 - 4*x379 - x53*x749 + x808 + x809 + x817) + x374**2*x774 + x375*x397 + x382**2*x99 + x383*x397 + x67*(x25*(-x1183 + x1196 - x1200 - x1310*x923 + x1310*x925 - x1343*x267 + x1352*x28 + x1355*x845 + x38*x741 + x667 + x842) + 2*x384*x390 + 2*x386*x392 + 2*x388*x395 + x43*(-x1329 + x1331 + x1344 + x854 + x894 - x902) + x51*(x1156*x1323 + x1310*x14*x32 - x1310*x893 - x1352*x26 + x1356*x798 - x213 + x69*x742 + x857 + x865 + x939 - x943) + x81*(x1324 + x1338 + x829) + x82*(x1323*x791 - x1350*x786 + x1356*x226 + x831) + x83*(x1177 - x1334*x918 + x1350*x788 + x1355*x190)) + x796*(-x122*x1351 - x1334*x947 + x1349*x193 - x1350*x945 + x795)), -x1309*x461 - x97*(x1312*x420 + x1336*x438 + x134*(x123*x1367 + x1238*x1358 + x1358*x999 + x1373) + x1384 + x67*(x1371 + x25*(-x1367*x267 + x1367*x353 + x1368) + x43*(x13*x1318*x17*x26*x44*x780 - x13*x1340 - x1366) + x51*(x1156*x1363 - x118*x1367 + x1370) + x81*(x12*x13*x1318*x17*x6*x780 - x1321*x957 - x1357) + x82*(-x13*x1351 + x1363*x791 + x1365) + x83*(x12*x13*x1318*x2*x4*x780 - x1358*x918 - x1362)) + x796*(-x1007*x1358 - x1008*x1358 - x1374*x957 - x1378)), -x1309*x498 - x97*(x1312*x473 + x1336*x480 + x134*(x123*x1389 + x1238*x1385 + x1385*x999 + x1395) + x1399 + x67*(x1393 + x25*(x1318*x16*x2*x28*x44*x780 - x1389*x267 - x1391) + x43*(x115*x1389 - x1340*x16 + x1390) + x51*(-x118*x1389 + x1318*x14*x16*x31*x44*x780 - x1392) + x81*(-x1021*x1321 + x1385*x828 + x1386) + x82*(x1321*x16*x791 - x1351*x16 + x1388) + x83*(x12*x1318*x16*x2*x4*x780 - x1385*x918 - x1387)) + x796*(-x1007*x1385 - x1008*x1385 - x1021*x1374 - x1397)), -x1309*x527 - x97*(x1312*x508 + x1336*x512 + x134*(x123*x1405 + x1238*x1400 + x1400*x999 + x1411 + 4*x815) + x1415 + x67*(x1410 + x25*(x1197 - x1405*x267 + x1405*x353 + x1408) + x43*(-x1*x1340 + x115*x1405 + x1407) + x51*(x1*x1318*x14*x31*x44*x780 + x1*x14*x186*x32*x34 - x118*x1405 - x1281 - x1409 - x856 - x863 - x940) + x81*(-x1063*x1321 + x1400*x828 + x1402) + x82*(x1*x1321*x791 - x1*x1351 + x1404 + x949) + x83*(x1*x12*x1318*x2*x4*x780 - x1204 - x1400*x918 - x1403)) + x796*(-x1007*x1400 - x1008*x1400 - x1063*x1374 - x1414 - 4*x436)), -x1309*x547 - x97*(x1336*x545 + x1418), -x1309*x565 - x97*(x1336*x563 + x1419), -x1309*x583 - x97*(x1336*x581 + x1420), -x1421*x86 - x97*(x1422*x76 + x466), -x1421*x150 - x97*(x1422*x148 + x673), -x1421*x177 - x97*(x1422*x175 + x756), -x1421*x249 - x97*(x1019 + x134*(x1006 + x1156*x1435 + x1429*x289 + x1433*x805 + x779*x798) + x1422*x214 + x1436*x196 + x67*(x25*(x1431 - x1432 + x993) + x43*(x115*x1430 + x13*x848 - x1429*x28 - x415 + x826*x845 + x851 + x928 - x984 + x986) + x51*(-x118*x1430 + x1429*x31 + x827*x855 + x997) + x81*(-x1424*x788 + x1425*x828 + x190*x826 + x968) + x82*(x12*x15*x780*x782*x9 + x12*x186*x827*x9 - x1426*x786 - x983) + x83*(x1427 - x1428 + x977) + x998) + x796*(-x1017 + x12*x186*x6*x778 - x1424*x787 - x1433*x789 - x1435*x791)), -x1421*x341 - x97*(x1249 + x134*(x115*x1445 + x1242 + x1440*x365 + x1441*x655) + x1422*x315 + x1436*x303 + x67*(x1237 + x25*(x1234 - x1440*x31 + x1442) + x43*(-x1232 - x1439*x17 + x18*x26*x44*x780*x782) + x51*(x1236 + x1432 + x1444) + x81*(x12*x18*x6*x780*x782 - x1217 - x1425*x789) + x82*(x1229 + x1428 + x1438) + x83*(x1225 + x1434*x789 - x1437*x791)) + x796*(-x1164*x1434 - x1247 - x1433*x828 - x1437*x904)), -x1421*x398 - x97*(x134*(x118*x1445 + x1373 + 2*x1442 + x1447*x410) + x1384 + x1422*x382 + x1436*x374 + x67*(x1371 + x25*(x1368 - x1443*x267 + x1447*x28) + x43*(-x1366 - x1431 - x1444) + x51*(x1156*x1426 + x1370 - x1447*x26) + x81*(-x1357 - x1427 - x1438) + x82*(x1365 + x1426*x791 - x1446*x786) + x83*(x12*x19*x4*x780*x782 - x1362 - x1434*x918)) + x796*(-x122*x1448 - x1378 - x1434*x947 - x1446*x945)), -x1421*x461 - x97*(x134*(4*x1000 + x123*x1449 + x1238*x1450 + x1450*x999 + x817) + x1382*x420*x797 + x1422*x420 + x1436*x438 + x420**2*x99 + x421*x460 + x438**2*x774 + x439*x460 + x67*(x25*(-x1449*x267 + x1449*x353 + x841) + x43*(-x1239 - x13*x1439 + x13*x17*x26*x44*x780*x782 - x1406) + 2*x440*x448 + 2*x442*x453 + 2*x444*x458 + x51*(x1156*x1476 - x118*x1449 + x1372 + x864) + x81*(-x1056 + x12*x13*x17*x6*x780*x782 + 2*x13*x17*x186*x424*x7 - x1401 - x1425*x957 + x1452*x1453*x17*x21*x7 - x1489 + x1491*x17*x21*x424 - x1492 - x1495 - x1496 - x222) + x82*(-x1010*x1466 + x1090 - x13*x1448 - x1459*x7 + x1476*x791 + x1490*x1500 + x1499 + x1501 + x1504 + x1506 + x831) + x83*(x1219*x1483 + x1434*x957 - x1450*x918 - x1460*x1481 + x1477 + x1479 + x1487)) + x796*(-x1007*x1450 - x1008*x1450 - x1010*x1456 + 4*x12*x13*x14*x186 - 4*x1215 + x14*x1451*x21*x424 - 4*x1457*x7 - x1458*x974 - 2*x1461 - x1462*x957 - 4*x1465 - 4*x1467 - x1475 - 4*x426)), -x1421*x498 - x97*(x134*(x123*x1529 + x1238*x1507 + x1507*x999 + x1532) + x1422*x473 + x1436*x480 + x1551 + x67*(x1531 + x25*(-x1529*x267 + x16*x2*x28*x44*x780*x782 - x995) + x43*(x115*x1529 - x1439*x16 + x1530) + x51*(x1044 + x1156*x1521 - x118*x1529) + x81*(-x1021*x1425 - x1456*x8 + x1478*x974 + x1507*x828 + x1515) + x82*(-x1448*x16 - x1477 + x1521*x791 + x1523 + x1528) + x83*(x12*x16*x2*x4*x780*x782 + x1452*x1453*x2*x21*x8 - x1478*x1516 - x1507*x918 - x1520)) + x796*(-x1007*x1507 - x1008*x1507 - x1021*x1462 - x1455*x1533 - 2*x1492 - 2*x1534 - x1545)), -x1421*x527 - x97*(x134*(x123*x1567 + x1238*x1558 + x1558*x999 + x1571) + x1422*x508 + x1436*x512 + x1583 + x67*(x1570 + x25*(-x1567*x267 + x1567*x353 + x986) + x43*(-x1*x1439 + x1*x17*x26*x44*x780*x782 - x1073) + x51*(x1*x14*x31*x44*x780*x782 - x118*x1567 - x1569) + x81*(x1*x12*x17*x6*x780*x782 - x1063*x1425 - x1479 - x1523 - x1556) + x82*(x1*x12*x14*x780*x782*x9 - x1*x1448 + x10*x14*x1452*x1453*x21 - x1461 - x1566) + x83*(x1063*x1434 - x1455*x1557 + x1534 - x1558*x918 + x1563)) + x796*(-x1007*x1558 - x1008*x1558 - x1063*x1462 - x1458*x1516 - x1459*x1572 - 2*x1501 - x1581)), -x1421*x547 - x97*(x1436*x545 + x1599 + x67*(x1594 + x81*(-x1105*x1456 + x1452*x1453*x17*x21*x530*x6 - x1536 - x1587) + x82*(-x1099*x1459 + x1456*x1588 + x1574 + x1593) + x83*(x1105*x1459 - x1478*x1588 + x1590)) + x796*(-x1099*x1595 - x1105*x1458 - 4*x1468 - x1588*x1596 - x1598)), -x1421*x565 - x97*(x1436*x563 + x1614 + x67*(x1609 + x81*(-x1116*x1456 + x1452*x1453*x17*x21*x549*x6 - x1603) + x82*(-x1109*x1459 + x1456*x1604 + x1608) + x83*(x1116*x1459 - x1478*x1604 + x1605)) + x796*(-x1109*x1595 - x1116*x1458 - x1596*x1604 - x1613)), -x1421*x583 - x97*(x1436*x581 + x1627 + x67*(x1623 + x81*(-x1126*x1456 + x1452*x1453*x17*x21*x567*x6 - x1617) + x82*(-x1120*x1459 + x1456*x1618 + x1622) + x83*(x1452*x1453*x2*x21*x4*x567 - x1478*x1618 - x1619)) + x796*(-x1120*x1595 - x1126*x1458 - x1596*x1618 - x1626)), -x1628*x86 - x97*(x1629*x76 + x501), -x150*x1628 - x97*(x148*x1629 + x675), -x1628*x177 - x97*(x1629*x175 + x758), -x1628*x249 - x97*(x1061 + x134*(x1051 + x1156*x1642 + x1636*x289 + x1640*x805) + x1629*x214 + x1643*x196 + x67*(x1047 + x25*(-x1043 - x1638 - x1639) + x43*(x1040 + x115*x1637 - x1636*x28) + x51*(-x1046 + x1139*x15*x31*x44*x780 - x118*x1637) + x81*(x1027 - x1631*x788 + x1632*x828) + x82*(-x1034 + x1139*x12*x15*x780*x9 - x1634*x786) + x83*(-x1030 - x1633 - x1635)) + x796*(-x1059 - x1631*x787 - x1640*x789 - x1642*x791)), -x1628*x341 - x97*(x1272 + x134*(x115*x1652 + x1162*x845 + x1266 + x1646*x365 + x1648*x655) + x1629*x315 + x1643*x303 + x67*(x1265 + x25*(x1045 + x1175*x855 + x1185*x16 - x1262 - x1646*x31 + x1649 + x839 - x886) + x43*(x1176*x798 + x1261 + x1646*x26 - x1647*x17) + x51*(-x1264 + x1638 - x1651) + x81*(x1139*x12*x18*x6*x780 + x1176*x12*x186*x6 - x1250 - x1632*x789) + x82*(-x1257 + x1633 - x1645) + x83*(x1175*x193 + x1256 + x1641*x789 - x1644*x791)) + x796*(x1161*x12*x186*x4 - x1164*x1641 - x1271 - x1640*x828 - x1644*x904)), -x1628*x398 - x97*(x134*(x118*x1652 + x1395 + 2*x1649 + x1654*x410) + x1399 + x1629*x382 + x1643*x374 + x67*(x1393 + x25*(x1139*x19*x28*x44*x780 - x1391 - x1650*x267) + x43*(x1390 + x1639 + x1651) + x51*(x1139*x14*x2*x31*x44*x780 - x1392 - x1654*x26) + x81*(x1386 + x1635 + x1645) + x82*(x1388 + x1634*x791 - x1653*x786) + x83*(x1139*x12*x19*x4*x780 - x1387 - x1641*x918)) + x796*(-x122*x1655 - x1397 - x1641*x947 - x1653*x945)), -x1628*x461 - x97*(x134*(x123*x1667 + x1238*x1658 + x1532 + x1658*x999) + x1551 + x1629*x420 + x1643*x438 + x67*(x1531 + x25*(x1139*x13*x2*x28*x44*x780 - x1667*x267 - x995) + x43*(x115*x1667 - x13*x1647 + x1530) + x51*(x1044 + x1156*x1666 - x118*x1667) + x81*(x1515 - x1632*x957 + x1657*x961 + x1658*x828 - x1660) + x82*(-x13*x1655 + x1522 + x1528 - x1663*x7 + x1665 + x1666*x791) + x83*(x1139*x12*x13*x2*x4*x780 - x1520 - x1658*x918 - x1662 + x1664)) + x796*(-x1007*x1658 - x1008*x1658 - x1010*x1659 - x1545 - x1668*x974 - 2*x1669 - x1670*x957)), -x1628*x498 - x97*(x134*(4*x1038 + x123*x1671 + x1238*x1672 + x1672*x999 + x817) + x1548*x473*x797 + x1629*x473 + x1643*x480 + x473**2*x99 + x474*x497 + x480**2*x774 + x481*x497 + x67*(x25*(x1139*x16*x2*x28*x44*x780 - x1394 - x1671*x267 - x992) + x43*(x1049 + x115*x1671 - x16*x1647 + x853) + 2*x482*x489 + 2*x484*x492 + 2*x486*x495 + x51*(x1139*x14*x16*x31*x44*x780 - x1178 - x118*x1671) + x81*(-x1021*x1632 + x1056 + x1540 + x1544 + x1610*x1682 - x1659*x8 + x1661*x974 + x1672*x828 - x1674*x1683 + x1685) + x82*(x1139*x12*x14*x16*x780*x9 - 2*x1219*x1510 + 2*x14*x16*x186*x4*x424*x9 - x16*x1655 - x1664 - x1679 - x1681) + x83*(x1139*x12*x16*x2*x4*x780 - x1291 + x1453*x1656*x2*x21*x8 - x1516*x1661 - x1589 + 2*x16*x186*x2*x424*x8 - x1672*x918 - x1686 + x1687*x2*x21*x424 - x1688 - x320 - x969)) + x796*(-x1007*x1672 - x1008*x1672 - x1021*x1670 - 4*x1024 + 4*x12*x16*x17*x186 - x1475 - 4*x1511 - x1533*x1657 - 2*x1660 + x1673*x17*x21*x424 - x1674*x1675 - 2*x1676 - 4*x1678 - 4*x476)), -x1628*x527 - x97*(x134*(x123*x1698 + x1238*x1689 + x1689*x999 + x1701) + x1629*x508 + x1643*x512 + x1714 + x67*(x1700 + x25*(-x1698*x267 + x1698*x353 + x1699) + x43*(-x1*x1647 + x1083 + x115*x1698) + x51*(x1*x1139*x14*x31*x44*x780 - x1048 - x118*x1698) + x81*(-x1063*x1632 + x1662 + x1679 + x1689*x828 + x1692) + x82*(x1*x1139*x12*x14*x780*x9 - x1*x1655 + x10*x14*x1453*x1656*x21 - x1669 - x1697) + x83*(x1063*x1641 - x1557*x1657 + x1676 - x1689*x918 + x1695)) + x796*(-x1007*x1689 - x1008*x1689 - x1063*x1670 - x1516*x1668 - x1572*x1663 - 2*x1665 - x1709)), -x1628*x547 - x97*(x1643*x545 + x1732 + x67*(x1723 + x81*(x1099*x1661 - x1105*x1659 + x1720) + x82*(-x1099*x1663 + x14*x1453*x1656*x21*x530*x9 - x1722) + x83*(x1453*x1656*x2*x21*x4*x530 - x1588*x1661 - x1721)) + x796*(-x1099*x1724 - x1105*x1668 - x1588*x1725 - x1730)), -x1628*x565 - x97*(x1643*x563 + x1739 + x67*(x1737 + x81*(x1109*x1661 - x1116*x1659 + x1538 + x1734) + x82*(-x1109*x1663 + x14*x1453*x1656*x21*x549*x9 - x1736) + x83*(x1453*x1656*x2*x21*x4*x549 - x1604*x1661 - x1703 - x1735)) + x796*(-x1109*x1724 - x1116*x1668 - 4*x1470 - x1604*x1725 - x1738)), -x1628*x583 - x97*(x1643*x581 + x1747 + x67*(x1744 + x81*(x1120*x1661 - x1126*x1659 + x1741) + x82*(-x1120*x1663 + x1618*x1659 + x1743) + x83*(x1453*x1656*x2*x21*x4*x567 - x1618*x1661 - x1742)) + x796*(-x1120*x1724 - x1126*x1668 - x1618*x1725 - x1746)), -x1748*x86 - x97*(x1749*x76 + x529), -x150*x1748 - x97*(x148*x1749 + x677), -x1748*x177 - x97*(x1749*x175 + x762), -x1748*x249 - x97*(x1095 + x134*(x1086 + x1156*x1762 + x1756*x289 + x1760*x805) + x1749*x214 + x1763*x196 + x67*(x1082 + x25*(x1078 + x1758 + x1759) + x43*(x1075 + x115*x1757 - x1756*x28) + x51*(-x1081 - x118*x1757 + x1317*x15*x31*x44*x780) + x81*(x1065 - x1751*x788 + x1752*x828) + x82*(-x1071 + x12*x1317*x15*x780*x9 - x1753*x786) + x83*(x1067 + x1754 + x1755)) + x796*(-x1093 - x1751*x787 - x1760*x789 - x1762*x791)), -x1748*x341 - x97*(x1295 + x134*(x115*x1772 + x1288 + x1766*x365 + x1768*x655) + x1749*x315 + x1763*x303 + x67*(x1286 + x25*(x1283 - x1766*x31 + x1769) + x43*(x1282 - x17*x1767 + x1766*x26) + x51*(-x1285 - x1759 - x1771) + x81*(x12*x1317*x18*x6*x780 - x1273 - x1752*x789) + x82*(-x1278 - x1755 - x1765) + x83*(x1276 + x1761*x789 - x1764*x791)) + x796*(-x1164*x1761 - x1294 - x1760*x828 - x1764*x904)), -x1748*x398 - x97*(x134*(x118*x1772 + x1349*x855 + x1411 + 2*x1769 + x1774*x410) + x1415 + x1749*x382 + x1763*x374 + x67*(x1410 + x25*(x1356*x845 + x1408 - x1770*x267 + x1774*x28) + x43*(x1407 - x1758 + x1771) + x51*(x1*x938 + x1083 + x1156*x1753 + x1190 + x1355*x798 - x1409 - x1774*x26 + x760 + x862) + x81*(x1402 - x1754 + x1765) + x82*(x1355*x226 + x1404 + x1753*x791 - x1773*x786) + x83*(x12*x1317*x19*x4*x780 + x12*x1356*x186*x4 - x1403 - x1761*x918)) + x796*(x12*x1348*x186*x9 - x122*x1775 - x1414 - x1761*x947 - x1773*x945)), -x1748*x461 - x97*(x134*(x123*x1785 + x1238*x1784 + x1571 + x1784*x999) + x1583 + x1749*x420 + x1763*x438 + x67*(x1570 + x25*(-x1785*x267 + x1785*x353 + x986) + x43*(-x1073 + x13*x1317*x17*x26*x44*x780 - x13*x1767) + x51*(-x118*x1785 + x13*x1317*x14*x31*x44*x780 - x1569) + x81*(x12*x13*x1317*x17*x6*x780 + x1453*x17*x1776*x21*x7 - x1522 - x1556 - x1752*x957 - x1779) + x82*(x12*x13*x1317*x14*x780*x9 - x13*x1775 + x14*x1453*x1776*x21*x6*x9 - x1566 - x1780*x7) + x83*(x1563 + x1761*x957 + x1781 - x1783 - x1784*x918)) + x796*(-x1007*x1784 - x1008*x1784 - x1010*x1778 - x1581 - x1786*x974 - 2*x1787 - x1788*x957)), -x1748*x498 - x97*(x134*(x123*x1791 + x1238*x1789 + x1701 + x1789*x999) + x1714 + x1749*x473 + x1763*x480 + x67*(x1700 + x25*(x1699 - x1791*x267 + x1791*x353) + x43*(x1083 + x115*x1791 - x16*x1767) + x51*(-x1048 - x118*x1791 + x1317*x14*x16*x31*x44*x780) + x81*(-x1021*x1752 + x1692 - x1778*x8 + x1782*x974 + x1789*x828) + x82*(x12*x1317*x14*x16*x780*x9 - x16*x1775 - x1697 - x1781 - x1790) + x83*(x1021*x1761 - x1516*x1782 + x1695 + x1780*x8 - x1789*x918)) + x796*(-x1007*x1789 - x1008*x1789 - x1021*x1788 - x1533*x1777 - x1709 - 2*x1779 - 2*x1792)), -x1748*x527 - x97*(x134*(4*x1079 + x123*x1793 + x1238*x1794 + x1794*x999 + x817) + x1711*x508*x797 + x1749*x508 + x1763*x512 + x508**2*x99 + x509*x526 + x512**2*x774 + x513*x526 + x67*(x25*(x1287 - x1793*x267 + x1793*x353 + x841) + x43*(-x1*x1767 + x115*x1793 + x853) + x51*(x1*x1317*x14*x31*x44*x780 - x1084 - x1178 - x118*x1793) + 2*x514*x520 + 2*x516*x522 + 2*x518*x524 + x81*(-x1063*x1752 + x1480*x1560 - x1683*x1693 + x1685 + x1783 + x1790 + x1794*x828) + x82*(2*x1*x10*x14*x186*x424 + x1*x12*x1317*x14*x780*x9 - x1*x1775 + x10*x14*x1453*x1776*x21 - x1090 + x14*x1797*x21*x424 - x1578 - x1580 - x1681 - x1787) + x83*(-2*x1*x1559 + x1063*x1761 + x1291 + x1487 - x1557*x1777 + x1584*x1796 + x1706 + x1708 + x1792 - x1794*x918)) + x796*(4*x1*x12*x186*x2 - 4*x1*x1579 - x1007*x1794 - x1008*x1794 - x1063*x1788 - 4*x1069 - 2*x1460*x1778 - x1475 - x1516*x1786 - 4*x1560*x960 - x1572*x1780 - x1675*x1693 + x1795*x2*x21*x424 - 4*x511)), -x1748*x547 - x97*(x1763*x545 + x1814 + x67*(x1805 + x81*(x1099*x1782 - x1105*x1778 + x1800) + x82*(-x1099*x1780 + x14*x1453*x1776*x21*x530*x9 - x1804) + x83*(x1105*x1780 - x1588*x1782 + x1801)) + x796*(-x1099*x1806 - x1105*x1786 - x1588*x1807 - x1811)), -x1748*x565 - x97*(x1763*x563 + x1823 + x67*(x1819 + x81*(-x1116*x1778 + x1453*x17*x1776*x21*x549*x6 - x1816) + x82*(-x1109*x1780 + x14*x1453*x1776*x21*x549*x9 - x1818) + x83*(x1116*x1780 - x1604<