d15fidb:=[
(1/3600*(a - 2)*r1 + 1/720*(-a + 2))*b2^2 + (1/720*(-a + 1)*r1 + 1/144*(a - 1))*b2 + 1/360*(-6*a + 7)*r1 + 1/72*(6*a - 7),
(((((1/172800*(5*a - 9)*r1 + 1/34560*(-3*a + 7))*t + (1/86400*r1 + 1/34560*(a - 1)))*b2^2 + ((1/17280*(-2*a + 3)*r1 + 1/17280*(-6*a + 7))*t + (1/34560*(3*a - 7)*r1 + 1/34560*(11*a - 17)))*b2 + ((1/8640*(-11*a + 15)*r1 + 1/8640*(9*a - 41))*t + (1/17280*(3*a - 19)*r1 + 1/17280*(11*a + 1))))*b3 + (((1/28800*(-a + 6)*r1 + 1/17280*(a - 6))*t + (1/43200*(2*a + 3)*r1 - 1/5760))*b2^2 + ((1/17280*(-a - 3)*r1 + 1/17280*(11*a - 15))*t + (-1/17280*a*r1 + 1/17280*(-8*a + 15)))*b2 + ((1/4320*(5*a - 21)*r1 + 1/864*(-a + 9))*t + (-1/576*r1 + 1/1728*(-2*a + 9)))))*b1 + (((1/7200*(2*a - 3)*r1*t + (1/57600*(-7*a + 9)*r1 + 1/11520*(-3*a + 5)))*b2^2 + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2))*b2 + (1/720*(-5*a + 9)*r1*t + (1/2880*(5*a - 3)*r1 + 1/576*(5*a - 11))))*b3 + ((1/21600*(7*a - 9)*r1*t + (1/86400*(20*a - 27)*r1 + 1/17280*(-8*a + 9)))*b2^2 + (1/2160*(-a + 3)*r1*t + (1/17280*(-7*a + 21)*r1 + 1/3456*(a - 3)))*b2 + (1/1080*(-13*a + 24)*r1*t + (1/4320*(-41*a + 75)*r1 + 1/864*(11*a - 21))))))*i + ((((1/172800*(5*a - 7)*r1 + 1/34560*(5*a - 7))*t + (1/172800*(-5*a + 11)*r1 + 1/17280*(-2*a + 3)))*b2^2 + ((1/17280*(a - 2)*r1 + 1/17280*(-7*a + 14))*t + (1/34560*(a - 1)*r1 + 1/34560*(7*a - 9)))*b2 + ((1/8640*(-5*a + 7)*r1 + 1/8640*(-37*a + 83))*t + (1/17280*(19*a - 25)*r1 + 1/17280*(67*a - 93))))*b3 + (((1/86400*(-4*a + 3)*r1 + 1/17280*(-8*a + 11))*t + (1/86400*(7*a - 14)*r1 + 1/4320*(a - 2)))*b2^2 + ((1/17280*(-5*a + 3)*r1 + 1/17280*(7*a + 7))*t + (1/17280*(a - 8)*r1 + 1/17280*(4*a - 1)))*b2 + ((1/1440*(a - 1)*r1 + 1/4320*(57*a - 73))*t + (1/8640*(-17*a + 28)*r1 + 1/8640*(-57*a + 98)))))*b1 + ((1/21600*(-2*a + 3)*r1*t + (1/172800*(37*a - 69)*r1 + 1/34560*(-7*a + 15)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(13*a - 15)*r1*t + (1/8640*(-73*a + 105)*r1 + 1/1728*(7*a - 15))))*b3 + (1/7200*(-a + 3)*r1*t + (1/28800*(3*a - 2)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(2*a - 3)*r1*t + 1/1440*(-13*a + 11)*r1 + 1/288*(3*a - 5),
(((((1/172800*(5*a - 7)*r1 + 1/34560*(5*a - 11))*t + (1/172800*(-10*a + 21)*r1 + 1/34560*(-3*a + 4)))*b2^2 + ((1/17280*(-2*a + 1)*r1 + 1/17280*(2*a + 1))*t + (1/34560*(-a - 3)*r1 + 1/34560*(3*a - 1)))*b2 + ((1/4320*(-7*a + 2)*r1 + 1/4320*(-14*a + 41))*t + (1/4320*(5*a - 18)*r1 + 1/4320*(12*a - 13))))*b3 + (((1/86400*(4*a - 9)*r1 + 1/17280*(2*a + 1))*t + (1/172800*(-9*a + 19)*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/5760*(-a + 1)*r1 + 1/17280*(17*a - 25))*t + (1/17280*(-4*a + 5)*r1 + 1/17280*(-a + 10)))*b2 + ((1/4320*(-11*a + 12)*r1 + 1/864*(3*a - 4))*t + (1/8640*(a - 25)*r1 + 1/1728*(a + 5)))))*b1 + (((1/7200*(2*a - 3)*r1*t + (1/57600*(11*a - 15)*r1 + 1/11520*(-a + 1)))*b2^2 + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2))*b2 + (1/720*(-5*a + 9)*r1*t + (1/1440*(-5*a + 6)*r1 + 1/288)))*b3 + ((1/21600*(-7*a + 9)*r1*t + (1/86400*(13*a - 18)*r1 - 1/17280*a))*b2^2 + (1/2160*(a - 3)*r1*t + (1/3456*(-a + 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/1080*(13*a - 24)*r1*t + (1/4320*(-28*a + 51)*r1 + 1/864*(-2*a + 3))))))*i + ((((1/172800*(5*a - 1)*r1 + 1/11520*(-a + 1))*t + (1/172800*(5*a - 2)*r1 - 1/34560))*b2^2 + ((1/17280*(a - 2)*r1 + 1/1920*(a - 2))*t + (1/34560*(-7*a + 13)*r1 + 1/34560*(-9*a + 13)))*b2 + ((1/4320*(-a + 2)*r1 + 1/1440*(4*a - 11))*t + (1/2160*(-2*a + 5)*r1 + 1/2160*(-3*a + 2))))*b3 + (((1/28800*(-a + 2)*r1 + 1/17280*(-5*a + 8))*t + (1/172800*(23*a - 31)*r1 + 1/34560*(5*a - 7)))*b2^2 + ((1/5760*(a - 3)*r1 + 1/17280*(19*a - 41))*t + (1/17280*(-9*a + 10)*r1 + 1/17280*(-12*a + 23)))*b2 + ((1/1440*(2*a - 7)*r1 + 1/4320*(54*a - 91))*t + (1/8640*(-49*a + 61)*r1 + 1/8640*(-59*a + 101)))))*b1 + ((1/21600*(-2*a + 3)*r1*t + (1/172800*(-41*a + 75)*r1 + 1/34560*(11*a - 21)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/17280*(14*a - 21)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(13*a - 15)*r1*t + (1/4320*(43*a - 60)*r1 + 1/864*(-10*a + 15))))*b3 + (1/7200*(a - 3)*r1*t + (1/28800*(4*a - 5)*r1 + 1/5760*(-2*a + 5)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(-2*a + 3)*r1*t + 1/1440*(-15*a + 14)*r1 + 1/288*(5*a - 8),
((((1/43200*(-2*a + 9)*r1 + 1/8640*(2*a - 5))*t + (1/86400*(-3*a + 16)*r1 + 1/17280*(-a - 2)))*b2^2 + ((1/2880*(-a + 1)*r1 + 1/8640*(-13*a + 29))*t + (1/17280*(13*a - 23)*r1 + 1/17280*(13*a - 19)))*b2 + ((1/4320*(2*a - 21)*r1 + 1/4320*(-48*a + 89))*t + (1/8640*(29*a - 68)*r1 + 1/8640*(13*a - 4))))*b1 + ((-1/5400*a*r1*t + (-1/43200*a*r1 + 1/8640*a))*b2^2 + (1/1080*(-a + 3)*r1*t + (1/8640*(-a + 3)*r1 + 1/1728*(a - 3)))*b2 + (1/1080*(a + 12)*r1*t + (1/4320*(23*a + 6)*r1 + 1/864*(-5*a - 6)))))*i + (((1/43200*(5*a - 6)*r1 + 1/8640*(7*a - 10))*t + (1/86400*(-20*a + 41)*r1 + 1/17280*(-6*a + 11)))*b2^2 + ((1/8640*(-a + 3)*r1 + 1/8640*(a - 11))*t + (1/17280*(-a + 5)*r1 + 1/17280*(-a + 1)))*b2 + ((1/4320*(-7*a + 30)*r1 + 1/4320*(-105*a + 112))*t + (1/8640*(66*a - 79)*r1 + 1/8640*(70*a - 167))))*b1 + (1/14400*(4*a - 13)*r1 + 1/2880)*b2^2 + (1/2880*(a - 1)*r1 + 1/576*(-a + 1))*b2 + -1/120*r1*t + 1/1440*(-14*a + 29)*r1 + 1/288*(-2*a + 1),
(((((1/172800*(-5*a + 9)*r1 + 1/34560*(3*a - 7))*t + (-1/86400*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/17280*(2*a - 3)*r1 + 1/17280*(6*a - 7))*t + (1/34560*(-3*a + 7)*r1 + 1/34560*(-11*a + 17)))*b2 + ((1/8640*(11*a - 15)*r1 + 1/8640*(-9*a + 41))*t + (1/17280*(-3*a + 19)*r1 + 1/17280*(-11*a - 1))))*b3 + (((1/28800*(a - 6)*r1 + 1/17280*(-a + 6))*t + (1/43200*(-2*a - 3)*r1 + 1/5760))*b2^2 + ((1/17280*(a + 3)*r1 + 1/17280*(-11*a + 15))*t + (1/17280*a*r1 + 1/17280*(8*a - 15)))*b2 + ((1/4320*(-5*a + 21)*r1 + 1/864*(a - 9))*t + (1/576*r1 + 1/1728*(2*a - 9)))))*b1 + (((1/7200*(2*a - 3)*r1*t + (1/57600*(-7*a + 9)*r1 + 1/11520*(-3*a + 5)))*b2^2 + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2))*b2 + (1/720*(-5*a + 9)*r1*t + (1/2880*(5*a - 3)*r1 + 1/576*(5*a - 11))))*b3 + ((1/21600*(7*a - 9)*r1*t + (1/86400*(20*a - 27)*r1 + 1/17280*(-8*a + 9)))*b2^2 + (1/2160*(-a + 3)*r1*t + (1/17280*(-7*a + 21)*r1 + 1/3456*(a - 3)))*b2 + (1/1080*(-13*a + 24)*r1*t + (1/4320*(-41*a + 75)*r1 + 1/864*(11*a - 21))))))*i + ((((1/172800*(-5*a + 7)*r1 + 1/34560*(-5*a + 7))*t + (1/172800*(5*a - 11)*r1 + 1/17280*(2*a - 3)))*b2^2 + ((1/17280*(-a + 2)*r1 + 1/17280*(7*a - 14))*t + (1/34560*(-a + 1)*r1 + 1/34560*(-7*a + 9)))*b2 + ((1/8640*(5*a - 7)*r1 + 1/8640*(37*a - 83))*t + (1/17280*(-19*a + 25)*r1 + 1/17280*(-67*a + 93))))*b3 + (((1/86400*(4*a - 3)*r1 + 1/17280*(8*a - 11))*t + (1/86400*(-7*a + 14)*r1 + 1/4320*(-a + 2)))*b2^2 + ((1/17280*(5*a - 3)*r1 + 1/17280*(-7*a - 7))*t + (1/17280*(-a + 8)*r1 + 1/17280*(-4*a + 1)))*b2 + ((1/1440*(-a + 1)*r1 + 1/4320*(-57*a + 73))*t + (1/8640*(17*a - 28)*r1 + 1/8640*(57*a - 98)))))*b1 + ((1/21600*(-2*a + 3)*r1*t + (1/172800*(37*a - 69)*r1 + 1/34560*(-7*a + 15)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/3456*(-2*a + 3)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(13*a - 15)*r1*t + (1/8640*(-73*a + 105)*r1 + 1/1728*(7*a - 15))))*b3 + (1/7200*(-a + 3)*r1*t + (1/28800*(3*a - 2)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(2*a - 3)*r1*t + 1/1440*(-13*a + 11)*r1 + 1/288*(3*a - 5),
(((1/28800*(7*a - 11)*r1 + 1/5760*(-7*a + 11))*b2^2 + (1/2880*(-a + 2)*r1 + 1/576*(a - 2))*b2 + (1/144*(-a + 2)*r1 + 1/144*(5*a - 10)))*b3 + ((-1/14400*a*r1 + 1/2880*a)*b2^2 + (1/2880*(-a + 3)*r1 + 1/576*(a - 3))*b2 + (1/720*(-2*a + 3)*r1 + 1/144*(2*a - 3))))*i + ((1/28800*(a - 3)*r1 + 1/5760*(-a + 3))*b2^2 + (1/2880*(2*a - 3)*r1 + 1/576*(-2*a + 3))*b2 + (1/720*a*r1 - 1/144*a))*b3 + (1/14400*(-2*a + 7)*r1 + 1/2880*(2*a - 7))*b2^2 + (1/2880*(-a + 1)*r1 + 1/576*(a - 1))*b2 + 1/720*(3*a - 10)*r1 + 1/144*(-3*a + 10),
((((1/43200*(-4*a + 3)*r1 + 1/8640*(-4*a + 9))*t + (1/9600*(a - 2)*r1 + 1/17280*(7*a - 12)))*b2^2 + ((1/2880*(-a + 1)*r1 + 1/8640*(19*a - 27))*t + (1/17280*(-5*a + 3)*r1 + 1/1920*(-a + 3)))*b2 + ((1/2160*(2*a - 3)*r1 + 1/2160*(47*a - 81))*t + (1/4320*(-17*a + 30)*r1 + 1/4320*(-52*a + 111))))*b1 + ((1/5400*a*r1*t + (1/43200*a*r1 - 1/8640*a))*b2^2 + (1/1080*(a - 3)*r1*t + (1/8640*(a - 3)*r1 + 1/1728*(-a + 3)))*b2 + (1/1080*(-a - 12)*r1*t + (1/2160*(11*a - 3)*r1 + 1/432*(-2*a + 3)))))*i + (((1/43200*(5*a - 12)*r1 + 1/8640*(-a + 4))*t + (1/86400*(-10*a + 7)*r1 + 1/17280))*b2^2 + ((1/2880*(-a + 3)*r1 + 1/8640*(-13*a + 23))*t + (1/17280*(11*a - 11)*r1 + 1/17280*(23*a - 43)))*b2 + ((1/1080*(-5*a + 12)*r1 + 7/2160)*t + (1/4320*(25*a - 22)*r1 + 1/4320*(20*a - 47))))*b1 + (1/14400*(4*a - 13)*r1 + 1/2880)*b2^2 + (1/2880*(a - 1)*r1 + 1/576*(-a + 1))*b2 + 1/120*r1*t + 1/720*(-7*a + 16)*r1 + 1/144*(-a - 1),
(((((1/172800*(-5*a + 7)*r1 + 1/34560*(-5*a + 11))*t + (1/172800*(10*a - 21)*r1 + 1/34560*(3*a - 4)))*b2^2 + ((1/17280*(2*a - 1)*r1 + 1/17280*(-2*a - 1))*t + (1/34560*(a + 3)*r1 + 1/34560*(-3*a + 1)))*b2 + ((1/4320*(7*a - 2)*r1 + 1/4320*(14*a - 41))*t + (1/4320*(-5*a + 18)*r1 + 1/4320*(-12*a + 13))))*b3 + (((1/86400*(4*a - 9)*r1 + 1/17280*(2*a + 1))*t + (1/172800*(-9*a + 19)*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/5760*(-a + 1)*r1 + 1/17280*(17*a - 25))*t + (1/17280*(-4*a + 5)*r1 + 1/17280*(-a + 10)))*b2 + ((1/4320*(-11*a + 12)*r1 + 1/864*(3*a - 4))*t + (1/8640*(a - 25)*r1 + 1/1728*(a + 5)))))*b1 + (((1/7200*(-2*a + 3)*r1*t + (1/57600*(-11*a + 15)*r1 + 1/11520*(a - 1)))*b2^2 + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2))*b2 + (1/720*(5*a - 9)*r1*t + (1/1440*(5*a - 6)*r1 - 1/288)))*b3 + ((1/21600*(-7*a + 9)*r1*t + (1/86400*(13*a - 18)*r1 - 1/17280*a))*b2^2 + (1/2160*(a - 3)*r1*t + (1/3456*(-a + 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/1080*(13*a - 24)*r1*t + (1/4320*(-28*a + 51)*r1 + 1/864*(-2*a + 3))))))*i + ((((1/172800*(-5*a + 1)*r1 + 1/11520*(a - 1))*t + (1/172800*(-5*a + 2)*r1 + 1/34560))*b2^2 + ((1/17280*(-a + 2)*r1 + 1/1920*(-a + 2))*t + (1/34560*(7*a - 13)*r1 + 1/34560*(9*a - 13)))*b2 + ((1/4320*(a - 2)*r1 + 1/1440*(-4*a + 11))*t + (1/2160*(2*a - 5)*r1 + 1/2160*(3*a - 2))))*b3 + (((1/28800*(-a + 2)*r1 + 1/17280*(-5*a + 8))*t + (1/172800*(23*a - 31)*r1 + 1/34560*(5*a - 7)))*b2^2 + ((1/5760*(a - 3)*r1 + 1/17280*(19*a - 41))*t + (1/17280*(-9*a + 10)*r1 + 1/17280*(-12*a + 23)))*b2 + ((1/1440*(2*a - 7)*r1 + 1/4320*(54*a - 91))*t + (1/8640*(-49*a + 61)*r1 + 1/8640*(-59*a + 101)))))*b1 + ((1/21600*(2*a - 3)*r1*t + (1/172800*(41*a - 75)*r1 + 1/34560*(-11*a + 21)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(-13*a + 15)*r1*t + (1/4320*(-43*a + 60)*r1 + 1/864*(10*a - 15))))*b3 + (1/7200*(a - 3)*r1*t + (1/28800*(4*a - 5)*r1 + 1/5760*(-2*a + 5)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(-2*a + 3)*r1*t + 1/1440*(-15*a + 14)*r1 + 1/288*(5*a - 8),
(((((1/172800*(-5*a + 7)*r1 + 1/34560*(-5*a + 11))*t + (1/172800*(10*a - 21)*r1 + 1/34560*(3*a - 4)))*b2^2 + ((1/17280*(2*a - 1)*r1 + 1/17280*(-2*a - 1))*t + (1/34560*(a + 3)*r1 + 1/34560*(-3*a + 1)))*b2 + ((1/4320*(7*a - 2)*r1 + 1/4320*(14*a - 41))*t + (1/4320*(-5*a + 18)*r1 + 1/4320*(-12*a + 13))))*b3 + (((1/86400*(-4*a + 9)*r1 + 1/17280*(-2*a - 1))*t + (1/172800*(9*a - 19)*r1 + 1/34560*(a - 1)))*b2^2 + ((1/5760*(a - 1)*r1 + 1/17280*(-17*a + 25))*t + (1/17280*(4*a - 5)*r1 + 1/17280*(a - 10)))*b2 + ((1/4320*(11*a - 12)*r1 + 1/864*(-3*a + 4))*t + (1/8640*(-a + 25)*r1 + 1/1728*(-a - 5)))))*b1 + (((1/7200*(2*a - 3)*r1*t + (1/57600*(11*a - 15)*r1 + 1/11520*(-a + 1)))*b2^2 + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2))*b2 + (1/720*(-5*a + 9)*r1*t + (1/1440*(-5*a + 6)*r1 + 1/288)))*b3 + ((1/21600*(-7*a + 9)*r1*t + (1/86400*(13*a - 18)*r1 - 1/17280*a))*b2^2 + (1/2160*(a - 3)*r1*t + (1/3456*(-a + 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/1080*(13*a - 24)*r1*t + (1/4320*(-28*a + 51)*r1 + 1/864*(-2*a + 3))))))*i + ((((1/172800*(-5*a + 1)*r1 + 1/11520*(a - 1))*t + (1/172800*(-5*a + 2)*r1 + 1/34560))*b2^2 + ((1/17280*(-a + 2)*r1 + 1/1920*(-a + 2))*t + (1/34560*(7*a - 13)*r1 + 1/34560*(9*a - 13)))*b2 + ((1/4320*(a - 2)*r1 + 1/1440*(-4*a + 11))*t + (1/2160*(2*a - 5)*r1 + 1/2160*(3*a - 2))))*b3 + (((1/28800*(a - 2)*r1 + 1/17280*(5*a - 8))*t + (1/172800*(-23*a + 31)*r1 + 1/34560*(-5*a + 7)))*b2^2 + ((1/5760*(-a + 3)*r1 + 1/17280*(-19*a + 41))*t + (1/17280*(9*a - 10)*r1 + 1/17280*(12*a - 23)))*b2 + ((1/1440*(-2*a + 7)*r1 + 1/4320*(-54*a + 91))*t + (1/8640*(49*a - 61)*r1 + 1/8640*(59*a - 101)))))*b1 + ((1/21600*(-2*a + 3)*r1*t + (1/172800*(-41*a + 75)*r1 + 1/34560*(11*a - 21)))*b2^2 + (1/2160*(2*a - 3)*r1*t + (1/17280*(14*a - 21)*r1 + 1/3456*(-2*a + 3)))*b2 + (1/2160*(13*a - 15)*r1*t + (1/4320*(43*a - 60)*r1 + 1/864*(-10*a + 15))))*b3 + (1/7200*(a - 3)*r1*t + (1/28800*(4*a - 5)*r1 + 1/5760*(-2*a + 5)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(-2*a + 3)*r1*t + 1/1440*(-15*a + 14)*r1 + 1/288*(5*a - 8),
((((1/43200*(4*a - 3)*r1 + 1/8640*(4*a - 9))*t + (1/9600*(-a + 2)*r1 + 1/17280*(-7*a + 12)))*b2^2 + ((1/2880*(a - 1)*r1 + 1/8640*(-19*a + 27))*t + (1/17280*(5*a - 3)*r1 + 1/1920*(a - 3)))*b2 + ((1/2160*(-2*a + 3)*r1 + 1/2160*(-47*a + 81))*t + (1/4320*(17*a - 30)*r1 + 1/4320*(52*a - 111))))*b1 + ((1/5400*a*r1*t + (1/43200*a*r1 - 1/8640*a))*b2^2 + (1/1080*(a - 3)*r1*t + (1/8640*(a - 3)*r1 + 1/1728*(-a + 3)))*b2 + (1/1080*(-a - 12)*r1*t + (1/2160*(11*a - 3)*r1 + 1/432*(-2*a + 3)))))*i + (((1/43200*(-5*a + 12)*r1 + 1/8640*(a - 4))*t + (1/86400*(10*a - 7)*r1 - 1/17280))*b2^2 + ((1/2880*(a - 3)*r1 + 1/8640*(13*a - 23))*t + (1/17280*(-11*a + 11)*r1 + 1/17280*(-23*a + 43)))*b2 + ((1/1080*(5*a - 12)*r1 - 7/2160)*t + (1/4320*(-25*a + 22)*r1 + 1/4320*(-20*a + 47))))*b1 + (1/14400*(4*a - 13)*r1 + 1/2880)*b2^2 + (1/2880*(a - 1)*r1 + 1/576*(-a + 1))*b2 + 1/120*r1*t + 1/720*(-7*a + 16)*r1 + 1/144*(-a - 1),
(((1/28800*(-7*a + 11)*r1 + 1/5760*(7*a - 11))*b2^2 + (1/2880*(a - 2)*r1 + 1/576*(-a + 2))*b2 + (1/144*(a - 2)*r1 + 1/144*(-5*a + 10)))*b3 + ((-1/14400*a*r1 + 1/2880*a)*b2^2 + (1/2880*(-a + 3)*r1 + 1/576*(a - 3))*b2 + (1/720*(-2*a + 3)*r1 + 1/144*(2*a - 3))))*i + ((1/28800*(-a + 3)*r1 + 1/5760*(a - 3))*b2^2 + (1/2880*(-2*a + 3)*r1 + 1/576*(2*a - 3))*b2 + (-1/720*a*r1 + 1/144*a))*b3 + (1/14400*(-2*a + 7)*r1 + 1/2880*(2*a - 7))*b2^2 + (1/2880*(-a + 1)*r1 + 1/576*(a - 1))*b2 + 1/720*(3*a - 10)*r1 + 1/144*(-3*a + 10),
(((((1/172800*(-5*a + 9)*r1 + 1/34560*(3*a - 7))*t + (-1/86400*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/17280*(2*a - 3)*r1 + 1/17280*(6*a - 7))*t + (1/34560*(-3*a + 7)*r1 + 1/34560*(-11*a + 17)))*b2 + ((1/8640*(11*a - 15)*r1 + 1/8640*(-9*a + 41))*t + (1/17280*(-3*a + 19)*r1 + 1/17280*(-11*a - 1))))*b3 + (((1/28800*(-a + 6)*r1 + 1/17280*(a - 6))*t + (1/43200*(2*a + 3)*r1 - 1/5760))*b2^2 + ((1/17280*(-a - 3)*r1 + 1/17280*(11*a - 15))*t + (-1/17280*a*r1 + 1/17280*(-8*a + 15)))*b2 + ((1/4320*(5*a - 21)*r1 + 1/864*(-a + 9))*t + (-1/576*r1 + 1/1728*(-2*a + 9)))))*b1 + (((1/7200*(-2*a + 3)*r1*t + (1/57600*(7*a - 9)*r1 + 1/11520*(3*a - 5)))*b2^2 + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2))*b2 + (1/720*(5*a - 9)*r1*t + (1/2880*(-5*a + 3)*r1 + 1/576*(-5*a + 11))))*b3 + ((1/21600*(7*a - 9)*r1*t + (1/86400*(20*a - 27)*r1 + 1/17280*(-8*a + 9)))*b2^2 + (1/2160*(-a + 3)*r1*t + (1/17280*(-7*a + 21)*r1 + 1/3456*(a - 3)))*b2 + (1/1080*(-13*a + 24)*r1*t + (1/4320*(-41*a + 75)*r1 + 1/864*(11*a - 21))))))*i + ((((1/172800*(-5*a + 7)*r1 + 1/34560*(-5*a + 7))*t + (1/172800*(5*a - 11)*r1 + 1/17280*(2*a - 3)))*b2^2 + ((1/17280*(-a + 2)*r1 + 1/17280*(7*a - 14))*t + (1/34560*(-a + 1)*r1 + 1/34560*(-7*a + 9)))*b2 + ((1/8640*(5*a - 7)*r1 + 1/8640*(37*a - 83))*t + (1/17280*(-19*a + 25)*r1 + 1/17280*(-67*a + 93))))*b3 + (((1/86400*(-4*a + 3)*r1 + 1/17280*(-8*a + 11))*t + (1/86400*(7*a - 14)*r1 + 1/4320*(a - 2)))*b2^2 + ((1/17280*(-5*a + 3)*r1 + 1/17280*(7*a + 7))*t + (1/17280*(a - 8)*r1 + 1/17280*(4*a - 1)))*b2 + ((1/1440*(a - 1)*r1 + 1/4320*(57*a - 73))*t + (1/8640*(-17*a + 28)*r1 + 1/8640*(-57*a + 98)))))*b1 + ((1/21600*(2*a - 3)*r1*t + (1/172800*(-37*a + 69)*r1 + 1/34560*(7*a - 15)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/3456*(2*a - 3)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(-13*a + 15)*r1*t + (1/8640*(73*a - 105)*r1 + 1/1728*(-7*a + 15))))*b3 + (1/7200*(-a + 3)*r1*t + (1/28800*(3*a - 2)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(2*a - 3)*r1*t + 1/1440*(-13*a + 11)*r1 + 1/288*(3*a - 5),
((((1/43200*(2*a - 9)*r1 + 1/8640*(-2*a + 5))*t + (1/86400*(3*a - 16)*r1 + 1/17280*(a + 2)))*b2^2 + ((1/2880*(a - 1)*r1 + 1/8640*(13*a - 29))*t + (1/17280*(-13*a + 23)*r1 + 1/17280*(-13*a + 19)))*b2 + ((1/4320*(-2*a + 21)*r1 + 1/4320*(48*a - 89))*t + (1/8640*(-29*a + 68)*r1 + 1/8640*(-13*a + 4))))*b1 + ((-1/5400*a*r1*t + (-1/43200*a*r1 + 1/8640*a))*b2^2 + (1/1080*(-a + 3)*r1*t + (1/8640*(-a + 3)*r1 + 1/1728*(a - 3)))*b2 + (1/1080*(a + 12)*r1*t + (1/4320*(23*a + 6)*r1 + 1/864*(-5*a - 6)))))*i + (((1/43200*(-5*a + 6)*r1 + 1/8640*(-7*a + 10))*t + (1/86400*(20*a - 41)*r1 + 1/17280*(6*a - 11)))*b2^2 + ((1/8640*(a - 3)*r1 + 1/8640*(-a + 11))*t + (1/17280*(a - 5)*r1 + 1/17280*(a - 1)))*b2 + ((1/4320*(7*a - 30)*r1 + 1/4320*(105*a - 112))*t + (1/8640*(-66*a + 79)*r1 + 1/8640*(-70*a + 167))))*b1 + (1/14400*(4*a - 13)*r1 + 1/2880)*b2^2 + (1/2880*(a - 1)*r1 + 1/576*(-a + 1))*b2 + -1/120*r1*t + 1/1440*(-14*a + 29)*r1 + 1/288*(-2*a + 1),
(((((1/172800*(5*a - 7)*r1 + 1/34560*(5*a - 11))*t + (1/172800*(-10*a + 21)*r1 + 1/34560*(-3*a + 4)))*b2^2 + ((1/17280*(-2*a + 1)*r1 + 1/17280*(2*a + 1))*t + (1/34560*(-a - 3)*r1 + 1/34560*(3*a - 1)))*b2 + ((1/4320*(-7*a + 2)*r1 + 1/4320*(-14*a + 41))*t + (1/4320*(5*a - 18)*r1 + 1/4320*(12*a - 13))))*b3 + (((1/86400*(-4*a + 9)*r1 + 1/17280*(-2*a - 1))*t + (1/172800*(9*a - 19)*r1 + 1/34560*(a - 1)))*b2^2 + ((1/5760*(a - 1)*r1 + 1/17280*(-17*a + 25))*t + (1/17280*(4*a - 5)*r1 + 1/17280*(a - 10)))*b2 + ((1/4320*(11*a - 12)*r1 + 1/864*(-3*a + 4))*t + (1/8640*(-a + 25)*r1 + 1/1728*(-a - 5)))))*b1 + (((1/7200*(-2*a + 3)*r1*t + (1/57600*(-11*a + 15)*r1 + 1/11520*(a - 1)))*b2^2 + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2))*b2 + (1/720*(5*a - 9)*r1*t + (1/1440*(5*a - 6)*r1 - 1/288)))*b3 + ((1/21600*(-7*a + 9)*r1*t + (1/86400*(13*a - 18)*r1 - 1/17280*a))*b2^2 + (1/2160*(a - 3)*r1*t + (1/3456*(-a + 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/1080*(13*a - 24)*r1*t + (1/4320*(-28*a + 51)*r1 + 1/864*(-2*a + 3))))))*i + ((((1/172800*(5*a - 1)*r1 + 1/11520*(-a + 1))*t + (1/172800*(5*a - 2)*r1 - 1/34560))*b2^2 + ((1/17280*(a - 2)*r1 + 1/1920*(a - 2))*t + (1/34560*(-7*a + 13)*r1 + 1/34560*(-9*a + 13)))*b2 + ((1/4320*(-a + 2)*r1 + 1/1440*(4*a - 11))*t + (1/2160*(-2*a + 5)*r1 + 1/2160*(-3*a + 2))))*b3 + (((1/28800*(a - 2)*r1 + 1/17280*(5*a - 8))*t + (1/172800*(-23*a + 31)*r1 + 1/34560*(-5*a + 7)))*b2^2 + ((1/5760*(-a + 3)*r1 + 1/17280*(-19*a + 41))*t + (1/17280*(9*a - 10)*r1 + 1/17280*(12*a - 23)))*b2 + ((1/1440*(-2*a + 7)*r1 + 1/4320*(-54*a + 91))*t + (1/8640*(49*a - 61)*r1 + 1/8640*(59*a - 101)))))*b1 + ((1/21600*(2*a - 3)*r1*t + (1/172800*(41*a - 75)*r1 + 1/34560*(-11*a + 21)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/17280*(-14*a + 21)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(-13*a + 15)*r1*t + (1/4320*(-43*a + 60)*r1 + 1/864*(10*a - 15))))*b3 + (1/7200*(a - 3)*r1*t + (1/28800*(4*a - 5)*r1 + 1/5760*(-2*a + 5)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(-2*a + 3)*r1*t + 1/1440*(-15*a + 14)*r1 + 1/288*(5*a - 8),
(((((1/172800*(5*a - 9)*r1 + 1/34560*(-3*a + 7))*t + (1/86400*r1 + 1/34560*(a - 1)))*b2^2 + ((1/17280*(-2*a + 3)*r1 + 1/17280*(-6*a + 7))*t + (1/34560*(3*a - 7)*r1 + 1/34560*(11*a - 17)))*b2 + ((1/8640*(-11*a + 15)*r1 + 1/8640*(9*a - 41))*t + (1/17280*(3*a - 19)*r1 + 1/17280*(11*a + 1))))*b3 + (((1/28800*(a - 6)*r1 + 1/17280*(-a + 6))*t + (1/43200*(-2*a - 3)*r1 + 1/5760))*b2^2 + ((1/17280*(a + 3)*r1 + 1/17280*(-11*a + 15))*t + (1/17280*a*r1 + 1/17280*(8*a - 15)))*b2 + ((1/4320*(-5*a + 21)*r1 + 1/864*(a - 9))*t + (1/576*r1 + 1/1728*(2*a - 9)))))*b1 + (((1/7200*(-2*a + 3)*r1*t + (1/57600*(7*a - 9)*r1 + 1/11520*(3*a - 5)))*b2^2 + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2))*b2 + (1/720*(5*a - 9)*r1*t + (1/2880*(-5*a + 3)*r1 + 1/576*(-5*a + 11))))*b3 + ((1/21600*(7*a - 9)*r1*t + (1/86400*(20*a - 27)*r1 + 1/17280*(-8*a + 9)))*b2^2 + (1/2160*(-a + 3)*r1*t + (1/17280*(-7*a + 21)*r1 + 1/3456*(a - 3)))*b2 + (1/1080*(-13*a + 24)*r1*t + (1/4320*(-41*a + 75)*r1 + 1/864*(11*a - 21))))))*i + ((((1/172800*(5*a - 7)*r1 + 1/34560*(5*a - 7))*t + (1/172800*(-5*a + 11)*r1 + 1/17280*(-2*a + 3)))*b2^2 + ((1/17280*(a - 2)*r1 + 1/17280*(-7*a + 14))*t + (1/34560*(a - 1)*r1 + 1/34560*(7*a - 9)))*b2 + ((1/8640*(-5*a + 7)*r1 + 1/8640*(-37*a + 83))*t + (1/17280*(19*a - 25)*r1 + 1/17280*(67*a - 93))))*b3 + (((1/86400*(4*a - 3)*r1 + 1/17280*(8*a - 11))*t + (1/86400*(-7*a + 14)*r1 + 1/4320*(-a + 2)))*b2^2 + ((1/17280*(5*a - 3)*r1 + 1/17280*(-7*a - 7))*t + (1/17280*(-a + 8)*r1 + 1/17280*(-4*a + 1)))*b2 + ((1/1440*(-a + 1)*r1 + 1/4320*(-57*a + 73))*t + (1/8640*(17*a - 28)*r1 + 1/8640*(57*a - 98)))))*b1 + ((1/21600*(2*a - 3)*r1*t + (1/172800*(-37*a + 69)*r1 + 1/34560*(7*a - 15)))*b2^2 + (1/2160*(-2*a + 3)*r1*t + (1/3456*(2*a - 3)*r1 + 1/3456*(2*a - 3)))*b2 + (1/2160*(-13*a + 15)*r1*t + (1/8640*(73*a - 105)*r1 + 1/1728*(-7*a + 15))))*b3 + (1/7200*(-a + 3)*r1*t + (1/28800*(3*a - 2)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/5760*(-7*a + 7)*r1 + 1/1152*(a - 1))*b2 + 1/360*(2*a - 3)*r1*t + 1/1440*(-13*a + 11)*r1 + 1/288*(3*a - 5)
];