Ошибка в Паскале

[Mystery08]

Помогите!!! Что здесь не так???

Program Metod_Gauss;
Const n=12;
A:array[1..n,1..n] of real=((1,0,0,0,0,0,0,0,0,0,0,0),(0,0,0,0,1,6,36,216,0,0,0,0),(0,0,0,0,0,0,0,0,1,14,196,2744),
(1,6,36,216,0,0,0,0,0,0,0,0),(0,0,0,0,1,14,196,2744,0,0,0,0),(0,0,0,0,0,0,0,0,1,20,400,8000),
(0,-1,-12,108,0,1,12,108,0,0,0,0),(0,0,0,0,0,-1,-28,-588,0,1,28,588),(0,0,-2,-36,0,0,2,36,0,0,0,0),
(0,0,0,0,0,0,-2,-84,0,0,2,84),(0,0,2,0,0,0,0,0,0,0,0,0),(0,0,0,0,0,0,0,0,0,0,2,120));
B:array[1..n] of real=(0.033,-0.118,-1.904,-0.118,-1.904,-5.616,0,0,0,0,0,0);
type matrtype=array[1..n,1..n+1] of real;
Var
i,j:byte;
matr:matrtype;
Procedure Gausse(var matr:matrtype; n:byte);
var i,j,k:byte;
begin
  for i:=1 to n-1 do
  for j:=i+1 to n do
  for k:=n+1 downto i do
  matr[j,k]:=matr[j,k]-matr[i,k]/matr[i,i]*matr[j,i];
  for i:=n downto 1 do
  begin
   for j:=i+1 to n do
   Matr[i,n+1]:=Matr[i,n+1]-Matr[i,j]*Matr[j,n+1];
   Matr[i,n+1]:=Matr[i,n+1]/Matr[i,i];
   end;
   end;
   begin
    Writeln('Решение системы из n линейных уравнений');
    for i:=1 to n do
    begin
    writeln('Вводим уравнение i');
    for j:=1 to n do
    begin
    write('A[',i,',',j,']=');
   end;
  Write('B[',i,']=');
  end;
  writeln('Cистема линейных уравнений');
  for i:=1 to n do
  begin
  for j:=1 to n do
  write(A[i,j]:5:3);
  writeln(B[i]:5:3);
  end;
  for i:=1 to n do
  begin
  for j:=1 to n do
  matr[i,j]:=A[i,j];
  matr[i,n+1]:=B[i];
  end;
  Gausse(matr,n);
  writeln('Решение системы уравнений');
  for i:=1 to n do
  writeln('C',i,'=',matr[i,n+1]:5:3);
  end.

выводит ошибку: вещественное деление на 0 (метод Гаусса.pas, строка 18)

[Mad_Cat]

дык так и есть)

matr[j,k]:=matr[j,k]-matr[i,k]/matr[i,i]*matr[j,i];

например элемент A[2][2]

[Mystery08]

А как мне исправить программу, чтоб выдавало без ошибки?)

[Izobara]

О Господи! Это Вы писали? Мда, это не исправить... Особенно константы - гениальное решение. Только вручную, только хардкор!

 {Метод Гаусса}
uses
  crt;

const
  n = 3;

var
  i, j, k: integer;
  a: array
  [1..n, 1..n] of
  real;
  b: array[1..n] of
  real;
  x: array[1..n] of
  real;
  h: real;

begin
  clrscr;
  for i := 1 to n do
  begin
    writeln
    ('Уравнение:', i);
    for j := 1 to n do
    begin
      write('Введите A(', i, ',', j, '):');
      readln(a[i, j]);
    end;
    write('Свободный член: B(', i, '):');
    readln(b[i]);
  end;
  for i := 1 to n do
  begin
    for j := i + 1 to n
    do
    begin
      a[j, i] := -a[j, i] / a
      [i, i];
      for k := i + 1 to n
        do 
        a[j, k] := a
        [j, k] + a[j, i] * a
        [i, k];
      b[j] := b[j] + a
      [j, i] * b[i];
    end;
  end;
  x[n] := b[n] / a
  [n, n];
  for i := n - 1
  downto 1 do
  begin
    h := b[i];
    for j := I + 1 to n
      do
      h := h - x[j] * a[i, j];
    x[i] := h / a[i, i];
  end;
  writeln('Решение системы:');
  for i := 1 to n do
    writeln('X (', i, ')=', x[i]:3:5);
  readln;
end.