I would like to evaluate the following integral:
$$I(x_0)=\int_{-1}^1 \left(\sum_{k=0}^n a_k \, T_k(x)\right) \, \mathrm{e}^{b(x-x_0)}\,\mathrm{d}x$$
with $a_k$ some constant coefficients comming from an Chebyshev interpolation of a function $f(x)$, $T_k(x)$ the Chebyshev polynomials, $b\in\mathbb{C}$, $n\in\mathbb{N}$ and $x_0>1$.
I assume, I can rewrite the integral as
$$I(x_0)=\mathrm{e}^{-b\,x_0}\sum_{k=0}^n a_k \left(\int_{-1}^1 T_k(x)\, \mathrm{e}^{b\,x}\,\mathrm{d}x\right)$$
In this question this integral is solved by a recurrence relation. I used this recurrence relation to evaluate the sum, but sadly the recurrence formular is not stable for large $n$. I also tryed to adopt the Clenshaw algorithm (Link) to overcome this problem, but still the evaluation diverges.
Is there a possibility to write this integral in explicit form?
Edit: I am wondering if replacing
$$\sum_{k=0}^n a_k \, T_k(x)$$ by the barycentric interpolation formula $$\sum_{j=0}^N a_j \, T_j(x) = \frac{\displaystyle \sum_{j=0}^N \frac{w_j}{x-x_j}f_j}{\displaystyle \sum_{j=0}^N \frac{w_j}{x-x_j}}, $$ where $$ w_j = \left\{ \begin{array}{cc} (-1)^j/2, & j=0\text{ or }j=N,\\ (-1)^j, & \text{otherwise} \end{array} \right. $$ is any better!