Mathematical modeling of an experimental facility for the study of the primary energy spectrum of cosmic rays at energies of 10 to the 15th to 10 to the 20th eV