Паскаль.Исправить ошибку!

[kikimka]

Проблема такая, не правильное считает матрицу C по формуле.С=В*В+Е
E-единичная матрица(a[i,j])
В-простая матрица(b[i,j])

uses crt;
cоnst n=5;
var
a,b,c:array [1..n,1..n] of integer;
i,j,k,s:integer;
begin
clrscr;
for i:=1 to n do
begin
for j:=1 to n do
begin
if i=j then a[i,j]:=1
else
a[i,j]:=0;
write(a[i,j]:3);
end;
readln;
end;
writeln;
begin
for i:=1 to n do
for j:=1 to n do
b[i,j]:=random(25);
end;
for i:=1 to n do
begin
for j:=1 to n do
write(b[i,j]:3);
readln;
end;
begin
for i:=1 to n do
for j:=1 to n do
begin
s:=0;
for k:=1 to n do
S:=S+b[i,k]*b[k,j]+a[i,j];
c[i,j]:=s;
end;
for i:=1 to n do
begin
for j:=1 to n do
write(c[i,j]:5);
writeln;
end;
readkey;
end;
end.

Матрицу В*В считает правильно а при прибавлении единичной матрицы считает не правильно, подскажите как исправить.Спасибо большое.
З.Ы.Искала,но не нашла

[Грымзик]

Зачем прибавление единичной матрицы в этот цикл загонять?
for k:=1 to n do
S:=S+b[i,k]*b[k,j]+a[i,j];
Надо

for i:=1 to n do
     for j:=1 to n do
     begin
          c[i,j]:=a[i,j];
          for k:=1 to n do
               c[i,j]:=c[i,j]+b[i,k]*b[k,j];
     end;

[kikimka]

Ура! Спасибо огромное, заработало!Терь все ясно!