Consider Laplace's equation in polar coordinates $$\frac{\partial ^2u}{\partial r^2}+\frac{1}{r}\frac{\partial u}{\partial r}+\frac{1}{r^2}\frac{\partial ^2u}{\partial \theta^2}=0$$ for $0<r<1$. I'm given the boundary condition $u(1,\theta)=\cos^2(\theta)$.
I know that the general solution is given by $$u(r,\theta)=b_0+\sum_{n=1}^{\infty}r^n\left(a_n\sin (n\theta) +b_n\cos(n\theta) \right)$$ with \begin{align*} b_0&=\frac{1}{2\pi}\int_0^{2\pi}\cos^2(\theta) \, d\theta \\ b_n&=\frac{1}{\pi}\int_0^{2\pi} \cos^2(\theta)\cos(n\theta) \, d\theta \\ a_n & = \frac{1}{\pi}\int_0^{2\pi}\cos^2(\theta)\sin(n\theta) \, d\theta \end{align*}
It is clear to me that here we have $a_n=0$ and $b_0=\frac{1}{2}$ but I'm not sure how to find an explicit form for $b_n$.
HINT:
Recall that $\cos^2(\theta)=\frac{1+\cos(2\theta)}{2}$ and exploit orthogonality $\int_0^{2\pi}\cos(n\theta)\cos(m\theta)\,d\theta=0$ for $m\ne n$.