var x,y,weight       : RealArrayType_TL;
    i,j,n,poly_order : integer;
    std_dev,a,b,c    : real;
    coeff,coeff_err  : CoeffType_TL;

begin
  clrscr;
  writeln('Durchfuehrung einer Polynomregression');
  writeln;

  poly_order:=3; { Ordnung des Polynoms }

  a:=1;          { Die Koeffizienten    }
  b:=2;
  c:=3;

  n:=20;         { Anzahl der Wertepaare }

  { ------ Erzeugung der Daten ------ }

  for i:=1 to n do
  begin
    x[i]:=i;                                { x                              }
    y[i]:=a+b*i+c*sqr(i)+random/10;         { f(x) = a + bx + c x*x + random }
    weight[i]:=1.0;                         { Einheitswichtung               }
  end;

  polynom_regression(x,y,weight,n,poly_order,coeff,coeff_err,std_dev);

  writeln;
  writeln('Die Standardabweichung des Fits : ',std_dev);
  writeln;
  writeln('Die Liste der Koeffizienten');
  writeln;
  writeln('  Index     Koeff.      Fehler');
  writeln;
  for i:=0 to poly_order do
  begin
    writeln(i:4,'   ',coeff[i]:10:7,'    ',coeff_err[i]:10:8);
  end;
end.