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 ] 

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



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

Сейчас этот форум просматривают: нет зарегистрированных пользователей


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

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