2(2cos(4x)+1)cos(x)=1,
2cos(x)(2(2cos2(2x)-1)+1)=1,
2cos(x)(4cos2(2x)-1)=1,
2cos(x)(2cos(2x)-1)(2cos(2x)+1)=1,
2cos(x)(2(2cos2(x)-1)-1)(2(2cos2(x)-1)+1)=1,
2cos(x)(4cos2(x)-3)(4cos2(x)-1)=1.
Let y=2cos(x) then y(y2-3)(y2-1)=1,
y5-4y3+3y-1=0=(y2+y-1)(y3-y2-2y+1).
Let f(y)=y5-4y3+3y-1, f'(y)=5y4-12y2+3 then, using Newton's iterative method:
yn+1=yn-f(yn)/f'(yn).
Let y0=0, then y1=-(-1)/3=⅓, y2=5/12, ..., y=2cos(x)=0.4450418679, so x=±3π/7+2πn, where n is an integer.
Let y0=1, then y1=¾, y2=73/111, ..., y=2cos(x)=0.6180339888, so x=±2π/5+2πn.
Let y0=2, then y1=13/7, ..., y=2cos(x)=1.801937736, so x=±π/7+2πn.
Let y0=-1, then y1=-5/4, ..., y=2cos(x)=-1.246979604, so x=±5π/7+2πn.
Let y0=-2, then y1=-9/5, ..., y=2cos(x)=-1.618033989, so x=±4π/5+2πn.
The substitutions for y0 appear to cover all 5 roots of the quintic equation, implying that the 5 generic solutions for x (in radians) represent all the solutions for x. The solutions are all rational fractions of π, so this suggests that another method is likely to exist which gives the solutions without involving solving equations of higher degree than 2. Two of the solutions for y involve the Golden Ratio which are the roots of the quadratic y2+y-1=0 (see above) and the other three solutions for y are roots of the cubic y3-y2-2y+1=0.
From the above factorisation of the quintic we can also write:
y2+y-1=0, that is, 4cos2(x)+2cos(x)-1=0, 2(cos(2x)+1)+2cos(x)-1=0,
2cos(2x)+2cos(x)+1=0, or 2cos(2x)+2cos(x)=-1.
y=(-1±√5)/2, hence cos(x)=(-1±√5)/4,
so x=2π/5+2πn or 2π/5+2πn and 2cos(4π/5)+2cos(2π/5)+1=0, or 2cos(8π/5)+2cos(4π/5)+1=0.
y3-y2-2y+1=0, that is, 8cos3(x)-4cos2(x)-4cos(x)+1=0,
4cos(x)(2cos2(x)-1)-(4cos2(x)-2+1)=0,
4cos(x)cos(2x)-(2cos(2x)+1)=0,
4cos(x)cos(2x)-2cos(2x)-1=0, or 2cos(2x)(2cos(x)-1)=1.
So 2cos(2π/7)(2cos(π/7)-1)=1, 2cos(6π/7)(2cos(3π/7)-1)=1, 2cos(10π/7)(2cos(5π/7)-1)=1.