Solve

$\displaystyle \bigtriangledown ^2 G = \delta (x - x_0)$

in the first quadrant $\displaystyle (x \geq 0, y \geq 0)$ with $\displaystyle G = 0$ on the boundaries in 2 dimensions.

Now I know the solution to this problem for the half-plane$\displaystyle (y \geq 0, - \infty < x < \infty)$ is determined by the fact that

$\displaystyle G = \frac{1}{2 \pi} \ln |r - r_0|$

where $\displaystyle r = \sqrt{(x - x_0)^2 + (y - y_0)^2}$

and $\displaystyle r_0 = \sqrt{(x - x_0)^2 + (y + y_0)^2}$

which leads to

$\displaystyle G = \frac{1}{4 \pi} \ln \Big((x - x_0)^2 + (y - y_0)^2 \Big) + \frac{1}{4 \pi} \ln \Big( (x - x_0)^2 + (y + y_0)^2 \Big)$

So to solve this problem, however, I think I need to add more image sources, so at a guess I would have...

$\displaystyle r = \sqrt{(x - x_0)^2 + (y - y_0)^2}$

$\displaystyle r_1 = \sqrt{(x + x_0)^2 + (y - y_0)^2}$

$\displaystyle r_2 = \sqrt{(x - x_0)^2 + (y + y_0)^2}$

$\displaystyle r_3 = \sqrt{(x + x_0)^2 + (y + y_0)^2}$

and then due to superposition, Green's function would be given by

$\displaystyle G = \frac{1}{2 \pi} \ln |r - r_0 - r_1 - r_3|$

or something to that effect. Am I on the right track with this?