The Fourier coefficients satisfy and
Thus the square wave has Fourier series
It equals away from its jump discontinuities, and converges to at the jumps. The assigned endpoint values do not affect its Fourier coefficients.
For the Poisson equation, use separation of variables . The radial equation is . Its bounded solution with zero Dirichlet boundary condition at is
Indeed the negative power is excluded at the centre, and the boundary fixes the coefficient of . No resonance at occurs because only odd modes are forced. Therefore
This is a Poisson equation on a disk with angular forcing. The series converges uniformly on the closed disk: for each term is bounded by a constant times . Its first spatial derivatives also converge uniformly. For instance the radial derivative and times the angular derivative of each term are bounded by constants times ; the harmonic term is smooth in Cartesian coordinates. It defines a function with zero boundary data and the required weak formulation. Applying the Laplacian to finite partial sums gives the Fourier partial sums of , which converge in ; passing to the limit against smooth compactly supported test functions verifies the Poisson equation. It holds classically off the jump rays and the centre, where the forcing is locally constant. It cannot be interpreted as a global classical equation with continuous second derivatives for discontinuous forcing. Uniqueness follows from the zero-boundary Dirichlet energy identity for the difference of two solutions.
Finally orthogonality gives . The radial integral is
Uniform convergence of and boundedness of justify integration term by term, so
The negative sign also agrees with . The original PDF's radial denominator is ; the converted TeX's is a transcription error.