I have the following Neumann problem, and I don't know what I'm missing:
Let $\Omega = D(0,1)\setminus\{0\}$. I want to find all $f \in C^2(\Omega)$, such that
$$\Delta_{\Bbb R^2} f = \frac{\partial^2 f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2} = r \cos\varphi = g(r,\varphi) \qquad (x,y) \in \Omega$$
with
$$\frac{\partial f}{\partial r} (1,\varphi) = 0.$$
So, what I did was to convert the Laplacian to polar coordinates and decompose $f(r,\varphi) = \sum_{n\in\Bbb Z} f_n(r) \cdot e^{in\varphi}$ as well as $g(r,\varphi) = \sum_{n\in\Bbb Z} g_n(r) \cdot e^{in\varphi}$ into its respective fourier series in order to get an equation for the coefficients $f_n(r)$ and $g_n(r)$. Obviously $g_n = 0$ except $g_{\pm 1} = \frac{r}{2}$. So I get
\begin{align} \text{for} \quad n \in \{-1,1\}&: \qquad \frac{1}{r} \frac{d}{dr} \left( r \frac{df_{\pm 1}(r)}{dr} \right) - \frac{1}{r^2} f_{\pm 1}(r) = \frac{r}{2} \\ \text{for} \quad n \notin \{-1,1\}, \,\, n \in \Bbb Z&: \qquad \frac{1}{r} \frac{d}{dr} \left( r \frac{df_{n}(r)}{dr} \right) - \frac{n^2}{r^2} f_{n}(r) = 0 \end{align}
Using the ansatz $f_n = c_n r^\alpha$, where $\alpha \in \Bbb R$, I get
\begin{align} \text{for} \quad n \in \{-1,1\}&: \qquad f_{\pm 1}(r) = \frac{1}{16} r^3 \\ \text{for} \quad n \notin \{-1,1\}, \,\, n \in \Bbb Z&: \qquad f_n(r) = c_n r^n \end{align}
and hence
$$f(r,\varphi) = \sum_{n\in\Bbb Z} f_n(r) \cdot e^{in\varphi} = \frac{1}{8} r^3 \cos \varphi + \sum_{n\in\Bbb Z \setminus \{-1,1\}} c_n r^n e^{in\varphi} \overset{\color{red}{?}}{=} \frac{1}{8} r^3 \cos \varphi + \frac{a_0}{2} + \sum_{\color{red}{n \geq 2}} r^n (A_n \cos n\varphi + B_n \sin n\varphi) \qquad a_0 = \text{const.}$$
The Neumann condition yields
$$\frac{\partial f}{\partial r} (1,\varphi) = {3 \over 8} \cos \varphi + \sum_{\color{red}{n \geq 2}} n (A_n \cos n\varphi + B_n \sin n\varphi) \overset{!}{=} 0$$
The crucial part is the summation index $n$. If I would have $n \geq 1$ it would immediately follow that $B_n = 0$ and $A_n = 0 \,\, \forall n$ except $A_1 = - {3 \over 8}$, which would imply
$$f(r,\varphi) = \frac{1}{8} r^3 \cos \varphi + \frac{a_0}{2} - \frac{3}{8} r \cos \varphi \qquad a_0 = \text{const.},$$
which solves the above stated problem. But is this really the only $f \in C^2(\Omega)$ solving the problem? With $n \geq 2$ the Neumann condition seems to have no solution. What am I missing? I suspect, that it should probably be $n \geq 1$, but I don't see how to get there. Or am I totally wrong? Any help would be appreciated.