Получены физические и математические модели и каскадно-вероятностные функции генерации пи- и мю-мезонов в атмосфере Земли, модель излома в спектре космического излучения. Разработан численный алгоритм расчетов дифференциальных и интегральных энергетических спектров вторичных пи-мезонов в области высоких энергий (109 - 1015 эВ). Установлено, что излом в спектре первичного космического излучения связан с флуктуациями и ошибками эксперимента.