Date: Oct 5, 2012 12:57 AM
Author: Brad Cooper
Subject: MuPAD verion 4 Problem with Newton Poly Interpolation

I set up a polynomial in T to perform Newton Interpolation

P := (N,T) -> sum(a[M]*product((T-t[k-1]),k=1..M),M=0..N);

Check is correct for 4 terms

P(3,T)

Get the table of data

t:=table(0=-0.9749279122,1=-0.7818314825,2=-0.4338837391,3=0,4=0.4338837391,5=0.7818314825,6=0.9749279122):

Check the polynomial for 4 terms now

P(3,T)

It is a mess. MuPAD substitutes the whole table each time instead of entries from the table.

I am using a table because the index can start at 0 and the formula then matches the one in the text.

Is there a workaround?

Cheers,
Brad