Hadrons serve as the "skeleton" of Extensive Air Shower (EAS), carrying critical information about cosmic ray composition and energy. Electron-Neutron Detector Array (ENDA) is capable of detecting the electromagnetic components generated by EAS as well as the secondary thermal
neutrons produced by hadrons. This report simulates the response of the ENDA detector and conducts an investigation into the lateral distribution of thermal neutrons using simulated data. By selecting various parameters, the fitting performance of different functional models for the thermal neutron’s lateral distribution is evaluated. The results indicate that a double exponential function exhibits superior fitting performance for the lateral distribution of thermal neutrons. These findings provide a basis for optimizing the modeling approach for thermal neutron lateral distributions.

