  In recent years we provided numerical methods based on pseudospectral
collocation for computing the Floquet multipliers of different types of delay
equations, with the goal of studying the stability of their periodic solutions.
The latest work of the series concerns the extension of these methods to a
piecewise approach, in order to take the properties of numerically computed
solutions into account. In this chapter we describe the MATLAB implementation
of this method and provide practical usage examples.
