d15fida:=[ ((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*(-27*a + 46)*r1 + 1/34560*(-5*a + 8))*t + (1/345600*(37*a - 66)*r1 + 1/69120*(19*a - 32)))*b2^2 + ((1/34560*(3*a - 5)*r1 + 1/34560*(71*a - 125))*t + (1/69120*(-23*a + 39)*r1 + 1/69120*(-91*a + 155)))*b2 + ((1/1728*(15*a - 26)*r1 + 1/8640*(-341*a + 590))*t + (1/17280*(5*a - 9)*r1 + 1/17280*(211*a - 365))))*b3 + (((1/172800*(79*a - 143)*r1 + 1/34560*(139*a - 239))*t + (1/345600*(-539*a + 923)*r1 + 1/23040*(-37*a + 65)))*b2^2 + ((1/17280*(-39*a + 68)*r1 + 1/17280*(-13*a + 22))*t + (1/34560*(44*a - 77)*r1 + 1/34560*(68*a - 117)))*b2 + ((1/1728*(-7*a + 13)*r1 + 1/8640*(-855*a + 1477))*t + (1/8640*(271*a - 466)*r1 + 1/8640*(320*a - 561)))))*b4 + ((((1/172800*a*r1 + 1/34560*(3*a - 2))*t + (1/345600*(-a + 10)*r1 + 1/69120*(-3*a + 4)))*b2^2 + ((1/34560*(a - 3)*r1 + 1/34560*(-3*a + 5))*t + (1/69120*(-a + 1)*r1 + 1/69120*(3*a + 5)))*b2 + ((1/8640*(-25*a + 42)*r1 + 1/8640*(153*a - 284))*t + (1/17280*(-17*a + 20)*r1 + 1/17280*(-123*a + 214))))*b3 + (((1/57600*(-17*a + 29)*r1 + 1/34560*(-79*a + 135))*t + (1/115200*(97*a - 169)*r1 + 1/69120*(71*a - 123)))*b2^2 + ((1/17280*(17*a - 30)*r1 + 1/5760*(5*a - 8))*t + (1/34560*(-20*a + 33)*r1 + 1/34560*(-40*a + 69)))*b2 + ((1/8640*(43*a - 63)*r1 + 1/8640*(499*a - 867))*t + (1/17280*(-319*a + 567)*r1 + 1/17280*(-439*a + 747))))))*b1 + ((((1/14400*(a - 2)*r1*t + (1/57600*(-4*a + 11)*r1 - 1/11520))*b2^2 + (1/2880*(-a + 1)*r1*t + (1/11520*(a - 1)*r1 + 1/2304*(a - 1)))*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(18*a - 35)*r1 + 1/576*(-14*a + 25))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/57600*(-a - 1)*r1 + 1/11520*(3*a - 5)))*b2^2 + (1/1440*(2*a - 3)*r1*t + (1/1920*(2*a - 3)*r1 + 1/1152*(-2*a + 3)))*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(-7*a + 11)*r1 + 1/576*(17*a - 29)))))*b4 + (((1/14400*(-a + 2)*r1*t + (1/57600*(4*a - 11)*r1 + 1/11520))*b2^2 + (1/2880*(a - 1)*r1*t + (1/11520*(-a + 1)*r1 + 1/2304*(-a + 1)))*b2 + (1/720*(-3*a + 4)*r1*t + (1/2880*(-12*a + 23)*r1 + 1/576*(8*a - 13))))*b3 + ((1/14400*(a - 1)*r1*t + (1/57600*(a + 1)*r1 + 1/11520*(-3*a + 5)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + (1/1920*(-2*a + 3)*r1 + 1/1152*(2*a - 3)))*b2 + (1/720*(9*a - 13)*r1*t + (1/576*(5*a - 7)*r1 + 1/576*(-11*a + 17)))))))*i + (((((1/172800*(4*a - 7)*r1 + 1/34560*(-22*a + 39))*t + (1/345600*(26*a - 43)*r1 + 1/69120*(24*a - 41)))*b2^2 + ((1/34560*(-17*a + 29)*r1 + 1/34560*(7*a - 15))*t + (1/69120*(17*a - 31)*r1 + 1/69120*(33*a - 55)))*b2 + ((1/8640*(86*a - 149)*r1 + 1/8640*(338*a - 585))*t + (1/17280*(-191*a + 331)*r1 + 1/17280*(-603*a + 1045))))*b3 + (((1/172800*(-79*a + 137)*r1 + 1/34560*(17*a - 27))*t + (1/345600*(-31*a + 53)*r1 + 1/69120*(21*a - 35)))*b2^2 + ((1/17280*(-24*a + 43)*r1 + 1/17280*(-250*a + 429))*t + (1/34560*(183*a - 314)*r1 + 1/34560*(215*a - 374)))*b2 + ((1/8640*(131*a - 227)*r1 + 1/8640*(167*a - 303))*t + (1/8640*(-71*a + 122)*r1 + 1/8640*(-226*a + 389)))))*b4 + ((((1/172800*(-2*a + 5)*r1 - 1/34560)*t + (1/345600*(2*a - 5)*r1 - 1/23040))*b2^2 + ((1/34560*(a - 1)*r1 + 1/34560*(9*a - 5))*t + (1/69120*(-a + 7)*r1 + 1/23040*(-3*a + 5)))*b2 + ((1/8640*(-40*a + 67)*r1 + 1/8640*(-114*a + 197))*t + (1/17280*(76*a - 133)*r1 + 1/1920*(26*a - 43))))*b3 + (((1/172800*(47*a - 85)*r1 + 1/34560*(-17*a + 31))*t + (1/345600*(23*a - 45)*r1 + 1/69120*(-5*a + 11)))*b2^2 + ((1/17280*(10*a - 17)*r1 + 1/17280*(116*a - 203))*t + (1/11520*(-29*a + 50)*r1 + 1/34560*(-91*a + 158)))*b2 + ((1/8640*(-71*a + 133)*r1 + 1/8640*(-3*a - 35))*t + (1/17280*(41*a - 51)*r1 + 1/17280*(153*a - 295))))))*b1 + (((1/43200*(4*a - 3)*r1*t + (1/172800*(-a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/1080*(-20*a + 33)*r1*t + (1/8640*(145*a - 246)*r1 + 1/1728*(5*a - 6))))*b3 + ((1/14400*(3*a - 5)*r1*t + (1/57600*(19*a - 27)*r1 + 1/11520*(-a + 1)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2)))*b2 + (1/180*(4*a - 7)*r1*t + (1/2880*(-83*a + 141)*r1 + 1/576*(-3*a + 5)))))*b4 + ((1/43200*(-4*a + 3)*r1*t + (1/172800*(a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/2160*(22*a - 39)*r1*t + (1/8640*(-73*a + 138)*r1 + 1/1728*(-5*a + 6))))*b3 + (1/14400*(-3*a + 5)*r1*t + (1/57600*(-19*a + 27)*r1 + 1/11520*(a - 1)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2)))*b2 + 1/720*(-7*a + 13)*r1*t + 1/2880*(65*a - 117)*r1 + 1/576*(-3*a + 7), ((((((1/57600*(-7*a + 11)*r1 + 1/34560*(-7*a + 13))*t + (1/86400*(9*a - 17)*r1 + 1/17280*(3*a - 5)))*b2^2 + (1/17280*(-31*a + 55)*t + (1/34560*(21*a - 37)*r1 + 1/34560*(21*a - 35)))*b2 + ((1/2880*(23*a - 40)*r1 + 1/8640*(497*a - 860))*t + (1/17280*(-381*a + 659)*r1 + 1/17280*(-447*a + 775))))*b3 + (((1/28800*(18*a - 31)*r1 + 1/17280*(-57*a + 98))*t + (1/172800*(21*a - 37)*r1 + 1/34560*(39*a - 67)))*b2^2 + ((1/1440*(-4*a + 7)*r1 + 1/8640*(-28*a + 47))*t + (1/17280*(31*a - 53)*r1 + 1/8640*(49*a - 86)))*b2 + ((1/2880*(-37*a + 63)*r1 + 1/8640*(425*a - 731))*t + (1/8640*(35*a - 62)*r1 + 1/8640*(-170*a + 293)))))*b4 + ((((1/172800*(3*a + 5)*r1 + 1/34560*(-3*a + 1))*t + (1/28800*a*r1 - 1/17280))*b2^2 + ((-1/8640*r1 + 1/17280*(3*a - 5))*t + (1/11520*(-a + 1)*r1 + 1/34560*(-3*a + 5)))*b2 + ((1/8640*(-27*a + 37)*r1 + 1/8640*(-201*a + 373))*t + (1/2880*(25*a - 47)*r1 + 1/8640*(93*a - 149))))*b3 + (((1/86400*(-26*a + 47)*r1 + 1/17280*(25*a - 44))*t + (1/172800*(a + 3)*r1 + 1/34560*(-17*a + 29)))*b2^2 + ((1/8640*(13*a - 23)*r1 + 1/8640*(25*a - 44))*t + (1/17280*(-23*a + 39)*r1 + 1/4320*(-15*a + 26)))*b2 + ((1/4320*(30*a - 53)*r1 + 1/4320*(-101*a + 178))*t + (1/960*(-3*a + 5)*r1 + 1/8640*(71*a - 133))))))*b1 + ((((1/14400*(a - 2)*r1*t + (1/57600*(5*a - 13)*r1 + 1/11520*(-a + 3)))*b2^2 + (1/2880*(-a + 1)*r1*t + 1/5760*(-a + 1)*r1)*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(-12*a + 25)*r1 + 1/576*(8*a - 15))))*b3 + ((1/14400*(a - 1)*r1*t + (-1/28800*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + 1/2880*(2*a - 3)*r1)*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(-a + 1)*r1 + 1/576*(11*a - 19)))))*b4 + (((1/14400*(-a + 2)*r1*t + (1/57600*(-5*a + 13)*r1 + 1/11520*(a - 3)))*b2^2 + (1/2880*(a - 1)*r1*t + 1/5760*(a - 1)*r1)*b2 + (1/720*(-3*a + 4)*r1*t + (1/2880*(9*a - 19)*r1 + 1/576*(-5*a + 9))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/28800*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(2*a - 3)*r1*t + 1/2880*(-2*a + 3)*r1)*b2 + (1/720*(-9*a + 13)*r1*t + (1/1440*(8*a - 11)*r1 + 1/288*(-a + 2)))))))*i + (((((1/172800*(7*a - 11)*r1 + 1/34560*(23*a - 41))*t + (1/172800*(-41*a + 73)*r1 + 1/11520*(-3*a + 5)))*b2^2 + ((1/17280*(-7*a + 11)*r1 + 1/8640*(-11*a + 20))*t + (1/34560*(19*a - 35)*r1 + 1/3840*(3*a - 5)))*b2 + ((1/8640*(56*a - 97)*r1 + 1/8640*(-136*a + 235))*t + (1/3456*(11*a - 19)*r1 + 1/5760*(-3*a + 5))))*b3 + (((1/86400*(-74*a + 127)*r1 + 1/17280*(-53*a + 92))*t + (1/172800*(169*a - 297)*r1 + 1/34560*(93*a - 161)))*b2^2 + ((1/8640*(-6*a + 11)*r1 + 1/8640*(105*a - 181))*t + (1/4320*(-5*a + 9)*r1 + 1/17280*(-105*a + 181)))*b2 + ((1/8640*(137*a - 233)*r1 + 1/8640*(469*a - 821))*t + (1/2880*(-56*a + 99)*r1 + 1/8640*(-417*a + 718)))))*b4 + ((((1/172800*(-a + 5)*r1 + 1/34560*(-a + 3))*t + (1/172800*(3*a - 5)*r1 + 1/34560*(a - 1)))*b2^2 + ((1/17280*(a + 1)*r1 - 1/4320*a)*t + (1/34560*(3*a + 1)*r1 + 1/34560*(-a - 5)))*b2 + ((1/8640*(-29*a + 43)*r1 + 1/8640*(23*a - 39))*t + (-1/2160*r1 + 1/2160*(4*a - 7))))*b3 + (((1/86400*(32*a - 55)*r1 + 1/17280*(23*a - 40))*t + (1/172800*(-67*a + 115)*r1 + 1/34560*(-43*a + 75)))*b2^2 + ((1/8640*(5*a - 8)*r1 + 1/4320*(-27*a + 46))*t + (1/8640*(5*a - 8)*r1 + 1/17280*(49*a - 87)))*b2 + ((1/864*(-6*a + 11)*r1 + 1/4320*(-113*a + 200))*t + (1/8640*(73*a - 115)*r1 + 1/2880*(71*a - 125))))))*b1 + (((1/43200*(4*a - 3)*r1*t + (1/172800*(5*a + 3)*r1 + 1/34560*(a - 3)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/4320*(a - 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/1080*(-20*a + 33)*r1*t + (1/8640*(-185*a + 312)*r1 + 1/1728*(35*a - 60))))*b3 + ((1/14400*(-3*a + 5)*r1*t + (1/28800*(8*a - 11)*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1440*(a - 2)*r1 + 1/576*(-a + 2)))*b2 + (1/180*(-4*a + 7)*r1*t + (1/2880*(-99*a + 169)*r1 + 1/576*(13*a - 23)))))*b4 + ((1/43200*(-4*a + 3)*r1*t + (1/172800*(-5*a - 3)*r1 + 1/34560*(-a + 3)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/4320*(-a + 3)*r1 + 1/3456*(a - 3)))*b2 + (1/2160*(22*a - 39)*r1*t + (1/8640*(95*a - 177)*r1 + 1/1728*(-17*a + 33))))*b3 + (1/14400*(3*a - 5)*r1*t + (1/28800*(-8*a + 11)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1440*(-a + 2)*r1 + 1/576*(a - 2)))*b2 + 1/720*(7*a - 13)*r1*t + 1/1440*(36*a - 65)*r1 + 1/288*(-5*a + 10), (((((1/14400*(-5*a + 8)*r1 + 1/8640*(-25*a + 44))*t + (1/86400*(35*a - 64)*r1 + 1/17280*(41*a - 70)))*b2^2 + ((1/1440*(-4*a + 7)*r1 + 1/4320*(28*a - 49))*t + (1/4320*(4*a - 7)*r1 + 1/4320*(-4*a + 7)))*b2 + ((1/1440*(26*a - 43)*r1 + 1/4320*(294*a - 521))*t + (1/8640*(-147*a + 266)*r1 + 1/8640*(-519*a + 886))))*b4 + (((1/43200*(9*a - 16)*r1 + 1/8640*(19*a - 32))*t + (1/28800*(-13*a + 22)*r1 + 1/17280*(-25*a + 44)))*b2^2 + ((1/4320*(5*a - 8)*r1 + 1/4320*(-17*a + 28))*t + (1/8640*(-a + 3)*r1 + 1/8640*(7*a - 13)))*b2 + ((1/4320*(-37*a + 71)*r1 + 1/4320*(-241*a + 395))*t + (1/8640*(125*a - 207)*r1 + 1/8640*(321*a - 575)))))*b1 + (((1/1800*(a - 2)*r1*t + (1/14400*(a - 2)*r1 + 1/2880*(-a + 2)))*b2^2 + (1/360*(-14*a + 25)*r1*t + (1/1440*(-97*a + 170)*r1 + 1/288*(25*a - 44))))*b4 + ((1/1800*(-a + 2)*r1*t + (1/14400*(-a + 2)*r1 + 1/2880*(a - 2)))*b2^2 + (1/360*(11*a - 19)*r1*t + (1/1440*(43*a - 77)*r1 + 1/288*(-13*a + 23))))))*i + ((((1/43200*(-32*a + 55)*r1 + 1/8640*(14*a - 25))*t + (1/86400*(17*a - 30)*r1 + 1/17280*(a - 2)))*b2^2 + ((1/4320*(9*a - 16)*r1 + 1/4320*(55*a - 94))*t + (1/2880*(-10*a + 17)*r1 + 1/8640*(-80*a + 139)))*b2 + ((1/4320*(47*a - 80)*r1 + 1/4320*(-393*a + 688))*t + (1/1728*(26*a - 45)*r1 + 1/8640*(218*a - 373))))*b4 + (((1/43200*(20*a - 37)*r1 + 1/8640*(-6*a + 11))*t + (1/86400*(-15*a + 22)*r1 + 1/5760*(-a + 2)))*b2^2 + ((1/4320*(-4*a + 7)*r1 + 1/864*(-4*a + 7))*t + (1/8640*(9*a - 16)*r1 + 1/1728*(7*a - 12)))*b2 + ((1/4320*(-39*a + 71)*r1 + 1/4320*(169*a - 281))*t + (1/8640*(-33*a + 61)*r1 + 1/8640*(-49*a + 81)))))*b1 + ((1/14400*(-8*a + 13)*r1 + 1/2880*(2*a - 3))*b2^2 + 1/720*(a - 2)*r1*b2 + (1/120*(-7*a + 12)*r1*t + (1/1440*(-20*a + 37)*r1 + 1/288*(-2*a + 3))))*b4 + (1/14400*(8*a - 13)*r1 + 1/2880*(-2*a + 3))*b2^2 + 1/720*(-a + 2)*r1*b2 + 1/120*(3*a - 5)*r1*t + 1/1440*(5*a - 13)*r1 + 1/288*(5*a - 9), ((((((1/172800*(27*a - 46)*r1 + 1/34560*(5*a - 8))*t + (1/345600*(-37*a + 66)*r1 + 1/69120*(-19*a + 32)))*b2^2 + ((1/34560*(-3*a + 5)*r1 + 1/34560*(-71*a + 125))*t + (1/69120*(23*a - 39)*r1 + 1/69120*(91*a - 155)))*b2 + ((1/1728*(-15*a + 26)*r1 + 1/8640*(341*a - 590))*t + (1/17280*(-5*a + 9)*r1 + 1/17280*(-211*a + 365))))*b3 + (((1/172800*(-79*a + 143)*r1 + 1/34560*(-139*a + 239))*t + (1/345600*(539*a - 923)*r1 + 1/23040*(37*a - 65)))*b2^2 + ((1/17280*(39*a - 68)*r1 + 1/17280*(13*a - 22))*t + (1/34560*(-44*a + 77)*r1 + 1/34560*(-68*a + 117)))*b2 + ((1/1728*(7*a - 13)*r1 + 1/8640*(855*a - 1477))*t + (1/8640*(-271*a + 466)*r1 + 1/8640*(-320*a + 561)))))*b4 + ((((-1/172800*a*r1 + 1/34560*(-3*a + 2))*t + (1/345600*(a - 10)*r1 + 1/69120*(3*a - 4)))*b2^2 + ((1/34560*(-a + 3)*r1 + 1/34560*(3*a - 5))*t + (1/69120*(a - 1)*r1 + 1/69120*(-3*a - 5)))*b2 + ((1/8640*(25*a - 42)*r1 + 1/8640*(-153*a + 284))*t + (1/17280*(17*a - 20)*r1 + 1/17280*(123*a - 214))))*b3 + (((1/57600*(17*a - 29)*r1 + 1/34560*(79*a - 135))*t + (1/115200*(-97*a + 169)*r1 + 1/69120*(-71*a + 123)))*b2^2 + ((1/17280*(-17*a + 30)*r1 + 1/5760*(-5*a + 8))*t + (1/34560*(20*a - 33)*r1 + 1/34560*(40*a - 69)))*b2 + ((1/8640*(-43*a + 63)*r1 + 1/8640*(-499*a + 867))*t + (1/17280*(319*a - 567)*r1 + 1/17280*(439*a - 747))))))*b1 + ((((1/14400*(a - 2)*r1*t + (1/57600*(-4*a + 11)*r1 - 1/11520))*b2^2 + (1/2880*(-a + 1)*r1*t + (1/11520*(a - 1)*r1 + 1/2304*(a - 1)))*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(18*a - 35)*r1 + 1/576*(-14*a + 25))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/57600*(-a - 1)*r1 + 1/11520*(3*a - 5)))*b2^2 + (1/1440*(2*a - 3)*r1*t + (1/1920*(2*a - 3)*r1 + 1/1152*(-2*a + 3)))*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(-7*a + 11)*r1 + 1/576*(17*a - 29)))))*b4 + (((1/14400*(-a + 2)*r1*t + (1/57600*(4*a - 11)*r1 + 1/11520))*b2^2 + (1/2880*(a - 1)*r1*t + (1/11520*(-a + 1)*r1 + 1/2304*(-a + 1)))*b2 + (1/720*(-3*a + 4)*r1*t + (1/2880*(-12*a + 23)*r1 + 1/576*(8*a - 13))))*b3 + ((1/14400*(a - 1)*r1*t + (1/57600*(a + 1)*r1 + 1/11520*(-3*a + 5)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + (1/1920*(-2*a + 3)*r1 + 1/1152*(2*a - 3)))*b2 + (1/720*(9*a - 13)*r1*t + (1/576*(5*a - 7)*r1 + 1/576*(-11*a + 17)))))))*i + (((((1/172800*(-4*a + 7)*r1 + 1/34560*(22*a - 39))*t + (1/345600*(-26*a + 43)*r1 + 1/69120*(-24*a + 41)))*b2^2 + ((1/34560*(17*a - 29)*r1 + 1/34560*(-7*a + 15))*t + (1/69120*(-17*a + 31)*r1 + 1/69120*(-33*a + 55)))*b2 + ((1/8640*(-86*a + 149)*r1 + 1/8640*(-338*a + 585))*t + (1/17280*(191*a - 331)*r1 + 1/17280*(603*a - 1045))))*b3 + (((1/172800*(79*a - 137)*r1 + 1/34560*(-17*a + 27))*t + (1/345600*(31*a - 53)*r1 + 1/69120*(-21*a + 35)))*b2^2 + ((1/17280*(24*a - 43)*r1 + 1/17280*(250*a - 429))*t + (1/34560*(-183*a + 314)*r1 + 1/34560*(-215*a + 374)))*b2 + ((1/8640*(-131*a + 227)*r1 + 1/8640*(-167*a + 303))*t + (1/8640*(71*a - 122)*r1 + 1/8640*(226*a - 389)))))*b4 + ((((1/172800*(2*a - 5)*r1 + 1/34560)*t + (1/345600*(-2*a + 5)*r1 + 1/23040))*b2^2 + ((1/34560*(-a + 1)*r1 + 1/34560*(-9*a + 5))*t + (1/69120*(a - 7)*r1 + 1/23040*(3*a - 5)))*b2 + ((1/8640*(40*a - 67)*r1 + 1/8640*(114*a - 197))*t + (1/17280*(-76*a + 133)*r1 + 1/1920*(-26*a + 43))))*b3 + (((1/172800*(-47*a + 85)*r1 + 1/34560*(17*a - 31))*t + (1/345600*(-23*a + 45)*r1 + 1/69120*(5*a - 11)))*b2^2 + ((1/17280*(-10*a + 17)*r1 + 1/17280*(-116*a + 203))*t + (1/11520*(29*a - 50)*r1 + 1/34560*(91*a - 158)))*b2 + ((1/8640*(71*a - 133)*r1 + 1/8640*(3*a + 35))*t + (1/17280*(-41*a + 51)*r1 + 1/17280*(-153*a + 295))))))*b1 + (((1/43200*(4*a - 3)*r1*t + (1/172800*(-a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/1080*(-20*a + 33)*r1*t + (1/8640*(145*a - 246)*r1 + 1/1728*(5*a - 6))))*b3 + ((1/14400*(3*a - 5)*r1*t + (1/57600*(19*a - 27)*r1 + 1/11520*(-a + 1)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2)))*b2 + (1/180*(4*a - 7)*r1*t + (1/2880*(-83*a + 141)*r1 + 1/576*(-3*a + 5)))))*b4 + ((1/43200*(-4*a + 3)*r1*t + (1/172800*(a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/2160*(22*a - 39)*r1*t + (1/8640*(-73*a + 138)*r1 + 1/1728*(-5*a + 6))))*b3 + (1/14400*(-3*a + 5)*r1*t + (1/57600*(-19*a + 27)*r1 + 1/11520*(a - 1)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2)))*b2 + 1/720*(-7*a + 13)*r1*t + 1/2880*(65*a - 117)*r1 + 1/576*(-3*a + 7), ((((1/28800*(a - 1)*r1 + 1/5760*(-a + 1))*b2^2 + (1/2880*(-a + 1)*r1 + 1/576*(a - 1))*b2 + (1/288*(4*a - 7)*r1 + 1/288*(-20*a + 35)))*b3 + ((1/4800*(a - 2)*r1 + 1/960*(-a + 2))*b2^2 + (1/480*(11*a - 19)*r1 + 1/96*(-11*a + 19))))*b4 + (((1/28800*(-a + 1)*r1 + 1/5760*(a - 1))*b2^2 + (1/2880*(a - 1)*r1 + 1/576*(-a + 1))*b2 + (1/1440*(-11*a + 17)*r1 + 1/288*(11*a - 17)))*b3 + ((1/4800*(-a + 2)*r1 + 1/960*(a - 2))*b2^2 + (1/240*(-a + 2)*r1 + 1/48*(a - 2)))))*i + (((1/9600*(a - 1)*r1 + 1/1920*(-a + 1))*b2^2 + (1/480*(-5*a + 8)*r1 + 1/96*(5*a - 8)))*b3 + ((1/14400*(5*a - 8)*r1 + 1/2880*(-5*a + 8))*b2^2 + (1/1440*(-a + 2)*r1 + 1/288*(a - 2))*b2 + (1/1440*(-33*a + 55)*r1 + 1/288*(33*a - 55))))*b4 + ((1/9600*(-a + 1)*r1 + 1/1920*(a - 1))*b2^2 + (1/480*(3*a - 5)*r1 + 1/96*(-3*a + 5)))*b3 + (1/14400*(-5*a + 8)*r1 + 1/2880*(5*a - 8))*b2^2 + (1/1440*(a - 2)*r1 + 1/288*(-a + 2))*b2 + 1/720*(9*a - 14)*r1 + 1/144*(-9*a + 14), (((((1/43200*(10*a - 17)*r1 + 1/8640*(-32*a + 55))*t + (1/86400*(85*a - 148)*r1 + 1/5760*(7*a - 12)))*b2^2 + ((1/2160*(3*a - 5)*r1 + 1/2160*(28*a - 49))*t + (1/864*(-4*a + 7)*r1 + 1/4320*(-23*a + 39)))*b2 + ((1/4320*(17*a - 29)*r1 + 1/4320*(303*a - 517))*t + (1/8640*(-191*a + 335)*r1 + 1/8640*(-313*a + 537))))*b4 + (((1/43200*(-2*a + 3)*r1 + 1/2880*(4*a - 7))*t + (1/86400*(-43*a + 72)*r1 + 1/17280*(-7*a + 12)))*b2^2 + ((1/4320*(-5*a + 9)*r1 + 1/4320*(-29*a + 51))*t + (1/8640*(23*a - 39)*r1 + 1/8640*(29*a - 51)))*b2 + ((1/360*(-a + 2)*r1 + 1/2160*(-61*a + 105))*t + (1/1440*(18*a - 29)*r1 + 1/4320*(56*a - 105)))))*b1 + (((1/1800*(-a + 2)*r1*t + (1/14400*(-a + 2)*r1 + 1/2880*(a - 2)))*b2^2 + (1/360*(14*a - 25)*r1*t + (1/1440*(-83*a + 145)*r1 + 1/288*(11*a - 19))))*b4 + ((1/1800*(a - 2)*r1*t + (1/14400*(a - 2)*r1 + 1/2880*(-a + 2)))*b2^2 + (1/360*(-11*a + 19)*r1*t + (1/720*(16*a - 29)*r1 + 1/144*(-a + 2))))))*i + ((((1/43200*(41*a - 70)*r1 + 1/8640*(31*a - 54))*t + (1/86400*(-131*a + 230)*r1 + 1/17280*(-37*a + 64)))*b2^2 + ((1/2160*(-3*a + 5)*r1 + 1/2160*(5*a - 9))*t + (1/8640*(-8*a + 13)*r1 + 1/8640*(10*a - 17)))*b2 + ((1/4320*(-83*a + 143)*r1 + 1/4320*(-151*a + 261))*t + (1/8640*(237*a - 415)*r1 + 1/8640*(171*a - 301))))*b4 + (((1/43200*(-15*a + 26)*r1 + 1/8640*(-13*a + 22))*t + (1/86400*(55*a - 96)*r1 + 1/17280*(13*a - 22)))*b2^2 + ((1/4320*(3*a - 5)*r1 + 1/864*(-3*a + 5))*t + (1/8640*(7*a - 12)*r1 + 1/1728*(a - 2)))*b2 + ((1/2160*(15*a - 29)*r1 + 1/1080*(17*a - 28))*t + (1/1440*(-18*a + 29)*r1 + 1/4320*(-24*a + 41)))))*b1 + ((1/14400*(-8*a + 13)*r1 + 1/2880*(2*a - 3))*b2^2 + 1/720*(a - 2)*r1*b2 + (1/120*(7*a - 12)*r1*t + (1/1440*(a + 1)*r1 + 1/288*(-23*a + 39))))*b4 + (1/14400*(8*a - 13)*r1 + 1/2880*(-2*a + 3))*b2^2 + 1/720*(-a + 2)*r1*b2 + 1/120*(-3*a + 5)*r1*t + 1/720*(-2*a + 1)*r1 + 1/144*(7*a - 12), ((((((1/57600*(7*a - 11)*r1 + 1/34560*(7*a - 13))*t + (1/86400*(-9*a + 17)*r1 + 1/17280*(-3*a + 5)))*b2^2 + (1/17280*(31*a - 55)*t + (1/34560*(-21*a + 37)*r1 + 1/34560*(-21*a + 35)))*b2 + ((1/2880*(-23*a + 40)*r1 + 1/8640*(-497*a + 860))*t + (1/17280*(381*a - 659)*r1 + 1/17280*(447*a - 775))))*b3 + (((1/28800*(18*a - 31)*r1 + 1/17280*(-57*a + 98))*t + (1/172800*(21*a - 37)*r1 + 1/34560*(39*a - 67)))*b2^2 + ((1/1440*(-4*a + 7)*r1 + 1/8640*(-28*a + 47))*t + (1/17280*(31*a - 53)*r1 + 1/8640*(49*a - 86)))*b2 + ((1/2880*(-37*a + 63)*r1 + 1/8640*(425*a - 731))*t + (1/8640*(35*a - 62)*r1 + 1/8640*(-170*a + 293)))))*b4 + ((((1/172800*(-3*a - 5)*r1 + 1/34560*(3*a - 1))*t + (-1/28800*a*r1 + 1/17280))*b2^2 + ((1/8640*r1 + 1/17280*(-3*a + 5))*t + (1/11520*(a - 1)*r1 + 1/34560*(3*a - 5)))*b2 + ((1/8640*(27*a - 37)*r1 + 1/8640*(201*a - 373))*t + (1/2880*(-25*a + 47)*r1 + 1/8640*(-93*a + 149))))*b3 + (((1/86400*(-26*a + 47)*r1 + 1/17280*(25*a - 44))*t + (1/172800*(a + 3)*r1 + 1/34560*(-17*a + 29)))*b2^2 + ((1/8640*(13*a - 23)*r1 + 1/8640*(25*a - 44))*t + (1/17280*(-23*a + 39)*r1 + 1/4320*(-15*a + 26)))*b2 + ((1/4320*(30*a - 53)*r1 + 1/4320*(-101*a + 178))*t + (1/960*(-3*a + 5)*r1 + 1/8640*(71*a - 133))))))*b1 + ((((1/14400*(-a + 2)*r1*t + (1/57600*(-5*a + 13)*r1 + 1/11520*(a - 3)))*b2^2 + (1/2880*(a - 1)*r1*t + 1/5760*(a - 1)*r1)*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(12*a - 25)*r1 + 1/576*(-8*a + 15))))*b3 + ((1/14400*(a - 1)*r1*t + (-1/28800*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + 1/2880*(2*a - 3)*r1)*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(-a + 1)*r1 + 1/576*(11*a - 19)))))*b4 + (((1/14400*(a - 2)*r1*t + (1/57600*(5*a - 13)*r1 + 1/11520*(-a + 3)))*b2^2 + (1/2880*(-a + 1)*r1*t + 1/5760*(-a + 1)*r1)*b2 + (1/720*(3*a - 4)*r1*t + (1/2880*(-9*a + 19)*r1 + 1/576*(5*a - 9))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/28800*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(2*a - 3)*r1*t + 1/2880*(-2*a + 3)*r1)*b2 + (1/720*(-9*a + 13)*r1*t + (1/1440*(8*a - 11)*r1 + 1/288*(-a + 2)))))))*i + (((((1/172800*(-7*a + 11)*r1 + 1/34560*(-23*a + 41))*t + (1/172800*(41*a - 73)*r1 + 1/11520*(3*a - 5)))*b2^2 + ((1/17280*(7*a - 11)*r1 + 1/8640*(11*a - 20))*t + (1/34560*(-19*a + 35)*r1 + 1/3840*(-3*a + 5)))*b2 + ((1/8640*(-56*a + 97)*r1 + 1/8640*(136*a - 235))*t + (1/3456*(-11*a + 19)*r1 + 1/5760*(3*a - 5))))*b3 + (((1/86400*(-74*a + 127)*r1 + 1/17280*(-53*a + 92))*t + (1/172800*(169*a - 297)*r1 + 1/34560*(93*a - 161)))*b2^2 + ((1/8640*(-6*a + 11)*r1 + 1/8640*(105*a - 181))*t + (1/4320*(-5*a + 9)*r1 + 1/17280*(-105*a + 181)))*b2 + ((1/8640*(137*a - 233)*r1 + 1/8640*(469*a - 821))*t + (1/2880*(-56*a + 99)*r1 + 1/8640*(-417*a + 718)))))*b4 + ((((1/172800*(a - 5)*r1 + 1/34560*(a - 3))*t + (1/172800*(-3*a + 5)*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/17280*(-a - 1)*r1 + 1/4320*a)*t + (1/34560*(-3*a - 1)*r1 + 1/34560*(a + 5)))*b2 + ((1/8640*(29*a - 43)*r1 + 1/8640*(-23*a + 39))*t + (1/2160*r1 + 1/2160*(-4*a + 7))))*b3 + (((1/86400*(32*a - 55)*r1 + 1/17280*(23*a - 40))*t + (1/172800*(-67*a + 115)*r1 + 1/34560*(-43*a + 75)))*b2^2 + ((1/8640*(5*a - 8)*r1 + 1/4320*(-27*a + 46))*t + (1/8640*(5*a - 8)*r1 + 1/17280*(49*a - 87)))*b2 + ((1/864*(-6*a + 11)*r1 + 1/4320*(-113*a + 200))*t + (1/8640*(73*a - 115)*r1 + 1/2880*(71*a - 125))))))*b1 + (((1/43200*(-4*a + 3)*r1*t + (1/172800*(-5*a - 3)*r1 + 1/34560*(-a + 3)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/4320*(-a + 3)*r1 + 1/3456*(a - 3)))*b2 + (1/1080*(20*a - 33)*r1*t + (1/8640*(185*a - 312)*r1 + 1/1728*(-35*a + 60))))*b3 + ((1/14400*(-3*a + 5)*r1*t + (1/28800*(8*a - 11)*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1440*(a - 2)*r1 + 1/576*(-a + 2)))*b2 + (1/180*(-4*a + 7)*r1*t + (1/2880*(-99*a + 169)*r1 + 1/576*(13*a - 23)))))*b4 + ((1/43200*(4*a - 3)*r1*t + (1/172800*(5*a + 3)*r1 + 1/34560*(a - 3)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/4320*(a - 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/2160*(-22*a + 39)*r1*t + (1/8640*(-95*a + 177)*r1 + 1/1728*(17*a - 33))))*b3 + (1/14400*(3*a - 5)*r1*t + (1/28800*(-8*a + 11)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1440*(-a + 2)*r1 + 1/576*(a - 2)))*b2 + 1/720*(7*a - 13)*r1*t + 1/1440*(36*a - 65)*r1 + 1/288*(-5*a + 10), ((((((1/57600*(7*a - 11)*r1 + 1/34560*(7*a - 13))*t + (1/86400*(-9*a + 17)*r1 + 1/17280*(-3*a + 5)))*b2^2 + (1/17280*(31*a - 55)*t + (1/34560*(-21*a + 37)*r1 + 1/34560*(-21*a + 35)))*b2 + ((1/2880*(-23*a + 40)*r1 + 1/8640*(-497*a + 860))*t + (1/17280*(381*a - 659)*r1 + 1/17280*(447*a - 775))))*b3 + (((1/28800*(-18*a + 31)*r1 + 1/17280*(57*a - 98))*t + (1/172800*(-21*a + 37)*r1 + 1/34560*(-39*a + 67)))*b2^2 + ((1/1440*(4*a - 7)*r1 + 1/8640*(28*a - 47))*t + (1/17280*(-31*a + 53)*r1 + 1/8640*(-49*a + 86)))*b2 + ((1/2880*(37*a - 63)*r1 + 1/8640*(-425*a + 731))*t + (1/8640*(-35*a + 62)*r1 + 1/8640*(170*a - 293)))))*b4 + ((((1/172800*(-3*a - 5)*r1 + 1/34560*(3*a - 1))*t + (-1/28800*a*r1 + 1/17280))*b2^2 + ((1/8640*r1 + 1/17280*(-3*a + 5))*t + (1/11520*(a - 1)*r1 + 1/34560*(3*a - 5)))*b2 + ((1/8640*(27*a - 37)*r1 + 1/8640*(201*a - 373))*t + (1/2880*(-25*a + 47)*r1 + 1/8640*(-93*a + 149))))*b3 + (((1/86400*(26*a - 47)*r1 + 1/17280*(-25*a + 44))*t + (1/172800*(-a - 3)*r1 + 1/34560*(17*a - 29)))*b2^2 + ((1/8640*(-13*a + 23)*r1 + 1/8640*(-25*a + 44))*t + (1/17280*(23*a - 39)*r1 + 1/4320*(15*a - 26)))*b2 + ((1/4320*(-30*a + 53)*r1 + 1/4320*(101*a - 178))*t + (1/960*(3*a - 5)*r1 + 1/8640*(-71*a + 133))))))*b1 + ((((1/14400*(a - 2)*r1*t + (1/57600*(5*a - 13)*r1 + 1/11520*(-a + 3)))*b2^2 + (1/2880*(-a + 1)*r1*t + 1/5760*(-a + 1)*r1)*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(-12*a + 25)*r1 + 1/576*(8*a - 15))))*b3 + ((1/14400*(a - 1)*r1*t + (-1/28800*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + 1/2880*(2*a - 3)*r1)*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(-a + 1)*r1 + 1/576*(11*a - 19)))))*b4 + (((1/14400*(-a + 2)*r1*t + (1/57600*(-5*a + 13)*r1 + 1/11520*(a - 3)))*b2^2 + (1/2880*(a - 1)*r1*t + 1/5760*(a - 1)*r1)*b2 + (1/720*(-3*a + 4)*r1*t + (1/2880*(9*a - 19)*r1 + 1/576*(-5*a + 9))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/28800*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(2*a - 3)*r1*t + 1/2880*(-2*a + 3)*r1)*b2 + (1/720*(-9*a + 13)*r1*t + (1/1440*(8*a - 11)*r1 + 1/288*(-a + 2)))))))*i + (((((1/172800*(-7*a + 11)*r1 + 1/34560*(-23*a + 41))*t + (1/172800*(41*a - 73)*r1 + 1/11520*(3*a - 5)))*b2^2 + ((1/17280*(7*a - 11)*r1 + 1/8640*(11*a - 20))*t + (1/34560*(-19*a + 35)*r1 + 1/3840*(-3*a + 5)))*b2 + ((1/8640*(-56*a + 97)*r1 + 1/8640*(136*a - 235))*t + (1/3456*(-11*a + 19)*r1 + 1/5760*(3*a - 5))))*b3 + (((1/86400*(74*a - 127)*r1 + 1/17280*(53*a - 92))*t + (1/172800*(-169*a + 297)*r1 + 1/34560*(-93*a + 161)))*b2^2 + ((1/8640*(6*a - 11)*r1 + 1/8640*(-105*a + 181))*t + (1/4320*(5*a - 9)*r1 + 1/17280*(105*a - 181)))*b2 + ((1/8640*(-137*a + 233)*r1 + 1/8640*(-469*a + 821))*t + (1/2880*(56*a - 99)*r1 + 1/8640*(417*a - 718)))))*b4 + ((((1/172800*(a - 5)*r1 + 1/34560*(a - 3))*t + (1/172800*(-3*a + 5)*r1 + 1/34560*(-a + 1)))*b2^2 + ((1/17280*(-a - 1)*r1 + 1/4320*a)*t + (1/34560*(-3*a - 1)*r1 + 1/34560*(a + 5)))*b2 + ((1/8640*(29*a - 43)*r1 + 1/8640*(-23*a + 39))*t + (1/2160*r1 + 1/2160*(-4*a + 7))))*b3 + (((1/86400*(-32*a + 55)*r1 + 1/17280*(-23*a + 40))*t + (1/172800*(67*a - 115)*r1 + 1/34560*(43*a - 75)))*b2^2 + ((1/8640*(-5*a + 8)*r1 + 1/4320*(27*a - 46))*t + (1/8640*(-5*a + 8)*r1 + 1/17280*(-49*a + 87)))*b2 + ((1/864*(6*a - 11)*r1 + 1/4320*(113*a - 200))*t + (1/8640*(-73*a + 115)*r1 + 1/2880*(-71*a + 125))))))*b1 + (((1/43200*(4*a - 3)*r1*t + (1/172800*(5*a + 3)*r1 + 1/34560*(a - 3)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/4320*(a - 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/1080*(-20*a + 33)*r1*t + (1/8640*(-185*a + 312)*r1 + 1/1728*(35*a - 60))))*b3 + ((1/14400*(-3*a + 5)*r1*t + (1/28800*(8*a - 11)*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1440*(a - 2)*r1 + 1/576*(-a + 2)))*b2 + (1/180*(-4*a + 7)*r1*t + (1/2880*(-99*a + 169)*r1 + 1/576*(13*a - 23)))))*b4 + ((1/43200*(-4*a + 3)*r1*t + (1/172800*(-5*a - 3)*r1 + 1/34560*(-a + 3)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/4320*(-a + 3)*r1 + 1/3456*(a - 3)))*b2 + (1/2160*(22*a - 39)*r1*t + (1/8640*(95*a - 177)*r1 + 1/1728*(-17*a + 33))))*b3 + (1/14400*(3*a - 5)*r1*t + (1/28800*(-8*a + 11)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1440*(-a + 2)*r1 + 1/576*(a - 2)))*b2 + 1/720*(7*a - 13)*r1*t + 1/1440*(36*a - 65)*r1 + 1/288*(-5*a + 10), (((((1/43200*(-10*a + 17)*r1 + 1/8640*(32*a - 55))*t + (1/86400*(-85*a + 148)*r1 + 1/5760*(-7*a + 12)))*b2^2 + ((1/2160*(-3*a + 5)*r1 + 1/2160*(-28*a + 49))*t + (1/864*(4*a - 7)*r1 + 1/4320*(23*a - 39)))*b2 + ((1/4320*(-17*a + 29)*r1 + 1/4320*(-303*a + 517))*t + (1/8640*(191*a - 335)*r1 + 1/8640*(313*a - 537))))*b4 + (((1/43200*(2*a - 3)*r1 + 1/2880*(-4*a + 7))*t + (1/86400*(43*a - 72)*r1 + 1/17280*(7*a - 12)))*b2^2 + ((1/4320*(5*a - 9)*r1 + 1/4320*(29*a - 51))*t + (1/8640*(-23*a + 39)*r1 + 1/8640*(-29*a + 51)))*b2 + ((1/360*(a - 2)*r1 + 1/2160*(61*a - 105))*t + (1/1440*(-18*a + 29)*r1 + 1/4320*(-56*a + 105)))))*b1 + (((1/1800*(-a + 2)*r1*t + (1/14400*(-a + 2)*r1 + 1/2880*(a - 2)))*b2^2 + (1/360*(14*a - 25)*r1*t + (1/1440*(-83*a + 145)*r1 + 1/288*(11*a - 19))))*b4 + ((1/1800*(a - 2)*r1*t + (1/14400*(a - 2)*r1 + 1/2880*(-a + 2)))*b2^2 + (1/360*(-11*a + 19)*r1*t + (1/720*(16*a - 29)*r1 + 1/144*(-a + 2))))))*i + ((((1/43200*(-41*a + 70)*r1 + 1/8640*(-31*a + 54))*t + (1/86400*(131*a - 230)*r1 + 1/17280*(37*a - 64)))*b2^2 + ((1/2160*(3*a - 5)*r1 + 1/2160*(-5*a + 9))*t + (1/8640*(8*a - 13)*r1 + 1/8640*(-10*a + 17)))*b2 + ((1/4320*(83*a - 143)*r1 + 1/4320*(151*a - 261))*t + (1/8640*(-237*a + 415)*r1 + 1/8640*(-171*a + 301))))*b4 + (((1/43200*(15*a - 26)*r1 + 1/8640*(13*a - 22))*t + (1/86400*(-55*a + 96)*r1 + 1/17280*(-13*a + 22)))*b2^2 + ((1/4320*(-3*a + 5)*r1 + 1/864*(3*a - 5))*t + (1/8640*(-7*a + 12)*r1 + 1/1728*(-a + 2)))*b2 + ((1/2160*(-15*a + 29)*r1 + 1/1080*(-17*a + 28))*t + (1/1440*(18*a - 29)*r1 + 1/4320*(24*a - 41)))))*b1 + ((1/14400*(-8*a + 13)*r1 + 1/2880*(2*a - 3))*b2^2 + 1/720*(a - 2)*r1*b2 + (1/120*(7*a - 12)*r1*t + (1/1440*(a + 1)*r1 + 1/288*(-23*a + 39))))*b4 + (1/14400*(8*a - 13)*r1 + 1/2880*(-2*a + 3))*b2^2 + 1/720*(-a + 2)*r1*b2 + 1/120*(-3*a + 5)*r1*t + 1/720*(-2*a + 1)*r1 + 1/144*(7*a - 12), ((((1/28800*(-a + 1)*r1 + 1/5760*(a - 1))*b2^2 + (1/2880*(a - 1)*r1 + 1/576*(-a + 1))*b2 + (1/288*(-4*a + 7)*r1 + 1/288*(20*a - 35)))*b3 + ((1/4800*(a - 2)*r1 + 1/960*(-a + 2))*b2^2 + (1/480*(11*a - 19)*r1 + 1/96*(-11*a + 19))))*b4 + (((1/28800*(a - 1)*r1 + 1/5760*(-a + 1))*b2^2 + (1/2880*(-a + 1)*r1 + 1/576*(a - 1))*b2 + (1/1440*(11*a - 17)*r1 + 1/288*(-11*a + 17)))*b3 + ((1/4800*(-a + 2)*r1 + 1/960*(a - 2))*b2^2 + (1/240*(-a + 2)*r1 + 1/48*(a - 2)))))*i + (((1/9600*(-a + 1)*r1 + 1/1920*(a - 1))*b2^2 + (1/480*(5*a - 8)*r1 + 1/96*(-5*a + 8)))*b3 + ((1/14400*(5*a - 8)*r1 + 1/2880*(-5*a + 8))*b2^2 + (1/1440*(-a + 2)*r1 + 1/288*(a - 2))*b2 + (1/1440*(-33*a + 55)*r1 + 1/288*(33*a - 55))))*b4 + ((1/9600*(a - 1)*r1 + 1/1920*(-a + 1))*b2^2 + (1/480*(-3*a + 5)*r1 + 1/96*(3*a - 5)))*b3 + (1/14400*(-5*a + 8)*r1 + 1/2880*(5*a - 8))*b2^2 + (1/1440*(a - 2)*r1 + 1/288*(-a + 2))*b2 + 1/720*(9*a - 14)*r1 + 1/144*(-9*a + 14), ((((((1/172800*(27*a - 46)*r1 + 1/34560*(5*a - 8))*t + (1/345600*(-37*a + 66)*r1 + 1/69120*(-19*a + 32)))*b2^2 + ((1/34560*(-3*a + 5)*r1 + 1/34560*(-71*a + 125))*t + (1/69120*(23*a - 39)*r1 + 1/69120*(91*a - 155)))*b2 + ((1/1728*(-15*a + 26)*r1 + 1/8640*(341*a - 590))*t + (1/17280*(-5*a + 9)*r1 + 1/17280*(-211*a + 365))))*b3 + (((1/172800*(79*a - 143)*r1 + 1/34560*(139*a - 239))*t + (1/345600*(-539*a + 923)*r1 + 1/23040*(-37*a + 65)))*b2^2 + ((1/17280*(-39*a + 68)*r1 + 1/17280*(-13*a + 22))*t + (1/34560*(44*a - 77)*r1 + 1/34560*(68*a - 117)))*b2 + ((1/1728*(-7*a + 13)*r1 + 1/8640*(-855*a + 1477))*t + (1/8640*(271*a - 466)*r1 + 1/8640*(320*a - 561)))))*b4 + ((((-1/172800*a*r1 + 1/34560*(-3*a + 2))*t + (1/345600*(a - 10)*r1 + 1/69120*(3*a - 4)))*b2^2 + ((1/34560*(-a + 3)*r1 + 1/34560*(3*a - 5))*t + (1/69120*(a - 1)*r1 + 1/69120*(-3*a - 5)))*b2 + ((1/8640*(25*a - 42)*r1 + 1/8640*(-153*a + 284))*t + (1/17280*(17*a - 20)*r1 + 1/17280*(123*a - 214))))*b3 + (((1/57600*(-17*a + 29)*r1 + 1/34560*(-79*a + 135))*t + (1/115200*(97*a - 169)*r1 + 1/69120*(71*a - 123)))*b2^2 + ((1/17280*(17*a - 30)*r1 + 1/5760*(5*a - 8))*t + (1/34560*(-20*a + 33)*r1 + 1/34560*(-40*a + 69)))*b2 + ((1/8640*(43*a - 63)*r1 + 1/8640*(499*a - 867))*t + (1/17280*(-319*a + 567)*r1 + 1/17280*(-439*a + 747))))))*b1 + ((((1/14400*(-a + 2)*r1*t + (1/57600*(4*a - 11)*r1 + 1/11520))*b2^2 + (1/2880*(a - 1)*r1*t + (1/11520*(-a + 1)*r1 + 1/2304*(-a + 1)))*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(-18*a + 35)*r1 + 1/576*(14*a - 25))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/57600*(-a - 1)*r1 + 1/11520*(3*a - 5)))*b2^2 + (1/1440*(2*a - 3)*r1*t + (1/1920*(2*a - 3)*r1 + 1/1152*(-2*a + 3)))*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(-7*a + 11)*r1 + 1/576*(17*a - 29)))))*b4 + (((1/14400*(a - 2)*r1*t + (1/57600*(-4*a + 11)*r1 - 1/11520))*b2^2 + (1/2880*(-a + 1)*r1*t + (1/11520*(a - 1)*r1 + 1/2304*(a - 1)))*b2 + (1/720*(3*a - 4)*r1*t + (1/2880*(12*a - 23)*r1 + 1/576*(-8*a + 13))))*b3 + ((1/14400*(a - 1)*r1*t + (1/57600*(a + 1)*r1 + 1/11520*(-3*a + 5)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + (1/1920*(-2*a + 3)*r1 + 1/1152*(2*a - 3)))*b2 + (1/720*(9*a - 13)*r1*t + (1/576*(5*a - 7)*r1 + 1/576*(-11*a + 17)))))))*i + (((((1/172800*(-4*a + 7)*r1 + 1/34560*(22*a - 39))*t + (1/345600*(-26*a + 43)*r1 + 1/69120*(-24*a + 41)))*b2^2 + ((1/34560*(17*a - 29)*r1 + 1/34560*(-7*a + 15))*t + (1/69120*(-17*a + 31)*r1 + 1/69120*(-33*a + 55)))*b2 + ((1/8640*(-86*a + 149)*r1 + 1/8640*(-338*a + 585))*t + (1/17280*(191*a - 331)*r1 + 1/17280*(603*a - 1045))))*b3 + (((1/172800*(-79*a + 137)*r1 + 1/34560*(17*a - 27))*t + (1/345600*(-31*a + 53)*r1 + 1/69120*(21*a - 35)))*b2^2 + ((1/17280*(-24*a + 43)*r1 + 1/17280*(-250*a + 429))*t + (1/34560*(183*a - 314)*r1 + 1/34560*(215*a - 374)))*b2 + ((1/8640*(131*a - 227)*r1 + 1/8640*(167*a - 303))*t + (1/8640*(-71*a + 122)*r1 + 1/8640*(-226*a + 389)))))*b4 + ((((1/172800*(2*a - 5)*r1 + 1/34560)*t + (1/345600*(-2*a + 5)*r1 + 1/23040))*b2^2 + ((1/34560*(-a + 1)*r1 + 1/34560*(-9*a + 5))*t + (1/69120*(a - 7)*r1 + 1/23040*(3*a - 5)))*b2 + ((1/8640*(40*a - 67)*r1 + 1/8640*(114*a - 197))*t + (1/17280*(-76*a + 133)*r1 + 1/1920*(-26*a + 43))))*b3 + (((1/172800*(47*a - 85)*r1 + 1/34560*(-17*a + 31))*t + (1/345600*(23*a - 45)*r1 + 1/69120*(-5*a + 11)))*b2^2 + ((1/17280*(10*a - 17)*r1 + 1/17280*(116*a - 203))*t + (1/11520*(-29*a + 50)*r1 + 1/34560*(-91*a + 158)))*b2 + ((1/8640*(-71*a + 133)*r1 + 1/8640*(-3*a - 35))*t + (1/17280*(41*a - 51)*r1 + 1/17280*(153*a - 295))))))*b1 + (((1/43200*(-4*a + 3)*r1*t + (1/172800*(a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/1080*(20*a - 33)*r1*t + (1/8640*(-145*a + 246)*r1 + 1/1728*(-5*a + 6))))*b3 + ((1/14400*(3*a - 5)*r1*t + (1/57600*(19*a - 27)*r1 + 1/11520*(-a + 1)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2)))*b2 + (1/180*(4*a - 7)*r1*t + (1/2880*(-83*a + 141)*r1 + 1/576*(-3*a + 5)))))*b4 + ((1/43200*(4*a - 3)*r1*t + (1/172800*(-a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/2160*(-22*a + 39)*r1*t + (1/8640*(73*a - 138)*r1 + 1/1728*(5*a - 6))))*b3 + (1/14400*(-3*a + 5)*r1*t + (1/57600*(-19*a + 27)*r1 + 1/11520*(a - 1)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2)))*b2 + 1/720*(-7*a + 13)*r1*t + 1/2880*(65*a - 117)*r1 + 1/576*(-3*a + 7), (((((1/14400*(5*a - 8)*r1 + 1/8640*(25*a - 44))*t + (1/86400*(-35*a + 64)*r1 + 1/17280*(-41*a + 70)))*b2^2 + ((1/1440*(4*a - 7)*r1 + 1/4320*(-28*a + 49))*t + (1/4320*(-4*a + 7)*r1 + 1/4320*(4*a - 7)))*b2 + ((1/1440*(-26*a + 43)*r1 + 1/4320*(-294*a + 521))*t + (1/8640*(147*a - 266)*r1 + 1/8640*(519*a - 886))))*b4 + (((1/43200*(-9*a + 16)*r1 + 1/8640*(-19*a + 32))*t + (1/28800*(13*a - 22)*r1 + 1/17280*(25*a - 44)))*b2^2 + ((1/4320*(-5*a + 8)*r1 + 1/4320*(17*a - 28))*t + (1/8640*(a - 3)*r1 + 1/8640*(-7*a + 13)))*b2 + ((1/4320*(37*a - 71)*r1 + 1/4320*(241*a - 395))*t + (1/8640*(-125*a + 207)*r1 + 1/8640*(-321*a + 575)))))*b1 + (((1/1800*(a - 2)*r1*t + (1/14400*(a - 2)*r1 + 1/2880*(-a + 2)))*b2^2 + (1/360*(-14*a + 25)*r1*t + (1/1440*(-97*a + 170)*r1 + 1/288*(25*a - 44))))*b4 + ((1/1800*(-a + 2)*r1*t + (1/14400*(-a + 2)*r1 + 1/2880*(a - 2)))*b2^2 + (1/360*(11*a - 19)*r1*t + (1/1440*(43*a - 77)*r1 + 1/288*(-13*a + 23))))))*i + ((((1/43200*(32*a - 55)*r1 + 1/8640*(-14*a + 25))*t + (1/86400*(-17*a + 30)*r1 + 1/17280*(-a + 2)))*b2^2 + ((1/4320*(-9*a + 16)*r1 + 1/4320*(-55*a + 94))*t + (1/2880*(10*a - 17)*r1 + 1/8640*(80*a - 139)))*b2 + ((1/4320*(-47*a + 80)*r1 + 1/4320*(393*a - 688))*t + (1/1728*(-26*a + 45)*r1 + 1/8640*(-218*a + 373))))*b4 + (((1/43200*(-20*a + 37)*r1 + 1/8640*(6*a - 11))*t + (1/86400*(15*a - 22)*r1 + 1/5760*(a - 2)))*b2^2 + ((1/4320*(4*a - 7)*r1 + 1/864*(4*a - 7))*t + (1/8640*(-9*a + 16)*r1 + 1/1728*(-7*a + 12)))*b2 + ((1/4320*(39*a - 71)*r1 + 1/4320*(-169*a + 281))*t + (1/8640*(33*a - 61)*r1 + 1/8640*(49*a - 81)))))*b1 + ((1/14400*(-8*a + 13)*r1 + 1/2880*(2*a - 3))*b2^2 + 1/720*(a - 2)*r1*b2 + (1/120*(-7*a + 12)*r1*t + (1/1440*(-20*a + 37)*r1 + 1/288*(-2*a + 3))))*b4 + (1/14400*(8*a - 13)*r1 + 1/2880*(-2*a + 3))*b2^2 + 1/720*(-a + 2)*r1*b2 + 1/120*(3*a - 5)*r1*t + 1/1440*(5*a - 13)*r1 + 1/288*(5*a - 9), ((((((1/57600*(-7*a + 11)*r1 + 1/34560*(-7*a + 13))*t + (1/86400*(9*a - 17)*r1 + 1/17280*(3*a - 5)))*b2^2 + (1/17280*(-31*a + 55)*t + (1/34560*(21*a - 37)*r1 + 1/34560*(21*a - 35)))*b2 + ((1/2880*(23*a - 40)*r1 + 1/8640*(497*a - 860))*t + (1/17280*(-381*a + 659)*r1 + 1/17280*(-447*a + 775))))*b3 + (((1/28800*(-18*a + 31)*r1 + 1/17280*(57*a - 98))*t + (1/172800*(-21*a + 37)*r1 + 1/34560*(-39*a + 67)))*b2^2 + ((1/1440*(4*a - 7)*r1 + 1/8640*(28*a - 47))*t + (1/17280*(-31*a + 53)*r1 + 1/8640*(-49*a + 86)))*b2 + ((1/2880*(37*a - 63)*r1 + 1/8640*(-425*a + 731))*t + (1/8640*(-35*a + 62)*r1 + 1/8640*(170*a - 293)))))*b4 + ((((1/172800*(3*a + 5)*r1 + 1/34560*(-3*a + 1))*t + (1/28800*a*r1 - 1/17280))*b2^2 + ((-1/8640*r1 + 1/17280*(3*a - 5))*t + (1/11520*(-a + 1)*r1 + 1/34560*(-3*a + 5)))*b2 + ((1/8640*(-27*a + 37)*r1 + 1/8640*(-201*a + 373))*t + (1/2880*(25*a - 47)*r1 + 1/8640*(93*a - 149))))*b3 + (((1/86400*(26*a - 47)*r1 + 1/17280*(-25*a + 44))*t + (1/172800*(-a - 3)*r1 + 1/34560*(17*a - 29)))*b2^2 + ((1/8640*(-13*a + 23)*r1 + 1/8640*(-25*a + 44))*t + (1/17280*(23*a - 39)*r1 + 1/4320*(15*a - 26)))*b2 + ((1/4320*(-30*a + 53)*r1 + 1/4320*(101*a - 178))*t + (1/960*(3*a - 5)*r1 + 1/8640*(-71*a + 133))))))*b1 + ((((1/14400*(-a + 2)*r1*t + (1/57600*(-5*a + 13)*r1 + 1/11520*(a - 3)))*b2^2 + (1/2880*(a - 1)*r1*t + 1/5760*(a - 1)*r1)*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(12*a - 25)*r1 + 1/576*(-8*a + 15))))*b3 + ((1/14400*(a - 1)*r1*t + (-1/28800*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + 1/2880*(2*a - 3)*r1)*b2 + (1/360*(3*a - 5)*r1*t + (1/2880*(-a + 1)*r1 + 1/576*(11*a - 19)))))*b4 + (((1/14400*(a - 2)*r1*t + (1/57600*(5*a - 13)*r1 + 1/11520*(-a + 3)))*b2^2 + (1/2880*(-a + 1)*r1*t + 1/5760*(-a + 1)*r1)*b2 + (1/720*(3*a - 4)*r1*t + (1/2880*(-9*a + 19)*r1 + 1/576*(5*a - 9))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/28800*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(2*a - 3)*r1*t + 1/2880*(-2*a + 3)*r1)*b2 + (1/720*(-9*a + 13)*r1*t + (1/1440*(8*a - 11)*r1 + 1/288*(-a + 2)))))))*i + (((((1/172800*(7*a - 11)*r1 + 1/34560*(23*a - 41))*t + (1/172800*(-41*a + 73)*r1 + 1/11520*(-3*a + 5)))*b2^2 + ((1/17280*(-7*a + 11)*r1 + 1/8640*(-11*a + 20))*t + (1/34560*(19*a - 35)*r1 + 1/3840*(3*a - 5)))*b2 + ((1/8640*(56*a - 97)*r1 + 1/8640*(-136*a + 235))*t + (1/3456*(11*a - 19)*r1 + 1/5760*(-3*a + 5))))*b3 + (((1/86400*(74*a - 127)*r1 + 1/17280*(53*a - 92))*t + (1/172800*(-169*a + 297)*r1 + 1/34560*(-93*a + 161)))*b2^2 + ((1/8640*(6*a - 11)*r1 + 1/8640*(-105*a + 181))*t + (1/4320*(5*a - 9)*r1 + 1/17280*(105*a - 181)))*b2 + ((1/8640*(-137*a + 233)*r1 + 1/8640*(-469*a + 821))*t + (1/2880*(56*a - 99)*r1 + 1/8640*(417*a - 718)))))*b4 + ((((1/172800*(-a + 5)*r1 + 1/34560*(-a + 3))*t + (1/172800*(3*a - 5)*r1 + 1/34560*(a - 1)))*b2^2 + ((1/17280*(a + 1)*r1 - 1/4320*a)*t + (1/34560*(3*a + 1)*r1 + 1/34560*(-a - 5)))*b2 + ((1/8640*(-29*a + 43)*r1 + 1/8640*(23*a - 39))*t + (-1/2160*r1 + 1/2160*(4*a - 7))))*b3 + (((1/86400*(-32*a + 55)*r1 + 1/17280*(-23*a + 40))*t + (1/172800*(67*a - 115)*r1 + 1/34560*(43*a - 75)))*b2^2 + ((1/8640*(-5*a + 8)*r1 + 1/4320*(27*a - 46))*t + (1/8640*(-5*a + 8)*r1 + 1/17280*(-49*a + 87)))*b2 + ((1/864*(6*a - 11)*r1 + 1/4320*(113*a - 200))*t + (1/8640*(-73*a + 115)*r1 + 1/2880*(-71*a + 125))))))*b1 + (((1/43200*(-4*a + 3)*r1*t + (1/172800*(-5*a - 3)*r1 + 1/34560*(-a + 3)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/4320*(-a + 3)*r1 + 1/3456*(a - 3)))*b2 + (1/1080*(20*a - 33)*r1*t + (1/8640*(185*a - 312)*r1 + 1/1728*(-35*a + 60))))*b3 + ((1/14400*(-3*a + 5)*r1*t + (1/28800*(8*a - 11)*r1 + 1/5760*(a - 2)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1440*(a - 2)*r1 + 1/576*(-a + 2)))*b2 + (1/180*(-4*a + 7)*r1*t + (1/2880*(-99*a + 169)*r1 + 1/576*(13*a - 23)))))*b4 + ((1/43200*(4*a - 3)*r1*t + (1/172800*(5*a + 3)*r1 + 1/34560*(a - 3)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/4320*(a - 3)*r1 + 1/3456*(-a + 3)))*b2 + (1/2160*(-22*a + 39)*r1*t + (1/8640*(-95*a + 177)*r1 + 1/1728*(17*a - 33))))*b3 + (1/14400*(3*a - 5)*r1*t + (1/28800*(-8*a + 11)*r1 + 1/5760*(-a + 2)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1440*(-a + 2)*r1 + 1/576*(a - 2)))*b2 + 1/720*(7*a - 13)*r1*t + 1/1440*(36*a - 65)*r1 + 1/288*(-5*a + 10), ((((((1/172800*(-27*a + 46)*r1 + 1/34560*(-5*a + 8))*t + (1/345600*(37*a - 66)*r1 + 1/69120*(19*a - 32)))*b2^2 + ((1/34560*(3*a - 5)*r1 + 1/34560*(71*a - 125))*t + (1/69120*(-23*a + 39)*r1 + 1/69120*(-91*a + 155)))*b2 + ((1/1728*(15*a - 26)*r1 + 1/8640*(-341*a + 590))*t + (1/17280*(5*a - 9)*r1 + 1/17280*(211*a - 365))))*b3 + (((1/172800*(-79*a + 143)*r1 + 1/34560*(-139*a + 239))*t + (1/345600*(539*a - 923)*r1 + 1/23040*(37*a - 65)))*b2^2 + ((1/17280*(39*a - 68)*r1 + 1/17280*(13*a - 22))*t + (1/34560*(-44*a + 77)*r1 + 1/34560*(-68*a + 117)))*b2 + ((1/1728*(7*a - 13)*r1 + 1/8640*(855*a - 1477))*t + (1/8640*(-271*a + 466)*r1 + 1/8640*(-320*a + 561)))))*b4 + ((((1/172800*a*r1 + 1/34560*(3*a - 2))*t + (1/345600*(-a + 10)*r1 + 1/69120*(-3*a + 4)))*b2^2 + ((1/34560*(a - 3)*r1 + 1/34560*(-3*a + 5))*t + (1/69120*(-a + 1)*r1 + 1/69120*(3*a + 5)))*b2 + ((1/8640*(-25*a + 42)*r1 + 1/8640*(153*a - 284))*t + (1/17280*(-17*a + 20)*r1 + 1/17280*(-123*a + 214))))*b3 + (((1/57600*(17*a - 29)*r1 + 1/34560*(79*a - 135))*t + (1/115200*(-97*a + 169)*r1 + 1/69120*(-71*a + 123)))*b2^2 + ((1/17280*(-17*a + 30)*r1 + 1/5760*(-5*a + 8))*t + (1/34560*(20*a - 33)*r1 + 1/34560*(40*a - 69)))*b2 + ((1/8640*(-43*a + 63)*r1 + 1/8640*(-499*a + 867))*t + (1/17280*(319*a - 567)*r1 + 1/17280*(439*a - 747))))))*b1 + ((((1/14400*(-a + 2)*r1*t + (1/57600*(4*a - 11)*r1 + 1/11520))*b2^2 + (1/2880*(a - 1)*r1*t + (1/11520*(-a + 1)*r1 + 1/2304*(-a + 1)))*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(-18*a + 35)*r1 + 1/576*(14*a - 25))))*b3 + ((1/14400*(-a + 1)*r1*t + (1/57600*(-a - 1)*r1 + 1/11520*(3*a - 5)))*b2^2 + (1/1440*(2*a - 3)*r1*t + (1/1920*(2*a - 3)*r1 + 1/1152*(-2*a + 3)))*b2 + (1/360*(-3*a + 5)*r1*t + (1/2880*(-7*a + 11)*r1 + 1/576*(17*a - 29)))))*b4 + (((1/14400*(a - 2)*r1*t + (1/57600*(-4*a + 11)*r1 - 1/11520))*b2^2 + (1/2880*(-a + 1)*r1*t + (1/11520*(a - 1)*r1 + 1/2304*(a - 1)))*b2 + (1/720*(3*a - 4)*r1*t + (1/2880*(12*a - 23)*r1 + 1/576*(-8*a + 13))))*b3 + ((1/14400*(a - 1)*r1*t + (1/57600*(a + 1)*r1 + 1/11520*(-3*a + 5)))*b2^2 + (1/1440*(-2*a + 3)*r1*t + (1/1920*(-2*a + 3)*r1 + 1/1152*(2*a - 3)))*b2 + (1/720*(9*a - 13)*r1*t + (1/576*(5*a - 7)*r1 + 1/576*(-11*a + 17)))))))*i + (((((1/172800*(4*a - 7)*r1 + 1/34560*(-22*a + 39))*t + (1/345600*(26*a - 43)*r1 + 1/69120*(24*a - 41)))*b2^2 + ((1/34560*(-17*a + 29)*r1 + 1/34560*(7*a - 15))*t + (1/69120*(17*a - 31)*r1 + 1/69120*(33*a - 55)))*b2 + ((1/8640*(86*a - 149)*r1 + 1/8640*(338*a - 585))*t + (1/17280*(-191*a + 331)*r1 + 1/17280*(-603*a + 1045))))*b3 + (((1/172800*(79*a - 137)*r1 + 1/34560*(-17*a + 27))*t + (1/345600*(31*a - 53)*r1 + 1/69120*(-21*a + 35)))*b2^2 + ((1/17280*(24*a - 43)*r1 + 1/17280*(250*a - 429))*t + (1/34560*(-183*a + 314)*r1 + 1/34560*(-215*a + 374)))*b2 + ((1/8640*(-131*a + 227)*r1 + 1/8640*(-167*a + 303))*t + (1/8640*(71*a - 122)*r1 + 1/8640*(226*a - 389)))))*b4 + ((((1/172800*(-2*a + 5)*r1 - 1/34560)*t + (1/345600*(2*a - 5)*r1 - 1/23040))*b2^2 + ((1/34560*(a - 1)*r1 + 1/34560*(9*a - 5))*t + (1/69120*(-a + 7)*r1 + 1/23040*(-3*a + 5)))*b2 + ((1/8640*(-40*a + 67)*r1 + 1/8640*(-114*a + 197))*t + (1/17280*(76*a - 133)*r1 + 1/1920*(26*a - 43))))*b3 + (((1/172800*(-47*a + 85)*r1 + 1/34560*(17*a - 31))*t + (1/345600*(-23*a + 45)*r1 + 1/69120*(5*a - 11)))*b2^2 + ((1/17280*(-10*a + 17)*r1 + 1/17280*(-116*a + 203))*t + (1/11520*(29*a - 50)*r1 + 1/34560*(91*a - 158)))*b2 + ((1/8640*(71*a - 133)*r1 + 1/8640*(3*a + 35))*t + (1/17280*(-41*a + 51)*r1 + 1/17280*(-153*a + 295))))))*b1 + (((1/43200*(-4*a + 3)*r1*t + (1/172800*(a + 6)*r1 + 1/34560*(5*a - 6)))*b2^2 + (1/8640*(-a + 3)*r1*t + (1/34560*(7*a - 21)*r1 + 1/6912*(-a + 3)))*b2 + (1/1080*(20*a - 33)*r1*t + (1/8640*(-145*a + 246)*r1 + 1/1728*(-5*a + 6))))*b3 + ((1/14400*(3*a - 5)*r1*t + (1/57600*(19*a - 27)*r1 + 1/11520*(-a + 1)))*b2^2 + (1/1440*(-a + 2)*r1*t + (1/1920*(a - 2)*r1 + 1/1152*(-a + 2)))*b2 + (1/180*(4*a - 7)*r1*t + (1/2880*(-83*a + 141)*r1 + 1/576*(-3*a + 5)))))*b4 + ((1/43200*(4*a - 3)*r1*t + (1/172800*(-a - 6)*r1 + 1/34560*(-5*a + 6)))*b2^2 + (1/8640*(a - 3)*r1*t + (1/34560*(-7*a + 21)*r1 + 1/6912*(a - 3)))*b2 + (1/2160*(-22*a + 39)*r1*t + (1/8640*(73*a - 138)*r1 + 1/1728*(5*a - 6))))*b3 + (1/14400*(-3*a + 5)*r1*t + (1/57600*(-19*a + 27)*r1 + 1/11520*(a - 1)))*b2^2 + (1/1440*(a - 2)*r1*t + (1/1920*(-a + 2)*r1 + 1/1152*(a - 2)))*b2 + 1/720*(-7*a + 13)*r1*t + 1/2880*(65*a - 117)*r1 + 1/576*(-3*a + 7) ];