“…The first way is a statistical method which unfolds the primary energy spectra from the corresponding integral equation set based on a detected EAS data set and the model of the EAS development in the atmosphere [5][6][7][8][9]. The second method is based on an eventby-event evaluation [2,[10][11][12][13] of the primary energy of the detected EAS with parameters q ≡ q(N e , N µ , N h , s, θ) using parametric E = f (q) [2,10,11,13] or non-parametric [12] energy estimator previously determined on the basis of shower simulations in the framework of a given model of EAS development.…”