2014 dxdy logo

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

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




Начать новую тему Ответить на тему
 
 Вычисление интеграла с заданной точностью в Maple
Сообщение06.03.2012, 15:19 


17/10/10
49
Здравствуйте!

Необходимо вычислить следующий интеграл:
$f(x) = \int_0^\infty e^{-v} \cos(xv) \sin(v^2) dv$

Использую следующие команды:

>integ:=Int(exp(-p*v)*sin(v^2), v=0..infinity);
>L:= laplace(exp(-p*v)*sin(v^2), v, p); # преобразование Лапласа для этого интеграла

>evalc(subs(p=1+I*x, integ)); # такой подстановкой можно получить в одном из двух слагаемых мой интеграл

тогда, делаю такую же подстановку и беру действительную часть:
>expr:=Re(subs(p = 1+I*x, L));

Но печаль в том, что если посчитать, например, в нуле первоначальный интеграл и полученный после преобразования, то ответы не совпадут.

В чем здесь может быть ошибка?

Посчитать надо с точностью до 11 знаков, поэтому обычными методами у меня это не получилось сделать.

 Профиль  
                  
 
 Re: Вычисление интеграла с заданной точностью в Maple
Сообщение11.03.2012, 18:13 
Заморожен


14/09/10
72
Следует набирать не >L:= laplace(exp(-p*v)*sin(v^2), v, p);, а >L:= laplace (sin(v^2), v, p);
По поводу числа знаков см. переменную Digits (Digits:= 12), параметр epsilon evalf/Int и параметр функции evalf.

 Профиль  
                  
 
 Re: Вычисление интеграла с заданной точностью в Maple
Сообщение12.03.2012, 20:19 


17/10/10
49
Спасибо, Andrew Gubarev. Да, косяк с функцией laplace понятен и устранен. Что же касается Digits, то не вычисляется ничего при Digit > 10, пишет в виде интегралов. Поэтому преобразование Лапласа и используется.

 Профиль  
                  
 
 Re: Вычисление интеграла с заданной точностью в Maple
Сообщение12.03.2012, 21:03 
Заморожен


14/09/10
72
Если нужна помощь, то приведите код (с комментариями), который не выполняется. Я постараюсь помочь. (Но придется набраться терпения. Я сразу не отвечу, если это будет не выходной.)

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

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



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

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


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

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