Обратная матрица

[ivan.tiran]

Пишу прогу по созданию обратной матрицы. За основу взял матрицу третьего порядка (потом переделаю под цикл). Появилась проблема, почему-то в конце b[i,j] выпадает таким же, как я задал в начальном цикле. В чём проблема?

		program prog;
		uses crt;
		var
		a:array[1..10,1..10] of real;
    b:array[1..10,1..10] of real;
		i,j,n:integer;
		BEGIN
  write('Vvedite razmernost kvadratnoy matrici:');
		read(n);
		        for i:=1 to n do
		        for j:=1 to n do
		        if i=j then
		        b[i,j]:=1
		        else
		        b[i,j]:=0;

		
		
            for i:=1 to n do
		        for j:=1 to n do
		        begin
		        writeln('Vvedite a[',i,',',j,']');
		        read(a[i,j]);
            end;
    

		a[1,1]:=a[1,1]/a[1,1];
		a[1,2]:=a[1,2]/a[1,1];
		a[1,3]:=a[1,3]/a[1,1];
		b[1,1]:=b[1,1]/a[1,1];
		b[1,2]:=b[1,2]/a[1,1];
		b[1,3]:=b[1,3]/a[1,1];


		a[2,1]:=a[2,1]-a[1,1]*a[2,1];
		a[2,2]:=a[2,2]-a[1,2]*a[2,1];
		a[2,3]:=a[2,3]-a[1,3]*a[2,1];
		b[2,1]:=b[2,1]-b[1,1]*a[2,1];
		b[2,2]:=b[2,2]-b[1,2]*a[2,1];
		b[2,3]:=b[2,3]-b[1,3]*a[2,1];

		a[3,1]:=a[3,1]-a[1,1]*a[3,1];
		a[3,2]:=a[3,2]-a[1,2]*a[3,1];
		a[3,3]:=a[3,3]-a[1,3]*a[3,1];
		b[3,1]:=b[3,1]-b[1,1]*a[3,1];
		b[3,2]:=b[3,2]-b[1,2]*a[3,1];
		b[3,3]:=b[3,3]-b[1,3]*a[3,1];

		a[2,1]:=a[2,1]/a[2,2];
		a[2,2]:=a[2,2]/a[2,2];
		a[2,3]:=a[2,3]/a[2,2];
		b[2,1]:=b[2,1]/a[2,2];
		b[2,2]:=b[2,2]/a[2,2];
		b[2,3]:=b[2,3]/a[2,2];

		a[3,1]:=a[3,1]-a[2,1]*a[3,2];
		a[3,2]:=a[3,2]-a[2,2]*a[3,2];
		a[3,3]:=a[3,3]-a[2,3]*a[3,2];
		b[3,1]:=b[3,1]-b[2,1]*a[3,2];
		b[3,2]:=b[3,2]-b[2,2]*a[3,2];
		b[3,3]:=b[3,3]-b[2,3]*a[3,2];

		a[3,1]:=a[3,1]/a[3,3];
		a[3,2]:=a[3,2]/a[3,3];
		a[3,3]:=a[3,3]/a[3,3];
		b[3,1]:=b[3,1]/a[3,3];
		b[3,2]:=b[3,2]/a[3,3];
		b[3,3]:=b[3,3]/a[3,3];

		a[2,1]:=a[2,1]-a[3,1]*a[2,3];
		a[2,2]:=a[2,2]-a[3,2]*a[2,3];
		a[2,3]:=a[2,3]-a[3,3]*a[2,3];
		b[2,1]:=b[2,1]-b[3,1]*a[2,3];
		b[2,2]:=b[2,2]-b[3,2]*a[2,3];
		b[2,3]:=b[2,3]-b[3,3]*a[2,3];

		a[1,1]:=a[1,1]-a[3,1]*a[1,3];
		a[1,2]:=a[1,2]-a[3,2]*a[1,3];
		a[1,3]:=a[1,3]-a[3,3]*a[1,3];
		b[1,1]:=b[1,1]-b[3,1]*a[1,3];
		b[1,2]:=b[1,2]-b[3,2]*a[1,3];
		b[1,3]:=b[1,3]-b[3,3]*a[1,3];

		a[1,1]:=a[1,1]-a[2,1]*a[1,2];
		a[1,2]:=a[1,2]-a[2,2]*a[1,2];
		a[1,3]:=a[1,3]-a[2,3]*a[1,2];
		b[1,1]:=b[1,1]-b[2,1]*a[1,2];
		b[1,2]:=b[1,2]-b[2,2]*a[1,2];
		b[1,3]:=b[1,3]-b[2,3]*a[1,2];

                             for i:=1 to n do
		                         for j:=1 to n do
                             writeln('b[',i,',',j,']=',b[i,j]);

  END.

Наверняка какая-нибудь ерунда, которую я в очередной раз не заметил)))

[ivan.tiran]

Кажется, уже разобрался сам, я просто изменяю a, а потом только делю b на а.