fit-MPHstar-method function

Fit method for mph class