Метод прогонки на Паскале

uses Crt;
type
 vector=array[1..100]of real;
 matrix=array[1..100,1..100]of integer;
const
 e:real=0.001; w=4;
var
  i,j,k,max,n,g,c:integer;
  ek,sum1,sum2:real;
  a,alpha,r:matrix;
  b,xc,xprev,tmp,x:vector;
 
  y,u,o,m,t,p,q:array[1..w] of real;
 
 
 begin
 
  writeln('vvedite poriadok matrici');
  readln(n);
  g:=n*(n-2);
 
  writeln('vvedite pervie ',n,' potencialov');
  for i:=1 to n do   {Zadaem pervie n potencialov}
   begin
    readln(b[i]);
    b[i]:=(-1)*b[i];
   end;
 
  writeln('vvedite poslednie ',n,' potencialov');
  for i:=(g-n+1) to g do   {Zadaem poslednie n potencialov}
   begin
    readln(b[i]);
    b[i]:=(-1)*b[i];
   end;
 
  for j:=0 to (n-3) do
   begin
    for i:=1 to n do
     begin
      c:=0;
      if i=1 then {Dlya granichnih tochek (sleva)}
       begin
         a[i+n*j,i+n*j]:=-3;
         a[i+n*j,i+n*j+1]:=1;
         a[i+n*j,i+n*j+n]:=1;
         if j<>0 then a[i+n*j,i+n*j-n]:=1;
       end;
      if i<>1 then c:=c+1;
 
      if i=n then   {Dlya granichnih tochek (sprava)}
       begin
        a[i+n*j,i+n*j]:=-3;
        a[i+n*j,i+n*j-1]:=1;
        if j<>0 then
         a[i+n*j,i+n*j-n]:=1;
        if j<>(n-1) then a[i+n*j,i+n*j+n]:=1;
       end;
      if i<>n then c:=c+1;
 
      if c=2 then           {Dlya ostalnih tochek}
       begin
        a[i+n*j,i+n*j]:=-4;
        a[i+n*j,i+n*j+1]:=1;
        a[i+n*j,i+n*j-1]:=1;
        if j<>0 then a[i+n*j,i+n*j-n]:=1;
        if j<>(n-1) then a[i+n*j,i+n*j+n]:=1;
       end;
     end;
   end;
 
   writeln('ishodnaya matrica');
   for i:=1 to g do
    begin
     for j:=1 to g do
     write(' ',a[i,j]:1);
     writeln;
    end;