% 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, f1=plot(X,P2,'g');

% 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;
  f2=plot(X,P3,'r--');
  
%Skriv ut info om plotten
  legend([f1,f2],'2-gradspolynom','3-gradspolynom')
  
% Prova 3D-ritning av tredjegradspolynomet
  figure(2), ribbon(X,P3,2), hold on
  plot(X,P3,'o')
  for i=10:10:80
	  plot3(X(i),P3(i),P3(i),'*')
  end
  xlabel('x-axeln')
  ylabel('y-axeln')
  zlabel('z-axeln')
  title('Tredjegradspolynomet  med ribbons')
