% bana1,  rutschbana i två stycken
  clear, clf
  xx=[0.5 1.5 2.5 3 4 5]';  yy=[3 1.5 1.5 1 1 0]'; 
  subplot(2,2,1), plot(xx,yy,'o', xx,yy),  axis equal
% Kvadratisk interpolation genom punkterna 1 - 3
  x=xx(1:3);  y=yy(1:3);   A=[ones(3,1)  x  x.^2],  c=A\y  
  X2=0.5:0.1:2.5;  P2=c(1)+c(2)*X2+c(3)*X2.^2; 
 
% Tredjegradspolynom genom punkterna 3 - 6
  x=xx(3:6);  y=yy(3:6);   A=[ones(4,1)  x  x.^2  x.^3],  c=A\y  
  X3=2.5:0.1:5;  P3=c(1)+c(2)*X3+c(3)*X3.^2+c(4)*X3.^3; 
  subplot(2,2,2), plot(xx,yy,'o', X2,P2, X3,P3),  axis equal
 
  
