неправильно считается интеграл

Форум пользователей пакета Mathematica

Модератор: Admin

alcrux
Сообщения: 1
Зарегистрирован: Сб май 07, 2011 5:53 pm

неправильно считается интеграл

Сообщение alcrux » Сб май 07, 2011 6:04 pm

Добрый день!

Пытаюсь взять интеграл

NIntegrate[BesselJ[170.5, x]*Exp[I*x/0.997]/x, {x, 1, 180}]

выдается
NIntegrate::ncvb: "NIntegrate failed to converge to prescribed accuracy after 9 recursive bisections in x near {x} = {1.01016115099766266328951136443947689258493483066558837890625000000}. NIntegrate obtained 1.81148762123506880154106336439441599950157535305498542718495716529*10^340+2.88546721151785992842753671985635930330482884570433052937468115537*10^340\ I and 7.92198220817303350795617785505479590777849989125122434274125828887`65.954589770191*^339 for the integral and error estimates"

и соврешенно неправильный ответ:
1.811487621235069*10^340 + 2.885467211517860*10^340 I

если же вместо 172.5 взять 172.4, то все нормально.

Как с этим бороться?

dims12
Сообщения: 37
Зарегистрирован: Чт янв 30, 2003 1:23 am

Сообщение dims12 » Пн май 09, 2011 2:08 pm

А откуда известно, каков правильный ответ?

VSI
Сообщения: 305
Зарегистрирован: Вт мар 15, 2005 6:40 pm
Откуда: Мариуполь

Re: неправильно считается интеграл

Сообщение VSI » Пн май 09, 2011 5:08 pm

alcrux писал(а):Добрый день!

Пытаюсь взять интеграл

NIntegrate[BesselJ[170.5, x]*Exp[I*x/0.997]/x, {x, 1, 180}]

Как с этим бороться?


Решить в Mathcad'e...

http://i.exponenta.ru/exponenta/2011/5/9/5780.jpg

Vladimir__
Сообщения: 1660
Зарегистрирован: Вс дек 21, 2003 8:52 am

Re: неправильно считается интеграл

Сообщение Vladimir__ » Вт июн 14, 2011 10:39 am

alcrux писал(а):Добрый день!

Пытаюсь взять интеграл...

Как с этим бороться?

Сообщить об этой ошибке в техподдержку по продукту Mathematica компании Wolfram.

_kolos
Сообщения: 2
Зарегистрирован: Вт янв 25, 2011 12:20 pm

Re: неправильно считается интеграл

Сообщение _kolos » Пт июн 24, 2011 2:29 pm

[quote="Vladimir__"][quote="alcrux"]Добрый день!

Пытаюсь взять интеграл...

Как с этим бороться?[/quote]
Сообщить об этой ошибке в техподдержку по продукту Mathematica компании Wolfram.[/quote]

А посмотреть свойства функций Бесселя полуцелого порядка?