The Neumann boundary conditions for Laplace's equation specify not the function φ itself on the boundary of D, but its normal derivative.
Solutions of Laplace's equation are called harmonic functions; they are all analytic within the domain where the equation is satisfied.