% bana5n.m,  femtegradspolynom med Newtons ansats
  x=[0.5 1.5 2.5 3 4 5]';  y=[3 1.5 1.5 1 1 0]'; 
  ak=ones(6,1);  A=ak;         % Bygg upp A kolumn för kolumn  
  for kol=2:6, ak=ak.*(x-x(kol-1)); A=[A ak]; end
  c=A\y        
  X=x(1):0.1:x(6);             % Beräkna P(X) med Horners schema 
  P=0;  for k=6:-1:1, P=P.*(X-x(k))+c(k); end 
  plot(x,y,'o', X,P), axis equal 


