Hello I'm trying to integrate the following function in MATLAB
but it's returing the wrong answer when I try something like
This is what I have tried so far:
fun = @(t,x,y) exp(1i.*(t.^4+x.*t.^2+y.*t));
P = @(x,y) integral(@(t)fun(t,x,y),-Inf,Inf);
P(1,1)
Any help appreciated and many thanks in advance

3 comentarios

madhan ravi
madhan ravi el 14 de Feb. de 2019
what’s the exact answer?
Torsten
Torsten el 15 de Feb. de 2019
exp(i*(t^4+x*t^2+y*t)) does not tend to 0 as | t| -> Inf. Thus your integral does not exist (at least in the usual sense).
Michael Devereux
Michael Devereux el 15 de Feb. de 2019
According to WolframAlpha the answer is 1.20759 + 0.601534 i
Keep in mind it's a complex exponential so there is a finite solution. This is know as the Pearcey Integral. I am more concerned that I have entered the formula incorrectly than the actual integral itself. Is this the correct way to approach the problem.

Iniciar sesión para comentar.

 Respuesta aceptada

Torsten
Torsten el 15 de Feb. de 2019

0 votos

format long
fun = @(t,x,y) exp(-t.^4 + 1i.*y.*t - x.*t.^2 + 1i*pi*0.125);
P = @(x,y) integral(@(t)fun(t,x*exp(-1i*pi*0.25),y*exp(1i*pi*0.125)),-Inf,Inf);
P(1,1)
Reference:
https://arxiv.org/pdf/1601.03615.pdf

Más respuestas (1)

Abhishek Hullur
Abhishek Hullur el 8 de Ag. de 2021

0 votos

. Evaluate around the rectangle with vertices

Categorías

Más información sobre Programming en Centro de ayuda y File Exchange.

Productos

Versión

R2018b

Etiquetas

Preguntada:

el 14 de Feb. de 2019

Respondida:

el 8 de Ag. de 2021

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by