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 information on this web page is distributed under the MIT license.