trigsimp takes forever on trig expressions with high degree polynomials

[meganly]

>>> from sympy import trigsimp, fu, symbols, sin
>>> x = symbols('x')
>>> trigsimp(sin((x+1)**9))

never finishes whereas

>>> fu(sin((x+1)**9))
sin((x+1)**9)

evaluates instantly. Is trigsimp is getting stuck trying to apply the trig sum formula?

[asmeurer]

It's stuck trying to factor this enormous expression

(4*(1 - 2*sin(9*x**7)**2)*(-4*(1 - 2*sin(21*x**6)**2)*((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) + 2*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*sin(21*x**6)*cos(21*x**6) + ((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) - 2*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*((1 - 2*sin(21*x**6)**2)**2 - 4*(1 - sin(21*x**6)**2)*sin(21*x**6)**2))*sin(9*x**7)*cos(9*x**7) + (4*(1 - 2*sin(21*x**6)**2)*((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) - 2*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*sin(21*x**6)*cos(21*x**6) + ((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) + 2*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*((1 - 2*sin(21*x**6)**2)**2 - 4*(1 - sin(21*x**6)**2)*sin(21*x**6)**2))*((1 - 2*sin(9*x**7)**2)**2 - 4*(1 - sin(9*x**7)**2)*sin(9*x**7)**2))*cos(9*x**8) + (-4*(1 - 2*sin(9*x**7)**2)*(4*(1 - 2*sin(21*x**6)**2)*((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) - 2*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*sin(21*x**6)*cos(21*x**6) + ((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) + 2*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*((1 - 2*sin(21*x**6)**2)**2 - 4*(1 - sin(21*x**6)**2)*sin(21*x**6)**2))*sin(9*x**7)*cos(9*x**7) + (-4*(1 - 2*sin(21*x**6)**2)*((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) + 2*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*sin(21*x**6)*cos(21*x**6) + ((1 - 2*sin(63*x**5)**2)*((1 - 2*sin(63*x**4)**2)*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) - 2*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4)) - 2*((1 - 2*sin(63*x**4)**2)*(4*(1 - 2*sin(21*x**3)**2)*(-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2)) + 2*(-4*(1 - 2*sin(21*x**3)**2)*(4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*sin(21*x**3)*cos(21*x**3) + (-4*(1 - 2*sin(9*x**2)**2)*((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*sin(x**9) + (sin(9*x)*cos(1) + sin(1)*cos(9*x))*cos(x**9))*sin(9*x**2)*cos(9*x**2) + ((-sin(1)*sin(9*x) + cos(1)*cos(9*x))*cos(x**9) - (sin(9*x)*cos(1) + sin(1)*cos(9*x))*sin(x**9))*((1 - 2*sin(9*x**2)**2)**2 - 4*(1 - sin(9*x**2)**2)*sin(9*x**2)**2))*((1 - 2*sin(21*x**3)**2)**2 - 4*(1 - sin(21*x**3)**2)*sin(21*x**3)**2))*sin(63*x**4)*cos(63*x**4))*sin(63*x**5)*cos(63*x**5))*((1 - 2*sin(21*x**6)**2)**2 - 4*(1 - sin(21*x**6)**2)*sin(21*x**6)**2))*((1 - 2*sin(9*x**7)**2)**2 - 4*(1 - sin(9*x**7)**2)*sin(9*x**7)**2))*sin(9*x**8)

I had thought that it was fu, not trigsimp, that had this issue.

[asmeurer]

I think this is the same issue as #7705, though with a much simpler repro.

[meganly]

Would including fu as part of the default in trigsimp, similar to how combined uses matching and groebner, help? In this case and other cases like #17633 fu does the right thing.

[asmeurer]

I'm not clear why fu isn't the default in trigsimp. I think there is a good reason for it, but I don't know if it is written down somewhere.

[meganly]

There are times when the current default (matching) works better than fu:

>>> trigsimp(sec(x)**2-1)
tan(x)**2
>>> fu(sec(x)**2-1)
sec(x)**2-1

It seems to me that trigsimp should try both matching and fu, and take the better of the two results.