Здравствуйте!
Необходимо вычислить следующий интеграл:
Использую следующие команды:
>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 знаков, поэтому обычными методами у меня это не получилось сделать.