решение задачи Дирихле методом Гаусса-Зейделя на Паскале. Объяснить код программы

[AleSanchez]

Ребят, помогите пожалуйста, завтра нужно преподавателю объяснить данную программу:
Это решение задачи Дирихле методом Гаусса-Зейделя, сам принцип я понимаю, но вот программный код вообще не шарю

Program lab2;
const n=20;
type matrix=array [0..n,0..n] of real;
var u,u0 :matrix;
i,j :integer;
s,eps,h,x,y :real;
begin
h:=1/n;
write('eps='); readln(eps);
for i:=1 to n-1 do
for j:=1 to n-1 do u0[i,j]:=0;
for j:=0 to n do
begin
y:=j*h;
u0[0,j]:=exp(y);
u0[n,j]:=exp(1+y);
end;

for i:=0 to n do
begin
x:=i*h;
u0[i,0]:=exp(x);
u0[i,n]:=exp(x+1);
end;

repeat
for i:=0 to n do
for j:=0 to n do
u[i,j]:=u0[i,j];


for i:=1 to n-1 do
begin
x:=i*h;
for j:=1 to n-1 do
begin
y:=j*h;
u[i,j]:=1/4*(u[i+1,j]+u[i-1,j]+u[i,j+1]+u[i,j-1]-2*exp(x+y)*sqr(h));
end;
end;

s:=0;
for i:=1 to n do
for j:=1 to n do
if s<abs(u[i,j]-u0[i,j]) then s:=abs(u[i,j]-u0[i,j]);
for i:=0 to n do
for j:=0 to n do u0[i,j]:=u[i,j];
until s<eps;

s:=0;
for i:=0 to n do
begin
x:=i*h;
for j:=0 to n do
begin
y:=j*h; if s<abs(u[i,j]-exp(x+y)) then s:=abs(u[i,j]-exp(x+y));
end;
end;
writeln ('s=',s);
readln
end.

Заранее спасибо!

________
1) Название темы по правилам форума должно адекватно отражать суть решаемой задачи/проблемы.
На первый раз я исправил.
В последующем, темы с подобным названием будут закрываться или удаляться,
а автор такой темы будет получать штрафные баллы.

2) Код нужно оформлять по правилам:
тегом [CODE]..[/СODE]
(это кнопочка на панели форматирования с решёточкой #)

Не забывайте об этом!

Модератор.