Let us consider the HEMs in the Figure by Telg and coworkers. Let us extract some data from the figure and start fitting to a BWF line and a Lorentzian line. Software Fityk (Wojdyr, 2010) is used. Fityk is a curve fitting and data analysis application, free and open source. It can perform a Levenberg Marquardt algorithm, among other best fit algorithms. The function of merit is the weighted sum of squared residuals (WSSR), also called chi square. Fityk has several built-in functions but allows to execute script with userdefined functions (see in the Appendix B how we can define the BWF in Fityk, and furthermore the q-Gaussian and q-BWF functions). Fig. : As in the Figure , here we use HEM data from the Figure by . The modes are fitted to q-BWF lines (we will discuss them in the last section of the review). The q-parameter is a parameter characterizing a generalization of the exponential function in the framework of the Costantino Tsallis statistics. Therefore, the parameter 1/q that we usually find in the definition of the BWF function is here substituted by -ξ. A small q-Gaussian is added to describe the shoulder of the G + peak. An even smaller q-Gaussian has been used in the high-frequency tail. The parameters q and ξ are given close to the respective q-BWF peaks.