Dear all!

Rational function interpolation is highly recommended for my purpose in the literature.

However, I encountered the following problem: some poles, i.e. roots of the denominator polynomial, fall into the interval being interpolated, whereas the interpolated function itself is continuous in the interval.

Need your advice, how to treat these poles?

I've found absolutely nothing anywhere illuminating this issue... :confused:

I have the same problem. Moreover, I've found that if I compute the raional function minimizing the L2 norm (sum of squares of differences) using, for example, Newthon's method, the result will be continuous (without poles) and closer to the original function than one of NR's algorithm. I used 10 iterations!

It's very hard to understand.