Hi BackEMF this is my Matlab code you can use quad instead of quade
%clc; clear all
lamda=1.55e-6;
k=2*pi/lamda;
a=10000;
r=1e-2;L=1000;
f=@(x)(1./(1+a*x.^2).*exp(i*k/(2*L)*(x.^2-2*x*r)));%% integral by numerical solution
numerical=quade(f,-inf,inf)
j=sqrt(-1) and a is positive
I solved this integral numerically and i found the exact result in both Matlab and Mathematica program but I need analytical solution. I tried residue theorem but result didn't match numeric solutions.I asked some mathematicians but they couldn't find true path...