> restart;
> with(LinearAlgebra):
> Newtonkoeff:=proc(x,f)
> #x ist der Vektor der Stuetzszellen
> #f ist der Vektor der Stuetzwerte
> # option trace;
> local y,n,i,k;
> n:=Dimension(x);
> y:=f;
> for k from 1 to n-1 do
>     for i from n by -1 to k+1 do
>        y[i]:=evalf((y[i]-y[i-1])/(x[i]-x[i-k]));
>     end do;
> end do;
> evalm(y);
> end proc:  
> 
> x:=Vector([0,1,3]);
> y:=Vector([1,3,2]);
> w:=Newtonkoeff(x,y);