 Two-boundary annular harmonic interpolation

For a harmonic function on the planar annulus $a<r<b$, the radial amplitude of a sine or cosine Fourier mode of positive integer index $n$ is a linear combination of $r^n$ and $r^{-n}$. If its values at $a,b$ are $A_n,B_n$, respectively, the displayed hyperbolic-sine formula gives that combination. It follows either by solving the two boundary equations or by putting $s=\log (r/a)$, which changes the radial equation to $R_{ss}-n^{2}R=0$. The zero angular mode is instead affine in $\log r$:Superposing the modes gives the harmonic function with suitable prescribed Dirichlet boundary conditions. For a fixed interior radius, high angular modes from either boundary are exponentially suppressed. Piecewise smooth boundary data are attained at continuity points, while jumps have the usual mean-value Fourier limit.