regional flood frequency analysis using L-moment

I have Annual maximum series flood of last 30 years of 7 gauging stations , With the use of package of hosking (http://lib.stat.cmu.edu/general/lmoments) convert the whole code in matlab find the quantile. Basically I need two things 1. Heterogeneity measure with Monte Carlo simulation with 500 simulation. 2. Best Fit Distribution along L-moment ratio diagram