Ортонормированная матрица. Pascal

[Punk_tir]

Нужно написать прогу на Паскале, вычисляющую, является ли данная квадратная матрица ортонормированной, т.е. такой, что скалярное произведение (сумма всех произведений n-нового элемента одной строки на n-ный элемент другой строки) любой строки на саму себя = 1, а на любую другую строку = 0. Желательно, чтобы прога продолжала вычислять до того момента, пока не найдет несоответствующую строку, а не проверяла абсолютно все сочетания. В принципе алгоритм я придумал, но он не работает (это только основной кусок):

begin
i:=1; flag:=true; skal:=0;
repeat
j:=1;
repeat
for k:=1 to n do
skal:=skal+matrix[i,k]*matrix[j,k];
if (i=j) and (skal<>1) or (i<>j) and (skal<>0) then flag:=false;
j:=j+1
until not(flag) or (j=n+1);
i:=i+1
until not (flag) or (i=n+1);

if flag=true then

begin

writeln ('Ортонормированна');

end

else

begin
writeln ('Не ортонормированна');

end
end

[n@sok]

Переменную skal надо обнулять перед каждым циклом, т.е.

begin
i:=1; flag:=true;
repeat
j:=1;
repeat
skal:=0;
for k:=1 to n do
skal:=skal+matrix[i,k]*matrix[j,k];

[Punk_tir]

Буду невероятно признателен за помощь!!!

[Punk_tir]

Спасибо за совет! Только опять не работает должным образом... Подскажите, пожалуйста, может, я еще что-нибудь пропустил?

[n@sok]

Возможно это обусловлено неточным представлением дробных чисел. Поэтому нужно заменить условие (scal<>0) на (abs(scal)>1e-12)

[Punk_tir]

Боюсь, что дело не в этом. Я ввожу матрицу
1 0
0 1

а прога выдает что она не ортонормированна =( Может, все-таки что-то не так с алгоритмом?

[n@sok]

Ошибка где-то в другом месте. Алгоритм работает

[Punk_tir]

Спасибо за помощь, помогло!!!

[асх]

матрица является ортонормированной если произведение матрицы на транспорированную матрицу равно единичной матрице.
т.е.

код на С:
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>


void main()
{
int i,j,k=0;
int o[5][5];
randomize();
//заполнение матрицы
for (i=0; i<5; i++)
{
for (j=0; j<5; j++)
{
o[i][j]=random(5);
printf("%d ", *(*(o+i)+j));
}
printf("\n");
}
printf("\n\n");

//вычисление ортонормированная ли матрица
for (i=0; i<5; i++)
{
for (j=0; j<5; j++)
{
if (o[i][j]*o[j][i]!=1)
{
printf("Matrica ne ortonormirovanaya");
getch();
return;
}


}
printf("\n");
}
printf("Matrica ortonormirovanaya");
getch();
}