Необходимо найти распределение

вторичных частиц по углу их разлета в лабораторной системе в реакции

. Распределение необходимо найти аналитически и численно методом Монте-Карло. Аналитическое выражение получилось, а вот в моделировании возникли проблемы с определением функции плотности вероятности для угла. Нормальное распределение плохо согласуется с аналитической кривой. В монографии Бюклинга и Каянти приводится такое выражение для розыгрыша угла -

, которое также посредственно согласуется с аналитической кривой. Метод Неймана также не подошел. Как должна выглядеть функция распределения плотности вероятности в этом случае?