Нахождение ортагональных векторов в пространстве, Среди n(n <= 10) векторов в пространстве (x;y;z) заданных своими к Опции

[latysh]

Среди n(n <= 10) векторов в пространстве (x;y;z) заданных своими координатами, имеется только 2 ортоганальных вектора.Написать программу нахождения ортоганальных векторов .Подготовить и использовать процедуру проверки ортоганальных векторов.
Вот чиво получилось,прога работает, о есть большие сомнения в формуле,если чиво неправильно,подскажите что исправить

program vektors;
var
x: array[1..10] of real;
y:array[1..10] of real;
z:array[1..10] of real;
k,j,i:integer;


BEGIN
writeln('Введите значение x');
for j:=1 to 10 do
read(x[j]);
writeln('Введите значение  y');
for j:=1 to 10 do
read(y[j]);
writeln('Введите значение z');
for j:=1 to 10 do
read(z[j]);
for k:=2 to 10 do
for i:=1 to 10 do
for j:=k to 10 do
if (x[j]*x[i]+y[j]*y[i]+z[j]*z[i])/(sqrt(x[j]*x[j]+y[j]*y[j]+z[j]*z[j])*sqrt(x[i]*x[i]+y[i]*y[i]+z[i]*z[i]))=0 then
writeln(x[j]:5:2,y[j]:5:2,z[j]:5:2 );
write('              ',x[i]:5:2,y[i]:5:2,z[i]:5:2);
end.

[Lapp]

имеется только 2 ортоганальных вектора.Написать программу нахождения ортоганальных векторов .Подготовить и использовать процедуру проверки ортоганальных векторов.
Вот чиво получилось,прога работает, о есть большие сомнения в формуле,если чиво неправильно,подскажите что исправить

latysh, я не совсем понял,

а. зачем тебе тройной цикл в проверке?
б. зачем делить на что-то, если ты проверяешь на равенство нулю?

И еще: избався от явного указания числа 10. Введи константу:

MaxN = 10;

Массивы заводи по максимуму:

x,y,z: array[1..MaxN]of double

, а в цикле проходись только до N (N вводится пользователем).
И на замечание Ozzи нужно обратить внимание.
Короче, я бы сделал как-то так:

function Ortog(i,j: integer): boolean;
const
  e=1e-12;
begin
  Ortog:=Abs(x[j]*x[i]+y[j]*y[i]+z[j]*z[i])<=e
end;

....
  i:=1;
  Found:=false
  while (i<=N)and not Found do begin
    j:=1+i;
    while j<=N do begin
      Found:=Ortog(i,j);
      Inc(j)
    end;
    Inc(i)
  end;
  if Found then WriteLn('Vectors ',i,' and ',j,' are orthogonal') else WriteLn('No orthogonal vectors were found')

Сообщение отредактировано: Lapp - 20.05.2011 0:19