d15fidc:=[
((1/7200*(-a + 1)*r1 + 1/1440*(a - 1))*b2^2 + (1/720*(-a + 2)*r1 + 1/144*(a - 2))*b2 + (1/360*(18*a - 31)*r1 + 1/72*(-18*a + 31)))*b4 + (1/7200*(-a + 1)*r1 + 1/1440*(a - 1))*b2^2 + (1/720*(-a + 2)*r1 + 1/144*(a - 2))*b2 + 1/180*(6*a - 11)*r1 + 1/36*(-6*a + 11),
((((((1/172800*(-43*a + 76)*r1 + 1/34560*(-7*a + 12))*t + (1/86400*(17*a - 29)*r1 + 1/34560*(10*a - 17)))*b2^2 + ((1/34560*(-7*a + 13)*r1 + 1/34560*(37*a - 63))*t + (1/34560*(-5*a + 8)*r1 + 1/34560*(-16*a + 29)))*b2 + ((1/8640*(193*a - 334)*r1 + 1/8640*(-97*a + 168))*t + (1/17280*(-176*a + 305)*r1 + 1/17280*(-148*a + 257))))*b3 + (((1/57600*(-49*a + 85)*r1 + 1/11520*(-17*a + 29))*t + (1/86400*(43*a - 75)*r1 + 1/34560*(67*a - 117)))*b2^2 + ((1/17280*(113*a - 198)*r1 + 1/5760*(11*a - 18))*t + (1/34560*(-89*a + 153)*r1 + 1/34560*(-383*a + 669)))*b2 + ((1/8640*(19*a - 33)*r1 + 1/2880*(97*a - 167))*t + (1/5760*(-23*a + 41)*r1 + 1/5760*(-147*a + 257)))))*b4 + ((((1/172800*(3*a + 2)*r1 + 1/34560*(-a + 2))*t + (1/172800*(a - 1)*r1 + 1/34560*a))*b2^2 + ((1/34560*(-3*a + 5)*r1 + 1/34560*(a + 1))*t + (1/34560*(a - 4)*r1 + 1/34560*(2*a - 3)))*b2 + ((1/8640*(69*a - 128)*r1 + 1/8640*(-55*a + 92))*t + (1/8640*(-28*a + 49)*r1 + 1/2160*(-5*a + 6))))*b3 + (((1/172800*(-87*a + 151)*r1 + 1/34560*(-31*a + 55))*t + (1/86400*(28*a - 49)*r1 + 1/34560*(41*a - 69)))*b2^2 + ((1/17280*(51*a - 88)*r1 + 1/17280*(19*a - 32))*t + (1/34560*(-39*a + 67)*r1 + 1/11520*(-63*a + 109)))*b2 + ((1/8640*(51*a - 79)*r1 + 1/8640*(183*a - 347))*t + (1/8640*(-43*a + 77)*r1 + 1/4320*(-97*a + 153))))))*b1 + ((((1/14400*(-a + 3)*r1*t + (1/57600*(4*a - 9)*r1 - 1/11520))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(8*a - 15)*r1*t + (1/576*(2*a - 3)*r1 + 1/576*(-14*a + 25))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(29*a - 51)*r1*t + (1/8640*(43*a - 75)*r1 + 1/1728*(-73*a + 129)))))*b4 + (((1/14400*(-a + 3)*r1*t + (1/57600*(4*a - 9)*r1 - 1/11520))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(5*a - 9)*r1*t + (1/2880*(4*a - 3)*r1 + 1/576*(-8*a + 13))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/1080*(19*a - 30)*r1*t + (1/8640*(97*a - 147)*r1 + 1/1728*(-55*a + 93)))))))*i + (((((1/172800*(-26*a + 47)*r1 + 1/34560*(10*a - 17))*t + (1/57600*(a - 2)*r1 + 1/17280*(-a + 2)))*b2^2 + ((1/34560*(31*a - 55)*r1 + 1/34560*(11*a - 19))*t + (1/11520*(-7*a + 12)*r1 + 1/34560*(-28*a + 47)))*b2 + ((1/8640*(-94*a + 163)*r1 + 1/8640*(-236*a + 409))*t + (1/1152*(15*a - 26)*r1 + 1/17280*(401*a - 694))))*b3 + (((1/57600*(-89*a + 155)*r1 + 1/34560*(-11*a + 17))*t + (1/172800*(91*a - 155)*r1 + 1/8640*(25*a - 44)))*b2^2 + ((1/17280*(-26*a + 45)*r1 + 1/17280*(-94*a + 161))*t + (1/6912*(11*a - 19)*r1 + 1/34560*(139*a - 241)))*b2 + ((1/8640*(335*a - 579)*r1 + 1/8640*(195*a - 331))*t + (1/17280*(-311*a + 535)*r1 + 1/17280*(-1275*a + 2231)))))*b4 + ((((1/172800*(-4*a + 9)*r1 + 1/34560)*t + (1/172800*(2*a - 7)*r1 + 1/34560*(a - 1)))*b2^2 + ((1/34560*(-a - 3)*r1 + 1/34560*(3*a - 7))*t + (1/34560*(-a + 2)*r1 + 1/34560*(-4*a + 1)))*b2 + ((1/8640*(-50*a + 87)*r1 + 1/8640*(-90*a + 151))*t + (1/4320*(26*a - 43)*r1 + 1/8640*(85*a - 148))))*b3 + (((1/172800*(-155*a + 269)*r1 + 1/34560*(-3*a + 5))*t + (1/172800*(55*a - 97)*r1 + 1/17280*(26*a - 45)))*b2^2 + ((1/17280*(-12*a + 19)*r1 + 1/17280*(-40*a + 71))*t + (1/34560*(21*a - 37)*r1 + 1/34560*(65*a - 111)))*b2 + ((1/8640*(199*a - 341)*r1 + 1/8640*(67*a - 121))*t + (1/4320*(-43*a + 77)*r1 + 1/8640*(-341*a + 573))))))*b1 + (((1/43200*(-a + 3)*r1*t + (1/172800*(-11*a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/6912*(-a + 3)*r1 + 1/6912*(-a + 3)))*b2 + (1/432*(-7*a + 12)*r1*t + (1/8640*(155*a - 258)*r1 + 1/1728*(-5*a + 6))))*b3 + ((1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(-23*a + 39)*r1*t + (1/2880*(55*a - 97)*r1 + 1/576*(3*a - 5)))))*b4 + ((1/43200*(-a + 3)*r1*t + (1/172800*(-11*a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/6912*(-a + 3)*r1 + 1/6912*(-a + 3)))*b2 + (1/2160*(-17*a + 33)*r1*t + (1/8640*(83*a - 150)*r1 + 1/1728*(-5*a + 6))))*b3 + (1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(-7*a + 12)*r1*t + 1/2880*(37*a - 73)*r1 + 1/576*(-3*a + 7),
((((((1/172800*(26*a - 47)*r1 + 1/34560*(-4*a + 7))*t + (1/345600*(-11*a + 17)*r1 + 1/23040*(-5*a + 9)))*b2^2 + ((1/34560*(23*a - 41)*r1 + 1/34560*(31*a - 49))*t + (1/34560*(-13*a + 22)*r1 + 1/11520*(-16*a + 29)))*b2 + ((1/4320*(-115*a + 199)*r1 + 1/4320*(-53*a + 92))*t + (1/8640*(94*a - 163)*r1 + 1/2880*(141*a - 244))))*b3 + (((1/57600*(-3*a + 5)*r1 + 1/34560*(-81*a + 139))*t + (1/172800*(102*a - 175)*r1 + 1/34560*(47*a - 82)))*b2^2 + ((1/17280*(-83*a + 144)*r1 + 1/17280*(41*a - 68))*t + (1/34560*(83*a - 145)*r1 + 1/34560*(49*a - 77)))*b2 + ((1/4320*(-17*a + 30)*r1 + 1/4320*(176*a - 301))*t + (1/8640*(-91*a + 157)*r1 + 1/8640*(-91*a + 161)))))*b4 + ((((1/172800*(-6*a + 1)*r1 + 1/11520)*t + (1/345600*(a - 11)*r1 + 1/69120*(5*a - 1)))*b2^2 + ((1/34560*(3*a - 5)*r1 + 1/11520*(a + 1))*t + (1/34560*(-a + 2)*r1 + 1/34560*(-4*a + 11)))*b2 + ((1/8640*(-90*a + 167)*r1 + 1/2880*(-20*a + 27))*t + (1/17280*(83*a - 133)*r1 + 1/17280*(335*a - 631))))*b3 + (((1/172800*(11*a - 13)*r1 + 1/34560*(-37*a + 63))*t + (1/172800*(42*a - 71)*r1 + 1/34560*(19*a - 34)))*b2^2 + ((1/17280*(-51*a + 88)*r1 + 1/17280*(13*a - 24))*t + (1/34560*(51*a - 89)*r1 + 1/34560*(57*a - 101)))*b2 + ((1/8640*(-37*a + 55)*r1 + 1/2880*(57*a - 103))*t + (1/8640*(-36*a + 65)*r1 + 1/8640*(-43*a + 62))))))*b1 + ((((1/14400*(-a + 3)*r1*t + (1/57600*(-5*a + 12)*r1 + 1/11520*(a - 2)))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(8*a - 15)*r1*t + (-1/1440*a*r1 + 1/288*(3*a - 5))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(-29*a + 51)*r1*t + (1/4320*(7*a - 12)*r1 + 1/864*(-22*a + 39)))))*b4 + (((1/14400*(-a + 3)*r1*t + (1/57600*(-5*a + 12)*r1 + 1/11520*(a - 2)))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(5*a - 9)*r1*t + (1/2880*(a - 6)*r1 + 1/576*(3*a - 4))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/1080*(-19*a + 30)*r1*t + (1/8640*(59*a - 87)*r1 + 1/1728*(-17*a + 33)))))))*i + (((((1/172800*(47*a - 84)*r1 + 1/34560*(9*a - 14))*t + (1/345600*(-47*a + 79)*r1 + 1/69120*(-37*a + 67)))*b2^2 + ((1/6912*(-5*a + 9)*r1 + 1/34560*(3*a - 7))*t + (1/34560*(8*a - 13)*r1 + 1/34560*(41*a - 74)))*b2 + ((1/4320*(-7*a + 12)*r1 + 1/2160*(-42*a + 73))*t + (1/8640*(34*a - 59)*r1 + 1/8640*(79*a - 136))))*b3 + (((1/57600*(117*a - 205)*r1 + 1/34560*(7*a - 11))*t + (1/172800*(-203*a + 350)*r1 + 1/34560*(-58*a + 101)))*b2^2 + ((1/1920*(-4*a + 7)*r1 + 1/17280*(-108*a + 187))*t + (1/34560*(95*a - 163)*r1 + 1/34560*(183*a - 317)))*b2 + ((1/1440*(-57*a + 100)*r1 + 1/4320*(-20*a + 29))*t + (1/1728*(39*a - 67)*r1 + 1/8640*(265*a - 469)))))*b4 + ((((1/172800*(3*a - 8)*r1 + 1/34560*(a + 2))*t + (1/115200*(-a + 1)*r1 + 1/69120*(-a + 7)))*b2^2 + ((1/34560*(3*a + 1)*r1 + 1/34560*(-a - 11))*t + (1/11520*r1 + 1/34560*(-7*a - 2)))*b2 + ((1/8640*(9*a - 16)*r1 + 1/8640*(-65*a + 98))*t + (1/5760*(5*a - 9)*r1 + 1/17280*(5*a - 23))))*b3 + (((1/172800*(155*a - 267)*r1 + 1/34560*(3*a - 7))*t + (1/172800*(-95*a + 166)*r1 + 1/34560*(-22*a + 37)))*b2^2 + ((1/17280*(-16*a + 27)*r1 + 1/17280*(-60*a + 107))*t + (1/34560*(47*a - 83)*r1 + 1/34560*(95*a - 157)))*b2 + ((1/2880*(-55*a + 91)*r1 + 1/8640*(-21*a + 53))*t + (1/8640*(99*a - 178)*r1 + 1/8640*(108*a - 179))))))*b1 + (((1/43200*(-a + 3)*r1*t + (1/172800*(10*a - 3)*r1 + 1/34560*(-4*a + 3)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/432*(-7*a + 12)*r1*t + (1/4320*(-95*a + 159)*r1 + 1/864*(20*a - 33))))*b3 + ((1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(23*a - 39)*r1*t + (1/1440*(39*a - 68)*r1 + 1/288*(-10*a + 17)))))*b4 + ((1/43200*(-a + 3)*r1*t + (1/172800*(10*a - 3)*r1 + 1/34560*(-4*a + 3)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/2160*(-17*a + 33)*r1*t + (1/8640*(-100*a + 183)*r1 + 1/1728*(22*a - 39))))*b3 + (1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(7*a - 12)*r1*t + 1/2880*(51*a - 97)*r1 + 1/576*(-17*a + 31),
(((((1/28800*(35*a - 61)*r1 + 1/17280*(-17*a + 31))*t + (1/172800*(-95*a + 163)*r1 + 1/6912*(a - 1)))*b2^2 + ((1/8640*(-41*a + 72)*r1 + 1/8640*(-51*a + 88))*t + (1/17280*(70*a - 121)*r1 + 1/17280*(116*a - 203)))*b2 + ((1/4320*(-43*a + 78)*r1 + 1/4320*(191*a - 344))*t + (1/8640*(6*a - 7)*r1 + 1/8640*(-286*a + 469))))*b4 + (((1/86400*(73*a - 127)*r1 + 1/5760*(-3*a + 5))*t + (1/172800*(-63*a + 107)*r1 + 1/34560*(-11*a + 19)))*b2^2 + ((1/8640*(-17*a + 28)*r1 + 1/8640*(-19*a + 36))*t + (1/17280*(30*a - 53)*r1 + 1/17280*(44*a - 71)))*b2 + ((1/216*(-3*a + 5)*r1 + 1/2160*(47*a - 75))*t + (1/4320*(18*a - 31)*r1 + 1/4320*(-12*a + 35)))))*b1 + (((1/10800*(5*a - 9)*r1*t + (1/86400*(5*a - 9)*r1 + 1/17280*(-5*a + 9)))*b2^2 + (1/1080*(-2*a + 3)*r1*t + (1/8640*(-2*a + 3)*r1 + 1/1728*(2*a - 3)))*b2 + (1/1080*(20*a - 33)*r1*t + (1/4320*(280*a - 489)*r1 + 1/864*(-64*a + 111))))*b4 + ((1/10800*(5*a - 9)*r1*t + (1/86400*(5*a - 9)*r1 + 1/17280*(-5*a + 9)))*b2^2 + (1/1080*(-2*a + 3)*r1*t + (1/8640*(-2*a + 3)*r1 + 1/1728*(2*a - 3)))*b2 + (1/1080*(11*a - 15)*r1*t + (1/2160*(59*a - 105)*r1 + 1/216*(-7*a + 12))))))*i + ((((1/86400*(61*a - 105)*r1 + 1/17280*(27*a - 47))*t + (1/172800*(-151*a + 265)*r1 + 1/34560*(-43*a + 73)))*b2^2 + ((1/8640*(44*a - 75)*r1 + 1/8640*(-40*a + 67))*t + (1/17280*(-25*a + 44)*r1 + 1/17280*(-35*a + 58)))*b2 + ((1/4320*(-166*a + 285)*r1 + 1/4320*(-128*a + 223))*t + (1/8640*(225*a - 394)*r1 + 1/8640*(513*a - 878))))*b4 + (((1/9600*(5*a - 9)*r1 + 1/17280*(19*a - 31))*t + (1/172800*(-95*a + 161)*r1 + 1/6912*(-7*a + 13)))*b2^2 + ((1/8640*(16*a - 27)*r1 + 1/1728*(-4*a + 7))*t + (1/17280*(-9*a + 16)*r1 + 1/3456*(a - 2)))*b2 + ((1/2160*(-47*a + 81)*r1 + 1/540*(-13*a + 22))*t + (1/4320*(68*a - 115)*r1 + 1/4320*(162*a - 293)))))*b1 + ((1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + (1/120*(7*a - 12)*r1*t + (1/1440*(55*a - 92)*r1 + 1/288*(-5*a + 8))))*b4 + (1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + 1/120*(3*a - 5)*r1*t + 1/360*(10*a - 17)*r1 + 1/144*(-a + 1),
((((((1/172800*(43*a - 76)*r1 + 1/34560*(7*a - 12))*t + (1/86400*(-17*a + 29)*r1 + 1/34560*(-10*a + 17)))*b2^2 + ((1/34560*(7*a - 13)*r1 + 1/34560*(-37*a + 63))*t + (1/34560*(5*a - 8)*r1 + 1/34560*(16*a - 29)))*b2 + ((1/8640*(-193*a + 334)*r1 + 1/8640*(97*a - 168))*t + (1/17280*(176*a - 305)*r1 + 1/17280*(148*a - 257))))*b3 + (((1/57600*(49*a - 85)*r1 + 1/11520*(17*a - 29))*t + (1/86400*(-43*a + 75)*r1 + 1/34560*(-67*a + 117)))*b2^2 + ((1/17280*(-113*a + 198)*r1 + 1/5760*(-11*a + 18))*t + (1/34560*(89*a - 153)*r1 + 1/34560*(383*a - 669)))*b2 + ((1/8640*(-19*a + 33)*r1 + 1/2880*(-97*a + 167))*t + (1/5760*(23*a - 41)*r1 + 1/5760*(147*a - 257)))))*b4 + ((((1/172800*(-3*a - 2)*r1 + 1/34560*(a - 2))*t + (1/172800*(-a + 1)*r1 - 1/34560*a))*b2^2 + ((1/34560*(3*a - 5)*r1 + 1/34560*(-a - 1))*t + (1/34560*(-a + 4)*r1 + 1/34560*(-2*a + 3)))*b2 + ((1/8640*(-69*a + 128)*r1 + 1/8640*(55*a - 92))*t + (1/8640*(28*a - 49)*r1 + 1/2160*(5*a - 6))))*b3 + (((1/172800*(87*a - 151)*r1 + 1/34560*(31*a - 55))*t + (1/86400*(-28*a + 49)*r1 + 1/34560*(-41*a + 69)))*b2^2 + ((1/17280*(-51*a + 88)*r1 + 1/17280*(-19*a + 32))*t + (1/34560*(39*a - 67)*r1 + 1/11520*(63*a - 109)))*b2 + ((1/8640*(-51*a + 79)*r1 + 1/8640*(-183*a + 347))*t + (1/8640*(43*a - 77)*r1 + 1/4320*(97*a - 153))))))*b1 + ((((1/14400*(-a + 3)*r1*t + (1/57600*(4*a - 9)*r1 - 1/11520))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(8*a - 15)*r1*t + (1/576*(2*a - 3)*r1 + 1/576*(-14*a + 25))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(29*a - 51)*r1*t + (1/8640*(43*a - 75)*r1 + 1/1728*(-73*a + 129)))))*b4 + (((1/14400*(-a + 3)*r1*t + (1/57600*(4*a - 9)*r1 - 1/11520))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(5*a - 9)*r1*t + (1/2880*(4*a - 3)*r1 + 1/576*(-8*a + 13))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/1080*(19*a - 30)*r1*t + (1/8640*(97*a - 147)*r1 + 1/1728*(-55*a + 93)))))))*i + (((((1/172800*(26*a - 47)*r1 + 1/34560*(-10*a + 17))*t + (1/57600*(-a + 2)*r1 + 1/17280*(a - 2)))*b2^2 + ((1/34560*(-31*a + 55)*r1 + 1/34560*(-11*a + 19))*t + (1/11520*(7*a - 12)*r1 + 1/34560*(28*a - 47)))*b2 + ((1/8640*(94*a - 163)*r1 + 1/8640*(236*a - 409))*t + (1/1152*(-15*a + 26)*r1 + 1/17280*(-401*a + 694))))*b3 + (((1/57600*(89*a - 155)*r1 + 1/34560*(11*a - 17))*t + (1/172800*(-91*a + 155)*r1 + 1/8640*(-25*a + 44)))*b2^2 + ((1/17280*(26*a - 45)*r1 + 1/17280*(94*a - 161))*t + (1/6912*(-11*a + 19)*r1 + 1/34560*(-139*a + 241)))*b2 + ((1/8640*(-335*a + 579)*r1 + 1/8640*(-195*a + 331))*t + (1/17280*(311*a - 535)*r1 + 1/17280*(1275*a - 2231)))))*b4 + ((((1/172800*(4*a - 9)*r1 - 1/34560)*t + (1/172800*(-2*a + 7)*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/34560*(a + 3)*r1 + 1/34560*(-3*a + 7))*t + (1/34560*(a - 2)*r1 + 1/34560*(4*a - 1)))*b2 + ((1/8640*(50*a - 87)*r1 + 1/8640*(90*a - 151))*t + (1/4320*(-26*a + 43)*r1 + 1/8640*(-85*a + 148))))*b3 + (((1/172800*(155*a - 269)*r1 + 1/34560*(3*a - 5))*t + (1/172800*(-55*a + 97)*r1 + 1/17280*(-26*a + 45)))*b2^2 + ((1/17280*(12*a - 19)*r1 + 1/17280*(40*a - 71))*t + (1/34560*(-21*a + 37)*r1 + 1/34560*(-65*a + 111)))*b2 + ((1/8640*(-199*a + 341)*r1 + 1/8640*(-67*a + 121))*t + (1/4320*(43*a - 77)*r1 + 1/8640*(341*a - 573))))))*b1 + (((1/43200*(-a + 3)*r1*t + (1/172800*(-11*a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/6912*(-a + 3)*r1 + 1/6912*(-a + 3)))*b2 + (1/432*(-7*a + 12)*r1*t + (1/8640*(155*a - 258)*r1 + 1/1728*(-5*a + 6))))*b3 + ((1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(-23*a + 39)*r1*t + (1/2880*(55*a - 97)*r1 + 1/576*(3*a - 5)))))*b4 + ((1/43200*(-a + 3)*r1*t + (1/172800*(-11*a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/6912*(-a + 3)*r1 + 1/6912*(-a + 3)))*b2 + (1/2160*(-17*a + 33)*r1*t + (1/8640*(83*a - 150)*r1 + 1/1728*(-5*a + 6))))*b3 + (1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(-7*a + 12)*r1*t + 1/2880*(37*a - 73)*r1 + 1/576*(-3*a + 7),
((((1/28800*(-a + 4)*r1 + 1/5760*(a - 4))*b2^2 + (1/5760*(-a + 1)*r1 + 1/1152*(a - 1))*b2 + (1/720*(11*a - 20)*r1 + 1/144*(-11*a + 20)))*b3 + ((1/28800*(5*a - 9)*r1 + 1/5760*(-5*a + 9))*b2^2 + (1/2880*(-2*a + 3)*r1 + 1/576*(2*a - 3))*b2 + (1/720*(-22*a + 39)*r1 + 1/144*(22*a - 39))))*b4 + (((1/28800*(-a + 4)*r1 + 1/5760*(a - 4))*b2^2 + (1/5760*(-a + 1)*r1 + 1/1152*(a - 1))*b2 + (1/1440*(13*a - 22)*r1 + 1/288*(-13*a + 22)))*b3 + ((1/28800*(5*a - 9)*r1 + 1/5760*(-5*a + 9))*b2^2 + (1/2880*(-2*a + 3)*r1 + 1/576*(2*a - 3))*b2 + (1/1440*(-17*a + 33)*r1 + 1/288*(17*a - 33)))))*i + (((1/28800*(-2*a + 3)*r1 + 1/5760*(2*a - 3))*b2^2 + (1/5760*(a - 3)*r1 + 1/1152*(-a + 3))*b2 + (1/720*(-5*a + 9)*r1 + 1/144*(5*a - 9)))*b3 + ((1/28800*(11*a - 17)*r1 + 1/5760*(-11*a + 17))*b2^2 + (1/2880*(-a + 2)*r1 + 1/576*(a - 2))*b2 + (1/720*(6*a - 11)*r1 + 1/144*(-6*a + 11))))*b4 + ((1/28800*(-2*a + 3)*r1 + 1/5760*(2*a - 3))*b2^2 + (1/5760*(a - 3)*r1 + 1/1152*(-a + 3))*b2 + (1/1440*(-4*a + 9)*r1 + 1/288*(4*a - 9)))*b3 + (1/28800*(11*a - 17)*r1 + 1/5760*(-11*a + 17))*b2^2 + (1/2880*(-a + 2)*r1 + 1/576*(a - 2))*b2 + 1/1440*(-3*a + 5)*r1 + 1/288*(3*a - 5),
(((((1/86400*(205*a - 351)*r1 + 1/17280*(19*a - 33))*t + (1/172800*(-185*a + 321)*r1 + 1/34560*(-137*a + 237)))*b2^2 + ((1/2880*(-7*a + 12)*r1 + 1/8640*(13*a - 24))*t + (1/5760*(4*a - 7)*r1 + 1/17280*(62*a - 111)))*b2 + ((1/2160*(-100*a + 171)*r1 + 1/1080*(-2*a + 3))*t + (1/4320*(68*a - 117)*r1 + 1/480*(31*a - 54))))*b4 + (((1/86400*(81*a - 139)*r1 + 1/17280*(7*a - 13))*t + (1/172800*(-61*a + 109)*r1 + 1/34560*(-61*a + 105)))*b2^2 + ((1/8640*(-9*a + 16)*r1 + 1/8640*(17*a - 28))*t + (-1/17280*r1 + 1/1920*(2*a - 3)))*b2 + ((1/4320*(-81*a + 145)*r1 + 1/4320*(3*a + 5))*t + (1/8640*(41*a - 73)*r1 + 1/8640*(277*a - 465)))))*b1 + (((1/10800*(-5*a + 9)*r1*t + (1/86400*(-5*a + 9)*r1 + 1/17280*(5*a - 9)))*b2^2 + (1/1080*(2*a - 3)*r1*t + (1/8640*(2*a - 3)*r1 + 1/1728*(-2*a + 3)))*b2 + (1/1080*(-20*a + 33)*r1*t + (1/1080*(65*a - 114)*r1 + 1/432*(-22*a + 39))))*b4 + ((1/10800*(-5*a + 9)*r1*t + (1/86400*(-5*a + 9)*r1 + 1/17280*(5*a - 9)))*b2^2 + (1/1080*(2*a - 3)*r1*t + (1/8640*(2*a - 3)*r1 + 1/1728*(-2*a + 3)))*b2 + (1/1080*(-11*a + 15)*r1*t + (1/4320*(107*a - 195)*r1 + 1/864*(-17*a + 33))))))*i + ((((1/86400*(7*a - 15)*r1 + 1/17280*(-31*a + 55))*t + (1/172800*(53*a - 95)*r1 + 1/34560*(21*a - 35)))*b2^2 + ((1/8640*(40*a - 69)*r1 + 1/8640*(40*a - 71))*t + (1/17280*(-35*a + 62)*r1 + 1/17280*(-185*a + 316)))*b2 + ((1/240*(3*a - 5)*r1 + 1/2160*(72*a - 127))*t + (1/4320*(-67*a + 118)*r1 + 1/4320*(-62*a + 107))))*b4 + (((1/86400*(-5*a + 13)*r1 + 1/17280*(-11*a + 19))*t + (1/172800*(25*a - 43)*r1 + 1/11520*(3*a - 5)))*b2^2 + ((1/8640*(28*a - 49)*r1 + 1/1728*(4*a - 7))*t + (1/17280*(-27*a + 46)*r1 + 1/1152*(-7*a + 12)))*b2 + ((1/4320*(33*a - 67)*r1 + 1/4320*(57*a - 83))*t + (1/8640*(-65*a + 103)*r1 + 1/2880*(-19*a + 41)))))*b1 + ((1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + (1/120*(-7*a + 12)*r1*t + (1/720*(17*a - 28)*r1 + 1/72*(4*a - 7))))*b4 + (1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + 1/120*(-3*a + 5)*r1*t + 1/1440*(31*a - 53)*r1 + 1/288*(7*a - 13),
((((((1/172800*(-26*a + 47)*r1 + 1/34560*(4*a - 7))*t + (1/345600*(11*a - 17)*r1 + 1/23040*(5*a - 9)))*b2^2 + ((1/34560*(-23*a + 41)*r1 + 1/34560*(-31*a + 49))*t + (1/34560*(13*a - 22)*r1 + 1/11520*(16*a - 29)))*b2 + ((1/4320*(115*a - 199)*r1 + 1/4320*(53*a - 92))*t + (1/8640*(-94*a + 163)*r1 + 1/2880*(-141*a + 244))))*b3 + (((1/57600*(-3*a + 5)*r1 + 1/34560*(-81*a + 139))*t + (1/172800*(102*a - 175)*r1 + 1/34560*(47*a - 82)))*b2^2 + ((1/17280*(-83*a + 144)*r1 + 1/17280*(41*a - 68))*t + (1/34560*(83*a - 145)*r1 + 1/34560*(49*a - 77)))*b2 + ((1/4320*(-17*a + 30)*r1 + 1/4320*(176*a - 301))*t + (1/8640*(-91*a + 157)*r1 + 1/8640*(-91*a + 161)))))*b4 + ((((1/172800*(6*a - 1)*r1 - 1/11520)*t + (1/345600*(-a + 11)*r1 + 1/69120*(-5*a + 1)))*b2^2 + ((1/34560*(-3*a + 5)*r1 + 1/11520*(-a - 1))*t + (1/34560*(a - 2)*r1 + 1/34560*(4*a - 11)))*b2 + ((1/8640*(90*a - 167)*r1 + 1/2880*(20*a - 27))*t + (1/17280*(-83*a + 133)*r1 + 1/17280*(-335*a + 631))))*b3 + (((1/172800*(11*a - 13)*r1 + 1/34560*(-37*a + 63))*t + (1/172800*(42*a - 71)*r1 + 1/34560*(19*a - 34)))*b2^2 + ((1/17280*(-51*a + 88)*r1 + 1/17280*(13*a - 24))*t + (1/34560*(51*a - 89)*r1 + 1/34560*(57*a - 101)))*b2 + ((1/8640*(-37*a + 55)*r1 + 1/2880*(57*a - 103))*t + (1/8640*(-36*a + 65)*r1 + 1/8640*(-43*a + 62))))))*b1 + ((((1/14400*(a - 3)*r1*t + (1/57600*(5*a - 12)*r1 + 1/11520*(-a + 2)))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(-8*a + 15)*r1*t + (1/1440*a*r1 + 1/288*(-3*a + 5))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(-29*a + 51)*r1*t + (1/4320*(7*a - 12)*r1 + 1/864*(-22*a + 39)))))*b4 + (((1/14400*(a - 3)*r1*t + (1/57600*(5*a - 12)*r1 + 1/11520*(-a + 2)))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(-5*a + 9)*r1*t + (1/2880*(-a + 6)*r1 + 1/576*(-3*a + 4))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/1080*(-19*a + 30)*r1*t + (1/8640*(59*a - 87)*r1 + 1/1728*(-17*a + 33)))))))*i + (((((1/172800*(-47*a + 84)*r1 + 1/34560*(-9*a + 14))*t + (1/345600*(47*a - 79)*r1 + 1/69120*(37*a - 67)))*b2^2 + ((1/6912*(5*a - 9)*r1 + 1/34560*(-3*a + 7))*t + (1/34560*(-8*a + 13)*r1 + 1/34560*(-41*a + 74)))*b2 + ((1/4320*(7*a - 12)*r1 + 1/2160*(42*a - 73))*t + (1/8640*(-34*a + 59)*r1 + 1/8640*(-79*a + 136))))*b3 + (((1/57600*(117*a - 205)*r1 + 1/34560*(7*a - 11))*t + (1/172800*(-203*a + 350)*r1 + 1/34560*(-58*a + 101)))*b2^2 + ((1/1920*(-4*a + 7)*r1 + 1/17280*(-108*a + 187))*t + (1/34560*(95*a - 163)*r1 + 1/34560*(183*a - 317)))*b2 + ((1/1440*(-57*a + 100)*r1 + 1/4320*(-20*a + 29))*t + (1/1728*(39*a - 67)*r1 + 1/8640*(265*a - 469)))))*b4 + ((((1/172800*(-3*a + 8)*r1 + 1/34560*(-a - 2))*t + (1/115200*(a - 1)*r1 + 1/69120*(a - 7)))*b2^2 + ((1/34560*(-3*a - 1)*r1 + 1/34560*(a + 11))*t + (-1/11520*r1 + 1/34560*(7*a + 2)))*b2 + ((1/8640*(-9*a + 16)*r1 + 1/8640*(65*a - 98))*t + (1/5760*(-5*a + 9)*r1 + 1/17280*(-5*a + 23))))*b3 + (((1/172800*(155*a - 267)*r1 + 1/34560*(3*a - 7))*t + (1/172800*(-95*a + 166)*r1 + 1/34560*(-22*a + 37)))*b2^2 + ((1/17280*(-16*a + 27)*r1 + 1/17280*(-60*a + 107))*t + (1/34560*(47*a - 83)*r1 + 1/34560*(95*a - 157)))*b2 + ((1/2880*(-55*a + 91)*r1 + 1/8640*(-21*a + 53))*t + (1/8640*(99*a - 178)*r1 + 1/8640*(108*a - 179))))))*b1 + (((1/43200*(a - 3)*r1*t + (1/172800*(-10*a + 3)*r1 + 1/34560*(4*a - 3)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/432*(7*a - 12)*r1*t + (1/4320*(95*a - 159)*r1 + 1/864*(-20*a + 33))))*b3 + ((1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(23*a - 39)*r1*t + (1/1440*(39*a - 68)*r1 + 1/288*(-10*a + 17)))))*b4 + ((1/43200*(a - 3)*r1*t + (1/172800*(-10*a + 3)*r1 + 1/34560*(4*a - 3)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/2160*(17*a - 33)*r1*t + (1/8640*(100*a - 183)*r1 + 1/1728*(-22*a + 39))))*b3 + (1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(7*a - 12)*r1*t + 1/2880*(51*a - 97)*r1 + 1/576*(-17*a + 31),
((((((1/172800*(-26*a + 47)*r1 + 1/34560*(4*a - 7))*t + (1/345600*(11*a - 17)*r1 + 1/23040*(5*a - 9)))*b2^2 + ((1/34560*(-23*a + 41)*r1 + 1/34560*(-31*a + 49))*t + (1/34560*(13*a - 22)*r1 + 1/11520*(16*a - 29)))*b2 + ((1/4320*(115*a - 199)*r1 + 1/4320*(53*a - 92))*t + (1/8640*(-94*a + 163)*r1 + 1/2880*(-141*a + 244))))*b3 + (((1/57600*(3*a - 5)*r1 + 1/34560*(81*a - 139))*t + (1/172800*(-102*a + 175)*r1 + 1/34560*(-47*a + 82)))*b2^2 + ((1/17280*(83*a - 144)*r1 + 1/17280*(-41*a + 68))*t + (1/34560*(-83*a + 145)*r1 + 1/34560*(-49*a + 77)))*b2 + ((1/4320*(17*a - 30)*r1 + 1/4320*(-176*a + 301))*t + (1/8640*(91*a - 157)*r1 + 1/8640*(91*a - 161)))))*b4 + ((((1/172800*(6*a - 1)*r1 - 1/11520)*t + (1/345600*(-a + 11)*r1 + 1/69120*(-5*a + 1)))*b2^2 + ((1/34560*(-3*a + 5)*r1 + 1/11520*(-a - 1))*t + (1/34560*(a - 2)*r1 + 1/34560*(4*a - 11)))*b2 + ((1/8640*(90*a - 167)*r1 + 1/2880*(20*a - 27))*t + (1/17280*(-83*a + 133)*r1 + 1/17280*(-335*a + 631))))*b3 + (((1/172800*(-11*a + 13)*r1 + 1/34560*(37*a - 63))*t + (1/172800*(-42*a + 71)*r1 + 1/34560*(-19*a + 34)))*b2^2 + ((1/17280*(51*a - 88)*r1 + 1/17280*(-13*a + 24))*t + (1/34560*(-51*a + 89)*r1 + 1/34560*(-57*a + 101)))*b2 + ((1/8640*(37*a - 55)*r1 + 1/2880*(-57*a + 103))*t + (1/8640*(36*a - 65)*r1 + 1/8640*(43*a - 62))))))*b1 + ((((1/14400*(-a + 3)*r1*t + (1/57600*(-5*a + 12)*r1 + 1/11520*(a - 2)))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(8*a - 15)*r1*t + (-1/1440*a*r1 + 1/288*(3*a - 5))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(-29*a + 51)*r1*t + (1/4320*(7*a - 12)*r1 + 1/864*(-22*a + 39)))))*b4 + (((1/14400*(-a + 3)*r1*t + (1/57600*(-5*a + 12)*r1 + 1/11520*(a - 2)))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(5*a - 9)*r1*t + (1/2880*(a - 6)*r1 + 1/576*(3*a - 4))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/1080*(-19*a + 30)*r1*t + (1/8640*(59*a - 87)*r1 + 1/1728*(-17*a + 33)))))))*i + (((((1/172800*(-47*a + 84)*r1 + 1/34560*(-9*a + 14))*t + (1/345600*(47*a - 79)*r1 + 1/69120*(37*a - 67)))*b2^2 + ((1/6912*(5*a - 9)*r1 + 1/34560*(-3*a + 7))*t + (1/34560*(-8*a + 13)*r1 + 1/34560*(-41*a + 74)))*b2 + ((1/4320*(7*a - 12)*r1 + 1/2160*(42*a - 73))*t + (1/8640*(-34*a + 59)*r1 + 1/8640*(-79*a + 136))))*b3 + (((1/57600*(-117*a + 205)*r1 + 1/34560*(-7*a + 11))*t + (1/172800*(203*a - 350)*r1 + 1/34560*(58*a - 101)))*b2^2 + ((1/1920*(4*a - 7)*r1 + 1/17280*(108*a - 187))*t + (1/34560*(-95*a + 163)*r1 + 1/34560*(-183*a + 317)))*b2 + ((1/1440*(57*a - 100)*r1 + 1/4320*(20*a - 29))*t + (1/1728*(-39*a + 67)*r1 + 1/8640*(-265*a + 469)))))*b4 + ((((1/172800*(-3*a + 8)*r1 + 1/34560*(-a - 2))*t + (1/115200*(a - 1)*r1 + 1/69120*(a - 7)))*b2^2 + ((1/34560*(-3*a - 1)*r1 + 1/34560*(a + 11))*t + (-1/11520*r1 + 1/34560*(7*a + 2)))*b2 + ((1/8640*(-9*a + 16)*r1 + 1/8640*(65*a - 98))*t + (1/5760*(-5*a + 9)*r1 + 1/17280*(-5*a + 23))))*b3 + (((1/172800*(-155*a + 267)*r1 + 1/34560*(-3*a + 7))*t + (1/172800*(95*a - 166)*r1 + 1/34560*(22*a - 37)))*b2^2 + ((1/17280*(16*a - 27)*r1 + 1/17280*(60*a - 107))*t + (1/34560*(-47*a + 83)*r1 + 1/34560*(-95*a + 157)))*b2 + ((1/2880*(55*a - 91)*r1 + 1/8640*(21*a - 53))*t + (1/8640*(-99*a + 178)*r1 + 1/8640*(-108*a + 179))))))*b1 + (((1/43200*(-a + 3)*r1*t + (1/172800*(10*a - 3)*r1 + 1/34560*(-4*a + 3)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/432*(-7*a + 12)*r1*t + (1/4320*(-95*a + 159)*r1 + 1/864*(20*a - 33))))*b3 + ((1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(23*a - 39)*r1*t + (1/1440*(39*a - 68)*r1 + 1/288*(-10*a + 17)))))*b4 + ((1/43200*(-a + 3)*r1*t + (1/172800*(10*a - 3)*r1 + 1/34560*(-4*a + 3)))*b2^2 + (1/4320*(a - 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/2160*(-17*a + 33)*r1*t + (1/8640*(-100*a + 183)*r1 + 1/1728*(22*a - 39))))*b3 + (1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(7*a - 12)*r1*t + 1/2880*(51*a - 97)*r1 + 1/576*(-17*a + 31),
(((((1/86400*(-205*a + 351)*r1 + 1/17280*(-19*a + 33))*t + (1/172800*(185*a - 321)*r1 + 1/34560*(137*a - 237)))*b2^2 + ((1/2880*(7*a - 12)*r1 + 1/8640*(-13*a + 24))*t + (1/5760*(-4*a + 7)*r1 + 1/17280*(-62*a + 111)))*b2 + ((1/2160*(100*a - 171)*r1 + 1/1080*(2*a - 3))*t + (1/4320*(-68*a + 117)*r1 + 1/480*(-31*a + 54))))*b4 + (((1/86400*(-81*a + 139)*r1 + 1/17280*(-7*a + 13))*t + (1/172800*(61*a - 109)*r1 + 1/34560*(61*a - 105)))*b2^2 + ((1/8640*(9*a - 16)*r1 + 1/8640*(-17*a + 28))*t + (1/17280*r1 + 1/1920*(-2*a + 3)))*b2 + ((1/4320*(81*a - 145)*r1 + 1/4320*(-3*a - 5))*t + (1/8640*(-41*a + 73)*r1 + 1/8640*(-277*a + 465)))))*b1 + (((1/10800*(-5*a + 9)*r1*t + (1/86400*(-5*a + 9)*r1 + 1/17280*(5*a - 9)))*b2^2 + (1/1080*(2*a - 3)*r1*t + (1/8640*(2*a - 3)*r1 + 1/1728*(-2*a + 3)))*b2 + (1/1080*(-20*a + 33)*r1*t + (1/1080*(65*a - 114)*r1 + 1/432*(-22*a + 39))))*b4 + ((1/10800*(-5*a + 9)*r1*t + (1/86400*(-5*a + 9)*r1 + 1/17280*(5*a - 9)))*b2^2 + (1/1080*(2*a - 3)*r1*t + (1/8640*(2*a - 3)*r1 + 1/1728*(-2*a + 3)))*b2 + (1/1080*(-11*a + 15)*r1*t + (1/4320*(107*a - 195)*r1 + 1/864*(-17*a + 33))))))*i + ((((1/86400*(-7*a + 15)*r1 + 1/17280*(31*a - 55))*t + (1/172800*(-53*a + 95)*r1 + 1/34560*(-21*a + 35)))*b2^2 + ((1/8640*(-40*a + 69)*r1 + 1/8640*(-40*a + 71))*t + (1/17280*(35*a - 62)*r1 + 1/17280*(185*a - 316)))*b2 + ((1/240*(-3*a + 5)*r1 + 1/2160*(-72*a + 127))*t + (1/4320*(67*a - 118)*r1 + 1/4320*(62*a - 107))))*b4 + (((1/86400*(5*a - 13)*r1 + 1/17280*(11*a - 19))*t + (1/172800*(-25*a + 43)*r1 + 1/11520*(-3*a + 5)))*b2^2 + ((1/8640*(-28*a + 49)*r1 + 1/1728*(-4*a + 7))*t + (1/17280*(27*a - 46)*r1 + 1/1152*(7*a - 12)))*b2 + ((1/4320*(-33*a + 67)*r1 + 1/4320*(-57*a + 83))*t + (1/8640*(65*a - 103)*r1 + 1/2880*(19*a - 41)))))*b1 + ((1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + (1/120*(-7*a + 12)*r1*t + (1/720*(17*a - 28)*r1 + 1/72*(4*a - 7))))*b4 + (1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + 1/120*(-3*a + 5)*r1*t + 1/1440*(31*a - 53)*r1 + 1/288*(7*a - 13),
((((1/28800*(a - 4)*r1 + 1/5760*(-a + 4))*b2^2 + (1/5760*(a - 1)*r1 + 1/1152*(-a + 1))*b2 + (1/720*(-11*a + 20)*r1 + 1/144*(11*a - 20)))*b3 + ((1/28800*(5*a - 9)*r1 + 1/5760*(-5*a + 9))*b2^2 + (1/2880*(-2*a + 3)*r1 + 1/576*(2*a - 3))*b2 + (1/720*(-22*a + 39)*r1 + 1/144*(22*a - 39))))*b4 + (((1/28800*(a - 4)*r1 + 1/5760*(-a + 4))*b2^2 + (1/5760*(a - 1)*r1 + 1/1152*(-a + 1))*b2 + (1/1440*(-13*a + 22)*r1 + 1/288*(13*a - 22)))*b3 + ((1/28800*(5*a - 9)*r1 + 1/5760*(-5*a + 9))*b2^2 + (1/2880*(-2*a + 3)*r1 + 1/576*(2*a - 3))*b2 + (1/1440*(-17*a + 33)*r1 + 1/288*(17*a - 33)))))*i + (((1/28800*(2*a - 3)*r1 + 1/5760*(-2*a + 3))*b2^2 + (1/5760*(-a + 3)*r1 + 1/1152*(a - 3))*b2 + (1/720*(5*a - 9)*r1 + 1/144*(-5*a + 9)))*b3 + ((1/28800*(11*a - 17)*r1 + 1/5760*(-11*a + 17))*b2^2 + (1/2880*(-a + 2)*r1 + 1/576*(a - 2))*b2 + (1/720*(6*a - 11)*r1 + 1/144*(-6*a + 11))))*b4 + ((1/28800*(2*a - 3)*r1 + 1/5760*(-2*a + 3))*b2^2 + (1/5760*(-a + 3)*r1 + 1/1152*(a - 3))*b2 + (1/1440*(4*a - 9)*r1 + 1/288*(-4*a + 9)))*b3 + (1/28800*(11*a - 17)*r1 + 1/5760*(-11*a + 17))*b2^2 + (1/2880*(-a + 2)*r1 + 1/576*(a - 2))*b2 + 1/1440*(-3*a + 5)*r1 + 1/288*(3*a - 5),
((((((1/172800*(43*a - 76)*r1 + 1/34560*(7*a - 12))*t + (1/86400*(-17*a + 29)*r1 + 1/34560*(-10*a + 17)))*b2^2 + ((1/34560*(7*a - 13)*r1 + 1/34560*(-37*a + 63))*t + (1/34560*(5*a - 8)*r1 + 1/34560*(16*a - 29)))*b2 + ((1/8640*(-193*a + 334)*r1 + 1/8640*(97*a - 168))*t + (1/17280*(176*a - 305)*r1 + 1/17280*(148*a - 257))))*b3 + (((1/57600*(-49*a + 85)*r1 + 1/11520*(-17*a + 29))*t + (1/86400*(43*a - 75)*r1 + 1/34560*(67*a - 117)))*b2^2 + ((1/17280*(113*a - 198)*r1 + 1/5760*(11*a - 18))*t + (1/34560*(-89*a + 153)*r1 + 1/34560*(-383*a + 669)))*b2 + ((1/8640*(19*a - 33)*r1 + 1/2880*(97*a - 167))*t + (1/5760*(-23*a + 41)*r1 + 1/5760*(-147*a + 257)))))*b4 + ((((1/172800*(-3*a - 2)*r1 + 1/34560*(a - 2))*t + (1/172800*(-a + 1)*r1 - 1/34560*a))*b2^2 + ((1/34560*(3*a - 5)*r1 + 1/34560*(-a - 1))*t + (1/34560*(-a + 4)*r1 + 1/34560*(-2*a + 3)))*b2 + ((1/8640*(-69*a + 128)*r1 + 1/8640*(55*a - 92))*t + (1/8640*(28*a - 49)*r1 + 1/2160*(5*a - 6))))*b3 + (((1/172800*(-87*a + 151)*r1 + 1/34560*(-31*a + 55))*t + (1/86400*(28*a - 49)*r1 + 1/34560*(41*a - 69)))*b2^2 + ((1/17280*(51*a - 88)*r1 + 1/17280*(19*a - 32))*t + (1/34560*(-39*a + 67)*r1 + 1/11520*(-63*a + 109)))*b2 + ((1/8640*(51*a - 79)*r1 + 1/8640*(183*a - 347))*t + (1/8640*(-43*a + 77)*r1 + 1/4320*(-97*a + 153))))))*b1 + ((((1/14400*(a - 3)*r1*t + (1/57600*(-4*a + 9)*r1 + 1/11520))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(-8*a + 15)*r1*t + (1/576*(-2*a + 3)*r1 + 1/576*(14*a - 25))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(29*a - 51)*r1*t + (1/8640*(43*a - 75)*r1 + 1/1728*(-73*a + 129)))))*b4 + (((1/14400*(a - 3)*r1*t + (1/57600*(-4*a + 9)*r1 + 1/11520))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(-5*a + 9)*r1*t + (1/2880*(-4*a + 3)*r1 + 1/576*(8*a - 13))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/1080*(19*a - 30)*r1*t + (1/8640*(97*a - 147)*r1 + 1/1728*(-55*a + 93)))))))*i + (((((1/172800*(26*a - 47)*r1 + 1/34560*(-10*a + 17))*t + (1/57600*(-a + 2)*r1 + 1/17280*(a - 2)))*b2^2 + ((1/34560*(-31*a + 55)*r1 + 1/34560*(-11*a + 19))*t + (1/11520*(7*a - 12)*r1 + 1/34560*(28*a - 47)))*b2 + ((1/8640*(94*a - 163)*r1 + 1/8640*(236*a - 409))*t + (1/1152*(-15*a + 26)*r1 + 1/17280*(-401*a + 694))))*b3 + (((1/57600*(-89*a + 155)*r1 + 1/34560*(-11*a + 17))*t + (1/172800*(91*a - 155)*r1 + 1/8640*(25*a - 44)))*b2^2 + ((1/17280*(-26*a + 45)*r1 + 1/17280*(-94*a + 161))*t + (1/6912*(11*a - 19)*r1 + 1/34560*(139*a - 241)))*b2 + ((1/8640*(335*a - 579)*r1 + 1/8640*(195*a - 331))*t + (1/17280*(-311*a + 535)*r1 + 1/17280*(-1275*a + 2231)))))*b4 + ((((1/172800*(4*a - 9)*r1 - 1/34560)*t + (1/172800*(-2*a + 7)*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/34560*(a + 3)*r1 + 1/34560*(-3*a + 7))*t + (1/34560*(a - 2)*r1 + 1/34560*(4*a - 1)))*b2 + ((1/8640*(50*a - 87)*r1 + 1/8640*(90*a - 151))*t + (1/4320*(-26*a + 43)*r1 + 1/8640*(-85*a + 148))))*b3 + (((1/172800*(-155*a + 269)*r1 + 1/34560*(-3*a + 5))*t + (1/172800*(55*a - 97)*r1 + 1/17280*(26*a - 45)))*b2^2 + ((1/17280*(-12*a + 19)*r1 + 1/17280*(-40*a + 71))*t + (1/34560*(21*a - 37)*r1 + 1/34560*(65*a - 111)))*b2 + ((1/8640*(199*a - 341)*r1 + 1/8640*(67*a - 121))*t + (1/4320*(-43*a + 77)*r1 + 1/8640*(-341*a + 573))))))*b1 + (((1/43200*(a - 3)*r1*t + (1/172800*(11*a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/6912*(a - 3)*r1 + 1/6912*(a - 3)))*b2 + (1/432*(7*a - 12)*r1*t + (1/8640*(-155*a + 258)*r1 + 1/1728*(5*a - 6))))*b3 + ((1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(-23*a + 39)*r1*t + (1/2880*(55*a - 97)*r1 + 1/576*(3*a - 5)))))*b4 + ((1/43200*(a - 3)*r1*t + (1/172800*(11*a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/6912*(a - 3)*r1 + 1/6912*(a - 3)))*b2 + (1/2160*(17*a - 33)*r1*t + (1/8640*(-83*a + 150)*r1 + 1/1728*(5*a - 6))))*b3 + (1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(-7*a + 12)*r1*t + 1/2880*(37*a - 73)*r1 + 1/576*(-3*a + 7),
(((((1/28800*(-35*a + 61)*r1 + 1/17280*(17*a - 31))*t + (1/172800*(95*a - 163)*r1 + 1/6912*(-a + 1)))*b2^2 + ((1/8640*(41*a - 72)*r1 + 1/8640*(51*a - 88))*t + (1/17280*(-70*a + 121)*r1 + 1/17280*(-116*a + 203)))*b2 + ((1/4320*(43*a - 78)*r1 + 1/4320*(-191*a + 344))*t + (1/8640*(-6*a + 7)*r1 + 1/8640*(286*a - 469))))*b4 + (((1/86400*(-73*a + 127)*r1 + 1/5760*(3*a - 5))*t + (1/172800*(63*a - 107)*r1 + 1/34560*(11*a - 19)))*b2^2 + ((1/8640*(17*a - 28)*r1 + 1/8640*(19*a - 36))*t + (1/17280*(-30*a + 53)*r1 + 1/17280*(-44*a + 71)))*b2 + ((1/216*(3*a - 5)*r1 + 1/2160*(-47*a + 75))*t + (1/4320*(-18*a + 31)*r1 + 1/4320*(12*a - 35)))))*b1 + (((1/10800*(5*a - 9)*r1*t + (1/86400*(5*a - 9)*r1 + 1/17280*(-5*a + 9)))*b2^2 + (1/1080*(-2*a + 3)*r1*t + (1/8640*(-2*a + 3)*r1 + 1/1728*(2*a - 3)))*b2 + (1/1080*(20*a - 33)*r1*t + (1/4320*(280*a - 489)*r1 + 1/864*(-64*a + 111))))*b4 + ((1/10800*(5*a - 9)*r1*t + (1/86400*(5*a - 9)*r1 + 1/17280*(-5*a + 9)))*b2^2 + (1/1080*(-2*a + 3)*r1*t + (1/8640*(-2*a + 3)*r1 + 1/1728*(2*a - 3)))*b2 + (1/1080*(11*a - 15)*r1*t + (1/2160*(59*a - 105)*r1 + 1/216*(-7*a + 12))))))*i + ((((1/86400*(-61*a + 105)*r1 + 1/17280*(-27*a + 47))*t + (1/172800*(151*a - 265)*r1 + 1/34560*(43*a - 73)))*b2^2 + ((1/8640*(-44*a + 75)*r1 + 1/8640*(40*a - 67))*t + (1/17280*(25*a - 44)*r1 + 1/17280*(35*a - 58)))*b2 + ((1/4320*(166*a - 285)*r1 + 1/4320*(128*a - 223))*t + (1/8640*(-225*a + 394)*r1 + 1/8640*(-513*a + 878))))*b4 + (((1/9600*(-5*a + 9)*r1 + 1/17280*(-19*a + 31))*t + (1/172800*(95*a - 161)*r1 + 1/6912*(7*a - 13)))*b2^2 + ((1/8640*(-16*a + 27)*r1 + 1/1728*(4*a - 7))*t + (1/17280*(9*a - 16)*r1 + 1/3456*(-a + 2)))*b2 + ((1/2160*(47*a - 81)*r1 + 1/540*(13*a - 22))*t + (1/4320*(-68*a + 115)*r1 + 1/4320*(-162*a + 293)))))*b1 + ((1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + (1/120*(7*a - 12)*r1*t + (1/1440*(55*a - 92)*r1 + 1/288*(-5*a + 8))))*b4 + (1/28800*(-19*a + 29)*r1 + 1/5760*(3*a - 5))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + 1/120*(3*a - 5)*r1*t + 1/360*(10*a - 17)*r1 + 1/144*(-a + 1),
((((((1/172800*(26*a - 47)*r1 + 1/34560*(-4*a + 7))*t + (1/345600*(-11*a + 17)*r1 + 1/23040*(-5*a + 9)))*b2^2 + ((1/34560*(23*a - 41)*r1 + 1/34560*(31*a - 49))*t + (1/34560*(-13*a + 22)*r1 + 1/11520*(-16*a + 29)))*b2 + ((1/4320*(-115*a + 199)*r1 + 1/4320*(-53*a + 92))*t + (1/8640*(94*a - 163)*r1 + 1/2880*(141*a - 244))))*b3 + (((1/57600*(3*a - 5)*r1 + 1/34560*(81*a - 139))*t + (1/172800*(-102*a + 175)*r1 + 1/34560*(-47*a + 82)))*b2^2 + ((1/17280*(83*a - 144)*r1 + 1/17280*(-41*a + 68))*t + (1/34560*(-83*a + 145)*r1 + 1/34560*(-49*a + 77)))*b2 + ((1/4320*(17*a - 30)*r1 + 1/4320*(-176*a + 301))*t + (1/8640*(91*a - 157)*r1 + 1/8640*(91*a - 161)))))*b4 + ((((1/172800*(-6*a + 1)*r1 + 1/11520)*t + (1/345600*(a - 11)*r1 + 1/69120*(5*a - 1)))*b2^2 + ((1/34560*(3*a - 5)*r1 + 1/11520*(a + 1))*t + (1/34560*(-a + 2)*r1 + 1/34560*(-4*a + 11)))*b2 + ((1/8640*(-90*a + 167)*r1 + 1/2880*(-20*a + 27))*t + (1/17280*(83*a - 133)*r1 + 1/17280*(335*a - 631))))*b3 + (((1/172800*(-11*a + 13)*r1 + 1/34560*(37*a - 63))*t + (1/172800*(-42*a + 71)*r1 + 1/34560*(-19*a + 34)))*b2^2 + ((1/17280*(51*a - 88)*r1 + 1/17280*(-13*a + 24))*t + (1/34560*(-51*a + 89)*r1 + 1/34560*(-57*a + 101)))*b2 + ((1/8640*(37*a - 55)*r1 + 1/2880*(-57*a + 103))*t + (1/8640*(36*a - 65)*r1 + 1/8640*(43*a - 62))))))*b1 + ((((1/14400*(a - 3)*r1*t + (1/57600*(5*a - 12)*r1 + 1/11520*(-a + 2)))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(-8*a + 15)*r1*t + (1/1440*a*r1 + 1/288*(-3*a + 5))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(-29*a + 51)*r1*t + (1/4320*(7*a - 12)*r1 + 1/864*(-22*a + 39)))))*b4 + (((1/14400*(a - 3)*r1*t + (1/57600*(5*a - 12)*r1 + 1/11520*(-a + 2)))*b2^2 + (1/3840*(-a + 1)*r1 + 1/2304*(a - 1))*b2 + (1/720*(-5*a + 9)*r1*t + (1/2880*(-a + 6)*r1 + 1/576*(-3*a + 4))))*b3 + ((1/21600*(4*a - 9)*r1*t + (1/172800*(-11*a + 27)*r1 + 1/34560*(5*a - 9)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/1080*(-19*a + 30)*r1*t + (1/8640*(59*a - 87)*r1 + 1/1728*(-17*a + 33)))))))*i + (((((1/172800*(47*a - 84)*r1 + 1/34560*(9*a - 14))*t + (1/345600*(-47*a + 79)*r1 + 1/69120*(-37*a + 67)))*b2^2 + ((1/6912*(-5*a + 9)*r1 + 1/34560*(3*a - 7))*t + (1/34560*(8*a - 13)*r1 + 1/34560*(41*a - 74)))*b2 + ((1/4320*(-7*a + 12)*r1 + 1/2160*(-42*a + 73))*t + (1/8640*(34*a - 59)*r1 + 1/8640*(79*a - 136))))*b3 + (((1/57600*(-117*a + 205)*r1 + 1/34560*(-7*a + 11))*t + (1/172800*(203*a - 350)*r1 + 1/34560*(58*a - 101)))*b2^2 + ((1/1920*(4*a - 7)*r1 + 1/17280*(108*a - 187))*t + (1/34560*(-95*a + 163)*r1 + 1/34560*(-183*a + 317)))*b2 + ((1/1440*(57*a - 100)*r1 + 1/4320*(20*a - 29))*t + (1/1728*(-39*a + 67)*r1 + 1/8640*(-265*a + 469)))))*b4 + ((((1/172800*(3*a - 8)*r1 + 1/34560*(a + 2))*t + (1/115200*(-a + 1)*r1 + 1/69120*(-a + 7)))*b2^2 + ((1/34560*(3*a + 1)*r1 + 1/34560*(-a - 11))*t + (1/11520*r1 + 1/34560*(-7*a - 2)))*b2 + ((1/8640*(9*a - 16)*r1 + 1/8640*(-65*a + 98))*t + (1/5760*(5*a - 9)*r1 + 1/17280*(5*a - 23))))*b3 + (((1/172800*(-155*a + 267)*r1 + 1/34560*(-3*a + 7))*t + (1/172800*(95*a - 166)*r1 + 1/34560*(22*a - 37)))*b2^2 + ((1/17280*(16*a - 27)*r1 + 1/17280*(60*a - 107))*t + (1/34560*(-47*a + 83)*r1 + 1/34560*(-95*a + 157)))*b2 + ((1/2880*(55*a - 91)*r1 + 1/8640*(21*a - 53))*t + (1/8640*(-99*a + 178)*r1 + 1/8640*(-108*a + 179))))))*b1 + (((1/43200*(a - 3)*r1*t + (1/172800*(-10*a + 3)*r1 + 1/34560*(4*a - 3)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/432*(7*a - 12)*r1*t + (1/4320*(95*a - 159)*r1 + 1/864*(-20*a + 33))))*b3 + ((1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(23*a - 39)*r1*t + (1/1440*(39*a - 68)*r1 + 1/288*(-10*a + 17)))))*b4 + ((1/43200*(a - 3)*r1*t + (1/172800*(-10*a + 3)*r1 + 1/34560*(4*a - 3)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/2160*(17*a - 33)*r1*t + (1/8640*(100*a - 183)*r1 + 1/1728*(-22*a + 39))))*b3 + (1/7200*(-2*a + 3)*r1*t + (1/57600*(5*a - 11)*r1 + 1/11520*(5*a - 7)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(7*a - 12)*r1*t + 1/2880*(51*a - 97)*r1 + 1/576*(-17*a + 31),
((((((1/172800*(-43*a + 76)*r1 + 1/34560*(-7*a + 12))*t + (1/86400*(17*a - 29)*r1 + 1/34560*(10*a - 17)))*b2^2 + ((1/34560*(-7*a + 13)*r1 + 1/34560*(37*a - 63))*t + (1/34560*(-5*a + 8)*r1 + 1/34560*(-16*a + 29)))*b2 + ((1/8640*(193*a - 334)*r1 + 1/8640*(-97*a + 168))*t + (1/17280*(-176*a + 305)*r1 + 1/17280*(-148*a + 257))))*b3 + (((1/57600*(49*a - 85)*r1 + 1/11520*(17*a - 29))*t + (1/86400*(-43*a + 75)*r1 + 1/34560*(-67*a + 117)))*b2^2 + ((1/17280*(-113*a + 198)*r1 + 1/5760*(-11*a + 18))*t + (1/34560*(89*a - 153)*r1 + 1/34560*(383*a - 669)))*b2 + ((1/8640*(-19*a + 33)*r1 + 1/2880*(-97*a + 167))*t + (1/5760*(23*a - 41)*r1 + 1/5760*(147*a - 257)))))*b4 + ((((1/172800*(3*a + 2)*r1 + 1/34560*(-a + 2))*t + (1/172800*(a - 1)*r1 + 1/34560*a))*b2^2 + ((1/34560*(-3*a + 5)*r1 + 1/34560*(a + 1))*t + (1/34560*(a - 4)*r1 + 1/34560*(2*a - 3)))*b2 + ((1/8640*(69*a - 128)*r1 + 1/8640*(-55*a + 92))*t + (1/8640*(-28*a + 49)*r1 + 1/2160*(-5*a + 6))))*b3 + (((1/172800*(87*a - 151)*r1 + 1/34560*(31*a - 55))*t + (1/86400*(-28*a + 49)*r1 + 1/34560*(-41*a + 69)))*b2^2 + ((1/17280*(-51*a + 88)*r1 + 1/17280*(-19*a + 32))*t + (1/34560*(39*a - 67)*r1 + 1/11520*(63*a - 109)))*b2 + ((1/8640*(-51*a + 79)*r1 + 1/8640*(-183*a + 347))*t + (1/8640*(43*a - 77)*r1 + 1/4320*(97*a - 153))))))*b1 + ((((1/14400*(a - 3)*r1*t + (1/57600*(-4*a + 9)*r1 + 1/11520))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(-8*a + 15)*r1*t + (1/576*(-2*a + 3)*r1 + 1/576*(14*a - 25))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(29*a - 51)*r1*t + (1/8640*(43*a - 75)*r1 + 1/1728*(-73*a + 129)))))*b4 + (((1/14400*(a - 3)*r1*t + (1/57600*(-4*a + 9)*r1 + 1/11520))*b2^2 + (1/3840*(a - 1)*r1 + 1/2304*(-a + 1))*b2 + (1/720*(-5*a + 9)*r1*t + (1/2880*(-4*a + 3)*r1 + 1/576*(8*a - 13))))*b3 + ((1/21600*(-4*a + 9)*r1*t + (1/172800*(-19*a + 45)*r1 + 1/34560*(13*a - 27)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/1080*(19*a - 30)*r1*t + (1/8640*(97*a - 147)*r1 + 1/1728*(-55*a + 93)))))))*i + (((((1/172800*(-26*a + 47)*r1 + 1/34560*(10*a - 17))*t + (1/57600*(a - 2)*r1 + 1/17280*(-a + 2)))*b2^2 + ((1/34560*(31*a - 55)*r1 + 1/34560*(11*a - 19))*t + (1/11520*(-7*a + 12)*r1 + 1/34560*(-28*a + 47)))*b2 + ((1/8640*(-94*a + 163)*r1 + 1/8640*(-236*a + 409))*t + (1/1152*(15*a - 26)*r1 + 1/17280*(401*a - 694))))*b3 + (((1/57600*(89*a - 155)*r1 + 1/34560*(11*a - 17))*t + (1/172800*(-91*a + 155)*r1 + 1/8640*(-25*a + 44)))*b2^2 + ((1/17280*(26*a - 45)*r1 + 1/17280*(94*a - 161))*t + (1/6912*(-11*a + 19)*r1 + 1/34560*(-139*a + 241)))*b2 + ((1/8640*(-335*a + 579)*r1 + 1/8640*(-195*a + 331))*t + (1/17280*(311*a - 535)*r1 + 1/17280*(1275*a - 2231)))))*b4 + ((((1/172800*(-4*a + 9)*r1 + 1/34560)*t + (1/172800*(2*a - 7)*r1 + 1/34560*(a - 1)))*b2^2 + ((1/34560*(-a - 3)*r1 + 1/34560*(3*a - 7))*t + (1/34560*(-a + 2)*r1 + 1/34560*(-4*a + 1)))*b2 + ((1/8640*(-50*a + 87)*r1 + 1/8640*(-90*a + 151))*t + (1/4320*(26*a - 43)*r1 + 1/8640*(85*a - 148))))*b3 + (((1/172800*(155*a - 269)*r1 + 1/34560*(3*a - 5))*t + (1/172800*(-55*a + 97)*r1 + 1/17280*(-26*a + 45)))*b2^2 + ((1/17280*(12*a - 19)*r1 + 1/17280*(40*a - 71))*t + (1/34560*(-21*a + 37)*r1 + 1/34560*(-65*a + 111)))*b2 + ((1/8640*(-199*a + 341)*r1 + 1/8640*(-67*a + 121))*t + (1/4320*(43*a - 77)*r1 + 1/8640*(341*a - 573))))))*b1 + (((1/43200*(a - 3)*r1*t + (1/172800*(11*a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/6912*(a - 3)*r1 + 1/6912*(a - 3)))*b2 + (1/432*(7*a - 12)*r1*t + (1/8640*(-155*a + 258)*r1 + 1/1728*(5*a - 6))))*b3 + ((1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + (1/720*(-23*a + 39)*r1*t + (1/2880*(55*a - 97)*r1 + 1/576*(3*a - 5)))))*b4 + ((1/43200*(a - 3)*r1*t + (1/172800*(11*a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/4320*(-a + 3)*r1*t + (1/6912*(a - 3)*r1 + 1/6912*(a - 3)))*b2 + (1/2160*(17*a - 33)*r1*t + (1/8640*(-83*a + 150)*r1 + 1/1728*(5*a - 6))))*b3 + (1/7200*(2*a - 3)*r1*t + (1/57600*(9*a - 17)*r1 + 1/11520*(a - 1)))*b2^2 + (1/5760*(-7*a + 14)*r1 + 1/1152*(a - 2))*b2 + 1/360*(-7*a + 12)*r1*t + 1/2880*(37*a - 73)*r1 + 1/576*(-3*a + 7)
];