Расчет матрицы и выведение графика

writeln('Введите коэффициенты матрицы нажимая ENTER:');
   for k:=1 to n do 
    for j:=1 to n do 
     begin
      write ('a[',k,',',j,']='); 
      read(a[k,j]); 
     end;
   writeln('Введите свободные члены нажимая ENTER:');
    for k:=1 to n do 
     begin
      write ('b[',k,']=');
      read (b[k]);
     end;
  writeln; 
  
 
  writeln('Матрица:'); 
  for i:=1 to n do 
   begin
    for j:=1 to n do
        write (a[i, j],' '); 
        write (b[i],' '); 
        writeln;
   end;
  writeln;
 
    for i:=1 to n-1 do 
     begin
      for k:=i+1 to n do 
       begin
if a[i,i]=0 then a[i,i]:=c; 
 
        delenie:= -1*(a[k,i] / a[i,i]); 
         for j:=1 to n do 
             a[k,j]:= a[k,j] + a[i,j] * delenie;
             b[k]:= b[k] + b[i] * delenie; 
       end;
    end;
    writeln('Результат прямого хода:'); 
    for i:=1 to n do
      begin
       for j:=1 to n do
        write (a[i, j],' '); 
        write (b[i],' '); 
        writeln;
      end;
    writeln;
    for i:=n downto 1 do 
     begin
      for j:= i+1 to n do
          b[i]:= b[i] - a[i,j] * x[j]; 
          x[i]:= b[i] / a[i,i]; 
     end;
 
     writeln('Ответ:'); 
    for i:=1 to n do
     writeln('X[',i:1,']=', x[i], '.'); 
   readln;
end.