A MATLAB function that fits the specified probability distribution function to mobile data consumption for using Monte Carlo methods for predicting which data plan is most likely to result in the lowest future cost. Applicable to single or shared data plans.