  Within the MNPBEM toolbox, we show how to simulate electron energy loss
spectroscopy (EELS) of plasmonic nanoparticles using a boundary element method
approach. The methodology underlying our approach closely follows the concepts
developed by Garcia de Abajo and coworkers [for a review see Rev. Mod. Phys.
82, 209 (2010)]. We introduce two Matlab classes that allow in combination with
our recently developed MNPBEM toolbox for a simple, robust, and efficient
computation of EEL spectra and maps. The classes are accompanied by a number of
demo programs for EELS simulation of metallic nanospheres, nanodisks, and
nanotriangles, and for electron trajectories passing by or penetrating through
the metallic nanoparticles. We also discuss how to compute electric fields
induced by the electron beam and cathodoluminescence.
