2014 dxdy logo

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

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


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


В этом разделе нельзя создавать новые темы.

Если Вы хотите задать новый вопрос, то не дописывайте его в существующую тему, а создайте новую в корневом разделе "Помогите решить/разобраться (М)".

Если Вы зададите новый вопрос в существующей теме, то в случае нарушения оформления или других правил форума Ваше сообщение и все ответы на него могут быть удалены без предупреждения.

Не ищите на этом форуме халяву, правила запрещают участникам публиковать готовые решения стандартных учебных задач. Автор вопроса обязан привести свои попытки решения и указать конкретные затруднения.

Обязательно просмотрите тему Правила данного раздела, иначе Ваша тема может быть удалена или перемещена в Карантин, а Вы так и не узнаете, почему.



Начать новую тему Ответить на тему
 
 int( dt/(exp(t^2)+a) , t=0..infinity ) ?
Сообщение03.05.2014, 23:45 


22/12/11
66
В работе встретился такой параметрический интеграл:
$  \int_0^\infty \frac{dt}{\exp(t^2)+a}, \; a\ge 0$.

Замены переменного ничего не дали, также как и дифференцирование по параметру, также как и обращение к системам символьной алгебры (Maple).

Тем не менее, интеграл - хорошо сходящийся при всех a, при a=0 равен известному значению ($\sqrt{\pi}/2$), в общем, явно должен как-то выражаться через спецфункции... Но как - пока ума не приложу...

Может у кого-то будут мысли, как к нему подступиться ?
Может есть какие-то регулярные методы на сей счет, но не излагаемые в обычных институтских курсах ?

 Профиль  
                  
 
 Re: int( dt/(exp(t^2)+a) , t=0..infinity ) ?
Сообщение04.05.2014, 00:04 
Заслуженный участник


25/02/08
2961
По определению полилогарифма $\[L{i_s}(z) = \frac{1}{{\Gamma (s)}}\int\limits_0^\infty  {\frac{{{t^{s - 1}}}}{{\frac{{{e^t}}}{z} - 1}}} dt\]$. В вашем интеграле делаем замену $\[{t^2} = \xi \]$(приведя его к виду выше) и получим результат $\[ - \frac{{\sqrt \pi  }}{{2a}}{{\mathop{\rm Li}\nolimits} _{\frac{1}{2}}}( - a)\]$.
P.S.Кстати, я не знаю как вы у систем компьютерной алгебры спрашивали (видимо, просто неправильно "попросили"), но Mathematica даёт такой же ответ (и уж уверен, что и Maple справится).

 Профиль  
                  
 
 Re: int( dt/(exp(t^2)+a) , t=0..infinity ) ?
Сообщение04.05.2014, 15:55 


22/12/11
66
Ms-dos4 в сообщении #858718 писал(а):
По определению полилогарифма $\[L{i_s}(z) = \frac{1}{{\Gamma (s)}}\int\limits_0^\infty  {\frac{{{t^{s - 1}}}}{{\frac{{{e^t}}}{z} - 1}}} dt\]$. В вашем интеграле делаем замену $\[{t^2} = \xi \]$(приведя его к виду выше) и получим результат $\[ - \frac{{\sqrt \pi  }}{{2a}}{{\mathop{\rm Li}\nolimits} _{\frac{1}{2}}}( - a)\]$.
P.S.Кстати, я не знаю как вы у систем компьютерной алгебры спрашивали (видимо, просто неправильно "попросили"), но Mathematica даёт такой же ответ (и уж уверен, что и Maple справится).

Спасибо! Да, все просто оказалось. Кстати, Вы не перепутали в определении полилогарифма, там кажется должно быть $z e^t$ вместо $\frac{e^t}{z}$?

Что касается Maple, то как ни странно, функция полилогарифма в нем есть, и для целых s он вычисляет интегралы типа $\int_0^\infty u^n/(a\exp(u)+1)$, выражая их через полилогарифмы, но для дробных s - уже нет!

Более того, при попытке ввести определение полилогарифма в качестве интеграла (причем с ограничением, что |x|<1 и x - действительное число) - Maple (Maple 17) выдает некое страшное выражение + полилогарифм вместо самого полилогарифма!

P.S: интересно, а свободные CAS - системы типа maxima - справились бы с этим интегралом ? Спрашиваю, поскольку использую в работе Maple (для символьных вычислений и проверок некоторых численных подпрограмм на фортране) и Fortran (для численных вычислений), и если уж использовать что-то дополнительное, то хотелось бы, чтобы бесплатное и относительно легко изучаемое...

 Профиль  
                  
 
 Re: int( dt/(exp(t^2)+a) , t=0..infinity ) ?
Сообщение04.05.2014, 20:46 
Заслуженный участник


25/02/08
2961
bme
1)Не перепутал
2)Ставьте Mathematica 9

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

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



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

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


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

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