  In this paper, we propose a numerical method to approximate the solution of
partial differential equations in irregular domains with no-flux boundary
conditions by means of spectral methods. The main features of this method are
its capability to deal with domains of arbitrary shape and its easy
implementation via Fast Fourier Transform routines. We discuss several examples
of practical interest and test the results against exact solutions and standard
numerical methods.
