Normal Derivative values at boundary

See my solution in this post.