  In \cite{AcostaBorthagaray}, a complete $n$-dimensional finite element
analysis of the homogeneous Dirichlet problem associated to a fractional
Laplacian was presented. Here we provide a comprehensive and simple 2D {\it
MATLAB}\textsuperscript{\textregistered} finite element code for such a
problem. The code is accompanied with a basic discussion of the theory relevant
in the context. The main program is written in about 80 lines and can be easily
modified to deal with other kernels as well as with time dependent problems.
The present work fills a gap by providing an input for a large number of
mathematicians and scientists interested in numerical approximations of
solutions of a large variety of problems involving nonlocal phenomena in
two-dimensional space.
