  We present a simple and efficient MATLAB implementation of the local
refinement for polygonal meshes. The purpose of this implementation is
primarily educational, especially on the teaching of adaptive virtual element
methods.
