On Tuesday 06 April 2010, Richard Fateman wrote:
> I suspect that replacing x^k by pow(x,k) in one place, will fix the
> interpolation program,
> where pow(a,b):= if (a=0)and(b=0) then 1 else a^b.
It looks like it does fix the problem with the power.
I have modified interpol.mac, but for my purpose it's terribly slow.
I killed my program after 20 minutes,
because doing it by my own I need only about 2 minutes
I usually vary the degree parameter for ratinterpol.
and select the nicest factored formula, which works for me.
Andre