Здравствуйте. Пытаюсь вычислить следующее:
![$$
N_\mu=
\int\limits_{E_{\min}}^{\infty} dE_\mu
\int\limits_{0}^{h_0} dh_2
\int\limits_{0}^{h_2} dh_1
\int\limits_{E_\mu}^{\frac{E_\mu}{\alpha^2}} dE_\pi
\; F(E_\pi, E_\mu, h_1, h_2, h_0, \theta)
$$ $$
N_\mu=
\int\limits_{E_{\min}}^{\infty} dE_\mu
\int\limits_{0}^{h_0} dh_2
\int\limits_{0}^{h_2} dh_1
\int\limits_{E_\mu}^{\frac{E_\mu}{\alpha^2}} dE_\pi
\; F(E_\pi, E_\mu, h_1, h_2, h_0, \theta)
$$](https://dxdy-03.korotkov.co.uk/f/2/5/1/251ca8a6c94bbddc8be1f5d46251165c82.png)
Где
![$$
F(E_\pi, E_\mu, h_1, h_2, h_0, \theta)=
\frac{m_\pi\;c}{\tau_\pi\;\cos\theta}\;
\frac{\alpha^2}{(1-\alpha^2)\;E_\mu}\;
\frac{E_\pi^{-\gamma-1}}{\rho(h_2)}\;
\times
$$ $$
F(E_\pi, E_\mu, h_1, h_2, h_0, \theta)=
\frac{m_\pi\;c}{\tau_\pi\;\cos\theta}\;
\frac{\alpha^2}{(1-\alpha^2)\;E_\mu}\;
\frac{E_\pi^{-\gamma-1}}{\rho(h_2)}\;
\times
$$](https://dxdy-01.korotkov.co.uk/f/8/3/7/8376b89796f9b76fa24c36093e77ab8c82.png)
![$$
\times Exp \left(
-\frac{h_1}{L_p\;\cos\theta}\;
-\frac{(h_2-h_1)}{l_\pi\;\cos\theta}
-\frac{m_\pi\;c}{\tau_\pi\; E_\pi\; \cos\theta} \int\limits_{h_1}^{h_2}\frac{dh}{\rho(h)}\;
-\frac{m_\mu\;c}{\tau_\mu} \int\limits_{h_0}^{h_2}\frac{dh}{\rho(h)\;(E_\mu\; \cos\theta)-a(h-h_2)}
\right)
$$ $$
\times Exp \left(
-\frac{h_1}{L_p\;\cos\theta}\;
-\frac{(h_2-h_1)}{l_\pi\;\cos\theta}
-\frac{m_\pi\;c}{\tau_\pi\; E_\pi\; \cos\theta} \int\limits_{h_1}^{h_2}\frac{dh}{\rho(h)}\;
-\frac{m_\mu\;c}{\tau_\mu} \int\limits_{h_0}^{h_2}\frac{dh}{\rho(h)\;(E_\mu\; \cos\theta)-a(h-h_2)}
\right)
$$](https://dxdy-03.korotkov.co.uk/f/6/5/9/659245ec0db63fbec718b591dc21c32a82.png)
![$$
\rho(h)=
\frac{h\;g}{Tp(h)\;R_0}
$$ $$
\rho(h)=
\frac{h\;g}{Tp(h)\;R_0}
$$](https://dxdy-02.korotkov.co.uk/f/9/0/6/906e8a40c9bf4c6c1892e230b7f1f12a82.png)
![$$
Tp(h)=228.5182 -0.0768 h+0.0001h^2
$$ $$
Tp(h)=228.5182 -0.0768 h+0.0001h^2
$$](https://dxdy-03.korotkov.co.uk/f/2/7/6/2763873f95a6e5b2f329e1dd6cb5f7f682.png)
Константы
![$$
h_0, E_{\min} , m_\mu ,m_\pi, \alpha, \gamma, L_p, l_\pi, c, \tau_\pi , \tau_\mu, g , R_0 , a , \theta
$$ $$
h_0, E_{\min} , m_\mu ,m_\pi, \alpha, \gamma, L_p, l_\pi, c, \tau_\pi , \tau_\mu, g , R_0 , a , \theta
$$](https://dxdy-02.korotkov.co.uk/f/5/f/a/5fac41e08e51a3bca2870c48ff51b19882.png)
Считаю в Wolfram Mathematica. Проблема в том, что при некоторых значениях придела Emin(конкретно 0.77) интеграл получается неразумным, отрицательным или равным нулю. С Wolfram Mathematica и численным интегрированием знаком не так давно. Пробовал считать по-разному. Напрямую, сначала символьно посчитать два последних интеграла, затем это выражение проинтегрировать по h2 с шагом по h0, превратить то что получилось в сплайн-функцию и затем проинтегрировать и ее по Emu с шагом по E. Не помогло. Перебирал различные правила и стратегии в Wolfram Mathematica(может не все), что так же не привело к разумным результатам. К тому же не ясно откуда в решении иногда появляется мнимая составляющая(хотя она есть и в примерах по интегрированию у Вольфрама).
Пожалуйста, подскажите в каком направлении двигаться по данному вопросу.