fac := [16 - 30*t - 300*t^2 + 880*t^3 + Dt*(17*t - 56*t^2 - 304*t^3 + 1024*t^4) + Dt^2*(3*t^2 - 12*t^3 - 48*t^4 + 192*t^5), -448 - 420*t + 3840*t^2 + 3006*t^3 + 81660*t^4 + 115488*t^5 + 442880*t^6 + 49152*t^7 + Dt*(-322*t - 420*t^2 + 1503*t^3 - 1512*t^4 + 30276*t^5 + 156672*t^6 + 812864*t^7 + 122880*t^8) + Dt^2*(-42*t^2 - 63*t^3 - 51*t^4 - 900*t^5 - 156*t^6 + 28224*t^7 + 187584*t^8 + 36864*t^9), 14 - 84*t - 156*t^2 + 1040*t^3 + Dt*(16*t - 76*t^2 - 224*t^3 + 1088*t^4) + Dt^2*(3*t^2 - 12*t^3 - 48*t^4 + 192*t^5), -80 + 220*t + 3616*t^2 - 17652*t^3 + 4008*t^4 + 292704*t^5 - 93152*t^6 - 8960*t^7 + 409600*t^8 + Dt*(-50*t + 70*t^2 + 2728*t^3 - 15168*t^4 - 15720*t^5 + 254832*t^6 - 73664*t^7 - 47360*t^8 + 901120*t^9) + Dt^2*(-6*t^2 + 3*t^3 + 384*t^4 - 2382*t^5 - 4392*t^6 + 39624*t^7 - 18816*t^8 - 36480*t^9 + 245760*t^10), 2 + Dt*t, -8 + 28*t + 112*t^2 - 160*t^3 + 320*t^4 + 512*t^5 + Dt*(-7*t + 20*t^2 + 96*t^3 - 256*t^4 + 128*t^5 + 1024*t^6) + Dt^2*(-t^2 + 2*t^3 + 16*t^4 - 48*t^5 + 256*t^7)];