How do we numerically integrate a rapidly decaying exponential function? A simple Gauss quadrature which is based on approximating the function by polynomial, I think will not work, since rapidly decaying exponential function will require more terms in the series. Is there a series approximation for the function?
Numerical integration - Gauss quadrature rule
2
$\begingroup$
integration
numerical-methods
approximation
-
0Did you try the variants of the basic Gauss quadrature like Gauss–Laguerre quadrature and Gauss–Hermite quadrature? – 2010-12-24
-
0I haven't applied them yet but aren't they of different form? Will they work for a rapidly decaying function? I have the integral of the form $$\int f(x)e^{-kx}$$, where k lies in the range of [10,50] – 2010-12-24
-
0@vijay: is $f(x)$ smooth, or does it have singular behavior (e.g. at the left endpoint)? – 2010-12-24
-
0@J.M. f(x) is smooth. Ok, I think I will be more clear here. The actual integration to be done goes like this,$$\int (\int f(v) e^{-kvx}dv \int g(v) e^{-kvx}dv)dx$$. Integrands with $f(v)$ & $g(v)$ can only be numerically evaluated. They have 10th order polynominal in $v$ in the denominator. – 2010-12-24
-
0...and how sure are you that the denominator polynomial has no zeroes within the interval of integration? – 2010-12-24
-
0I have plotted the function for $x=0$ (the denominator has no $x$ anyway). Both $f(v)$ & $g(v)$ decays for $v>2$. And the functions are smooth. – 2010-12-24
-
0I'm able integrate using simpson's rule. But the number of function evaluation is more. – 2010-12-24