My approach was to set $\zeta + \zeta^6 = a_1x + a_2x^2 + a_3x^3$ and then solve for the constants.
I know there is a better way to solve this using the cyclotomic polynomials, but I don't understand how to do that. If anyone is willing to explain I would appreciate it.
You know that $\zeta^7=1$, so $\zeta$ is a root of $X^6+X^5+X^4+X^3+X^2+1$. Moreover $\zeta^6=\zeta^{-1}$. Consider $$ \zeta^6+\zeta^5+\zeta^4+\zeta^3+\zeta^2+\zeta+1=0 $$ and divide by $\zeta^3$: $$ \zeta^3+\frac{1}{\zeta^3}+\zeta^2+\frac{1}{\zeta^2}+\zeta+\frac{1}{\zeta}+1=0 $$ Observe that $$ x^3+y^3=(x+y)^3-3xy(x+y) $$ and $$ x^2+y^2=(x+y)^2-2xy $$ Use with $x=\zeta$ and $y=\zeta^{-1}$.