% polynom, polynom genom givna punkter
  clear, clf
  x=[0  4  6.5]';  y=[1.42  6.18  4.75]'; 
  A=[ones(size(x)) x x.^2], c=A\y
  X=0:0.1:8;  P2=c(1)+c(2)*X+c(3)*X.^2;
  stem(x,y), hold on, plot(X,P2,'g')
  title('Andragradspolynom genom tre punkter, heldraget'), pause(1)
% Lägg till en punkt:
  x=[x; 8];  y=[y; 5.55]; stem(x,y)    
  A=[ones(size(x)) x x.^2 x.^3], cc=A\y
  P3=cc(1)+cc(2)*X+cc(3)*X.^2+cc(4)*X.^3;
  plot(X,P3,'r--')
  xlabel('Tredjegradspolynom genom fyra punkter, streckat')
% figure(2), ribbon(X,P3,1)
