function testPolyfit()

  % Parameter für eine Parabel
  a = 1;
  b = 2; 
  c = -3;
  p1 = [a,b,c]
  n = 15;
  % Parabel: x- und y-Werte
  x = linspace( -2, 1, n );
  y = a*x.^2 + b*x + c;
  % Verrauschung um Werte zwischen [-0.5, 0,5] 
  rd = rand( 1, n ) - 0.5; 
  y = y + rd;
  % Interpolation mit Parable
  p = polyfit( x, y, 2 )
  yp = polyval( p, x );
  % Zeichnen
  hold on 
  plot( x, y,  '*r', x, yp, '-b' );
  legend( 'Original', 'Interpolation' )
  hold off

end
