Построить аппроксимирующий полином степени m=2 методом МНК

Program chetire; 
uses CRT;
var n,m:integer;
 
function f(x: real): real;
begin
  f := sqrt(abs(x));
end;
 
procedure MNK(n,m:integer);
Var x,x1,y: array [1..100] of real; a: array [1..100, 1..100] of real;
i,j,z,x0,k:integer; s,sum,sumf,r,p, g,pmax,max:real;
begin 
//Постройка полинома степени m
  write('x0='); //Начало интервала
  readln(x0);
  x[1]:=x0;
  write('s=');
  readln(s);
  for i:=2 to n do
    x[i]:=x[i-1]+s; //Заполнение массива Х узлов с равномерным шагом
  for i:=1 to n do
    x1[i]:=x[i];
  for i:=1 to n do
    y[i]:=f(x[i]); //Заполнение массива У
  for i:=1 to m+1 do //Заполнение полинома
    for j:=1 to m+1 do
      //a[строка,столбец]
      begin
      sum:=0; sumf:=0;
      for z:=1 to n do begin//Подсчёт суммы (элемента)
        if j=1 then a[i,m+2]:=a[i,m+2]+(power(x[z],(i-2+j))*y[z]); //Подстчёт столбца свободных членов
        sum:=sum+power(x[z],(i-2+j));
        end;
      a[i,j]:=sum;    
      end;
 
   
for i:=1 to n do begin //Вывод полученной матрицы
  write(x[i]:6:2,' ');
  end;
 
writeln;
 
 
for i:=1 to n do begin //Вывод полученной матрицы
  write(y[i]:6:2,' ');
  end;
 
writeln;  
  
writeln('Получена матрица');    
for i:=1 to m+1 do begin //Вывод полученной матрицы
  for j:=1 to m+2 do write(a[i,j]:6:2,' ');
  writeln;
  end;
 
//Решение матрицы методом Гаусса
for k:=1 to m+1 do begin
//Единицы
a[k,m+2]:=a[k,m+2]/a[k,k];
for j:=m+1 downto k do a[k,j]:=a[k,j]/a[k,k];
//Прямой ход
for j:=k+1 to m+1 do begin
  r:=a[j,k];
   for i:=k to m+1 do a[j,i]:=a[j,i]-r*a[k,i];
   a[j,m+2]:=a[j,m+2]-r*a[k,m+2];
  end;
end;
//Обратный ход
for k:=m+1 downto 1 do
begin
r:=0;
for j:=k+1 to m+1 do
begin
g:=a[k,j]*x[j];
r:=r+g;
end;
x[k]:=(a[k,m+2]-r)/a[k,k];
end;
write('f(x):');
for i:=1 to m+1 do
write('x[',i-1, ']*', x[i], ' + ');
writeln;
 
//Подсчёт погрешности
for i:=1 to n do begin
  p:=0;
  for j:=1 to m+1 do
    p:=p+power(x1[i],j-1)*x[j];
    write('p[',i, ']=', p, ' ');
    
  end;
writeln;  
writeln('Количество узлов сетки =',n,' Степень полинома =',m);
writeln;
end;
 
begin
n:=7; m:=2;
MNK(n,m);
end.