Illustration of m-code for estimating an ELD parametrization for univariate Monte Carlo data. The computer program works by loading the Monte Carlo univariate data as an input ydata = [y₁, … , y_M] ^T where M = 10 000 is the number of Monte Carlo simulation events and it then uses this data to determine the extended lambda distribution parameters a, b, c and d so that the probability density function can be analytically specified.