MatLDA function

Fits the $J$-class penalized matrix-normal model for a single pair of tuning parameters.