I have the Laplace equation $$\nabla^2 \phi(x,y) = 0$$ where $\phi$ is a function defined on the domain $\mathbb{R^2}\setminus(C_1 \cup C_2)$.
$C_1$ and $C_2$ are two circles of radius $r = 0.5$ and whose centers are $(0,1)$ and $(0,-1)$.
The boundary conditions are $\phi = 0$ on the bottom circumference, $\phi = 1$ on the upper circumference, and $\phi \rightarrow 0.5$ at infinity.
For reference, this problem is related to the calculation of the potential field in a two-wire transmission line. My question is the following: is it solvable analytically and, if so, how does one do it? I suppose this necessitates some sort of expansion into orthogonal functions, but what coordinate system would be the most appropriate?
Hint: The function $$ \phi(x,y)=\frac{1}{2}+k\ln\left(\frac{x^2+(y+a)^2}{x^2+(y-a)^2}\right) $$ is a solution to the Laplace equation satisfying $\phi\to\frac{1}{2}$ at infinity and such that the level curves $\phi(x,y)=c$ are circumferences with center on the $y$-axis. To solve your problem, find $a$ and $k$ such that $\phi=1$ on the upper circumference (by symmetry, $\phi=0$ on the bottom circumference).