Add Matlab implementation of Hughes-Hartoggs