Собираем Ошибки! Опции

[Altair]

В этой теме сообщайте обо всех ошибках (опечатках), найденных Вами в FAQ'е.

В этой теме премодерация отключена.
Ваше сообщение появиться не позднее чем через минуту после отправки!

[IUnknown]

Хм. По мотивам этой темы [Слау методом простых итерации. Не решаемое? - Lapp] выясняется, что в FAQ-е лежит не всегда корректно работающая программа решения систем методом Крамера. Нельзя в ней подобным образом находить дискриминант. Исправьте ту функцию Det. Скажем, вот на такую:

" (Показать/Скрыть)

Эта функция выдрана из одной моей старой программы, она была написана лет 10 назад, и жалоб на нее не было. А приведенная в FAQ-е функция на матрице из первого поста отрабатывает неверно, из-за чего результат также получается неправильным: <3, 0, 1> вместо <3, 1, 1>...

function det(var p; const size: integer): double;

   function minusOne (n : integer): integer;
   begin
      minusOne := (1 - 2 * Byte (Odd (n)));
   end;
   function get_addr(i, j : integer;
            const n: integer) : integer;
   begin
      get_addr := pred (i) * n + j
   end;

   type vector = array[1 .. n * n] of double;

var
   my_p : vector absolute p;
   pp : ^vector;
   s : double;
   i, j, curr: integer;
begin
   s := 0.0;
   if size = 2 then
   begin
      det := my_p[1]*my_p[4] - my_p[2]*my_p[3]; exit
   end;

   for i := 1 to size do
   begin
      GetMem (pp, Sqr (pred (size)) * sizeof (double));
      curr := 1;
      for j := 1 to size do
         if j <> i then
         begin
            move (my_p[get_addr (j, 2, size)],
                  pp^[get_addr (curr, 1, pred (size))],
                  pred(size) * sizeof(double));
            inc (curr);
         end;

      s := s + minusOne (succ(i)) * my_p[get_addr (i, 1, size)] *
               det(pp^, pred (size));
      FreeMem (pp, Sqr (pred (size)) * sizeof (double))
   end;
   det := s
end;

Сообщение отредактировано: Lapp - 12.05.2011 23:48

[Lapp]

Хм. ...

Спасибо, Володь. Я сделал (не проверял).
Не возражаешь, если я перенесу твой мессадж в тему про ошибки в FAQ?