region, a FreeFem++ code which sets up an interesting mesh with four internal boxes where the solution has a jump, and writes the data to files, using the ffmatlib() interface, for subsequent graphics processing by MATLAB or Octave.
The computer code and data files described and made available on this web page are distributed under the MIT license