A MatLab routine to extract Principal Components of Discrimination