Hi, I'll post the question and how far I got with it.

Using the method of images, construct the Green's function of the Neumann problem for Laplace's equation in the half-space . Hence solve:

and evaluate the resulting integral on the z-axis.

=========

So the first thing I did was to write down the fundamental solution to Laplace's equation in 3D, and take off the image points in the lower half-space (z<0). i.e.

Now to solve, we use Green's identity.

Integrating over the upper half-space leaves just on one side (since vanishes for ).

For the other side, we can write since the outward normal is going in the opposite direction to increasing z.

We can also say at (as it's a Neumann problem..?)

We now evaluate this resulting integral on the z-axis to get:

I get stuck at this point. In the solutions, the above has been transformed into polars. So on the z-axis, G becomes (leaving the constant outside the integral),

=

And the goes to , as usual

The next line of the solution is:

I don't understand what's happened to the - is it equal to 1? Something to do with the boundary conditions maybe?

After this step, it looks pretty straightforward, just normal integration.

Any help is much appreciated!

Thanks.