% bana5.m, femtegradspolynom genom alla punkterna
  x=[0.5 1.5 2.5 3 4 5]';  y=[3 1.5 1.5 1 1 0]'; 
  A=[ones(6,1) x  x.^2  x.^3  x.^4  x.^5],  c=A\y;  koeff=c'
  X=x(1):0.1:x(6);   
  P=c(1)+c(2)*X+c(3)*X.^2+c(4)*X.^3+c(5)*X.^4+c(6)*X.^5; 
  % Beräkningen av P kan också göras med Horners schema:
  % P=0;  for k=6:-1:1, P=P.*X+c(k); end  
  subplot(2,2,4), plot(x,y,'o', X,P),  axis equal

 
  