I'm having trouble finding out how to set up Neumann (or, rather, "Robin") boundary conditions for a diffusion-type PDE. More specifically, I have a scalar function f(\boldsymbol{x}, t) where \boldsymbol{x} is n-dimensional vector space with some boundary region defined by A(\boldsymbol{x})=0...