a=2;b=3;n=7;
h=(b-a)/(n-1);x=[];y=[];
//
for i=1:n
  x(i)=a+(i-1)*h;
  y(i)=x(i)^4;
end
for i=1:n-1
  for j=n:-1:i+1
    y(j)=(y(j)-y(j-1))/(i*h);
  end
end
//
t=2.5;p=y(n);
for i=n-1:-1:1
  p=p*(t-x(i))+y(i);
end
