Added new heterogeneity function, adapted from Tzeng YS, Albert A (2009) J. Appl. Physiol. Automated some section to speed up process and made algorithm more general. Some problems still exist in replicating results from original implementation