<div>The code, fibriprogram_v5.m, manages the optimisation process for fitting the gaussian curves to the frequency versus D dataset by calling on sub-routines generation of a gaussian function (gaussiansingle_v5.m), linear combination of the gaussian functions (gaussiancomposite_v5.m), SA (simulanneal_v5.m) and the metropolis rule (metrop_v1.m). </div><div>Similarly, software Matlab version 6.1 was used for running the code.</div