How to find the parameter values of Hyper-Erlang distribution in Matlab

Hi,

MLE and fitdist functions in Matlab gives the parameters values of given data set from various distributions such as exponential, gamma etc. But, its not applicable for Erlang distribution. I want to get the parameter values of Hyper-Erlang distribution. Kindly help me to know the method for finding the parameter values of Hyper-Erlang distribution.