procedure normal_matrix(var feld:matrix;xwerte,wicht:werte;anzahl,polygrad:integer);

{ Berechnung der Normal-Matrix der Polynom-Regression  }

var i,j,k:integer;
    sum:real;

begin
  mat_set(feld,polygrad+1,0);
  for i:=1 to polygrad+1 do
     begin
       for j:=i to polygrad+1 do
          begin
            sum:=0.0;
            for k:=1 to anzahl do
               sum:=sum+wicht[k]*potenz(xwerte[k],i+j-2);
            feld[i,j]:=sum;feld[j,i]:=sum;
          end;
     end;
end;

