You may wish to review the section "Advanced Numerical Integration" in Mathematica. One method is "DoubleExponentialOscillatory".

If you post the numeric values of all the parameters, and a range for x, I'll try to numerically integrate it and create an interpolation function f(x) for the chosen range.