2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


Посмотреть правила форума



Начать новую тему Ответить на тему
 
 Численное интегрирование нетривиальной функции
Сообщение17.06.2016, 08:40 


24/05/16
2
Здравствуйте. Пытаюсь вычислить следующее:


$$
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)
$$

Где
$$ 
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 
$$

$$ 
\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)
$$

$$ 
\rho(h)=
\frac{h\;g}{Tp(h)\;R_0}
$$

$$ 
Tp(h)=228.5182 -0.0768 h+0.0001h^2
$$


Константы
$$ 
h_0, E_{\min} , m_\mu ,m_\pi, \alpha, \gamma, L_p, l_\pi, c, \tau_\pi , \tau_\mu, g , R_0 , a , \theta
$$


Считаю в Wolfram Mathematica. Проблема в том, что при некоторых значениях придела Emin(конкретно 0.77) интеграл получается неразумным, отрицательным или равным нулю. С Wolfram Mathematica и численным интегрированием знаком не так давно. Пробовал считать по-разному. Напрямую, сначала символьно посчитать два последних интеграла, затем это выражение проинтегрировать по h2 с шагом по h0, превратить то что получилось в сплайн-функцию и затем проинтегрировать и ее по Emu с шагом по E. Не помогло. Перебирал различные правила и стратегии в Wolfram Mathematica(может не все), что так же не привело к разумным результатам. К тому же не ясно откуда в решении иногда появляется мнимая составляющая(хотя она есть и в примерах по интегрированию у Вольфрама).
Пожалуйста, подскажите в каком направлении двигаться по данному вопросу.

 Профиль  
                  
 
 Re: Численное интегрирование нетривиальной функции
Сообщение17.06.2016, 23:17 
Заслуженный участник
Аватара пользователя


15/10/08
12801
Направление: наймите специалиста.

P.S. Потому как задача явно на пару порядков превосходит ваши возможности.

P.P.S. Я не специалист.

 Профиль  
                  
 
 Re: Численное интегрирование нетривиальной функции
Сообщение22.06.2016, 06:43 


24/05/16
2
Подскажите, как найти специалиста который может помочь (естественно не безвозмездно).

 Профиль  
                  
 
 Re: Численное интегрирование нетривиальной функции
Сообщение22.06.2016, 11:57 
Заслуженный участник


25/02/11
1803
Раз уж вы в математике считали, разместили бы полный код для NIntegrate с каким-нибудь набором параметров, для которого интеграл получается неразумным.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 4 ] 

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: mihaild


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group