Метод Крамера надо доделать .

Program ll;
uses crt;
procedure Kramer(var k2,k3,k4,k6,k7,k8,k10,k11,k12,k14,k15,k16:real);
var x,y,z,D,Dx,Dy,Dz:real;
begin
D:=(k2*k7*k12)+(k4*k6*k11)+(k3*k8*k10)-(k4*k7*k10)-(k2*k8*k11)-(k3*k6*k12);
 if D<>0 then
 Dx:=(k14*k7*k11)+(k8*k10*k15)+(k6*k11*k16)-(k7*k10*k16)-(k15*k6*k12)-(k14*k8*k11);
 Dy:=(k2*k15*k12)+(k3*k16*k10)+(k4*k14*k11)-(K4*k15*k10)-(k3*k14*k12)-(k2*k16*k11);
 Dz:=(k2*k7*k16)+(k3*k8*k14)+(k4*k6*k15)-(k4*k7*k14)-(k8*k15*k2)-(k3*k6*k16);
 x:=Dx/D;
 y:=Dy/D;
 z:=Dz/D;
 writeln('x=',x:1:2);
 writeln('y=',y:1:2);
 writeln('z=',z:1:2);
 end;
else
showmessage('net reshenij');
Const
n=3;
m=4;
var a:array[1..n,1..m] of integer;
i,j:integer;
Begin
clrscr;
writeln('vvedute znacheniy dlay masuva');
For i:=1 to 3 do
For j:=1 to 4 do
Begin
readln(a[i,j]);
end;
for i:=1 to n do
begin
for j:=1 to m do
write(a[i,j]:4,'');
writeln('');
readln
end;
Kramer
end.

uses crt;
type
 Tmass=array[1..20] of real;
 Tmatrix=array[1..20] of Tmass;
procedure Per(k,n:integer;var a:Tmatrix;var p:integer);
var z:Real;j,i:integer;
begin
z:=abs(a[k,k]);
i:=k;
p:=0;
for j:=k+1 to n do
  begin
   if abs(a[j,k])>z then
    begin
     z:=abs(a[j,k]);
     i:=j;
     p:=p+1;{счетчик перестановок}
    end;
  end;
if i>k then
for j:=k to n do
  begin
   z:=a[i,j];
   a[i,j]:=a[k,j];
   a[k,j]:=z;
  end;
end;
function Znak(p:integer):integer;
begin
if p mod 2=0 then
Znak:=1 else Znak:=-1;
end;
procedure Opr(n:integer;a:tmatrix;var det:real);
var k,i,j,p:integer;r:real;
begin
det:=1.0;
for k:=1 to n do
  begin
   if a[k,k]=0 then Per(k,n,a,p);{перестановка строк}
   det:=znak(p)*det*a[k,k];{вычисление определителя}
   for j:=k+1 to n do {пересчет коэффициентов}
    begin
      r:=a[j,k]/a[k,k];
      for i:=k to n do
       begin
        a[j,i]:=a[j,i]-r*a[k,i];
       end;
    end;
  end;
end;
var a:Tmatrix;
    c:array[1..20] of Tmatrix;
    b,x:Tmass;
    det,det1:real;
    n,k,j,i:integer;
begin
clrscr;
write('Порядок системы n=');
readln(n);
writeln('Введите коэффициенты системы:');
for i:=1 to n do
for j:=1 to n do
read(a[i,j]);
readln;
writeln('Введите свободные члены:');
for i:=1 to n do
read(b[i]);
readln;
clrscr;
writeln('Расширенная матрица системы:');
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j]:7:2);
  write(b[i]:9:2);
  writeln;
 end;
Opr(n,a,det);{определитель системы}
for i:=1 to n do
 begin
  for k:=1 to n do
   begin
    for j:=1 to n do
    c[i][k,j]:=a[k,j];
    c[i][k,i]:=b[k];
   end;
  Opr(n,c[i],det1);
  if(det=0)and(det1=0) then
    begin
     writeln('Система не определена!');
     readln;
     exit;
    end;
  if(det=0)and(det1<>0) then
    begin
     writeln('Система не имеет решений!');
     readln;
     exit;
    end;
  x[i]:=det1/det;
 end;
writeln('Корни сиcтемы:');
for i:=1 to n do
writeln('x',i,'=',x[i]:7:3);
readln
end.

writeln('x=',x[1]:7:3);
writeln('y=',x[2]:7:3);
writeln('z=',x[3]:7:3);

Program ll;
uses crt;
procedure Kramer;
var x,y,z,D,Dx,Dy,Dz:real;
    k2,k3,k4,k6,k7,k8,k10,k11,k12,k14,k15,k16:real;
begin
D:=(k2*k7*k12)+(k4*k6*k11)+(k3*k8*k10)-(k4*k7*k10)-(k2*k8*k11)-(k3*k6*k12);
if D<>0 then
   begin
    Dx:=(k14*k7*k11)+(k8*k10*k15)+(k6*k11*k16)-(k7*k10*k16)-(k15*k6*k12)-(k14*k8*k11);
    Dy:=(k2*k15*k12)+(k3*k16*k10)+(k4*k14*k11)-(K4*k15*k10)-(k3*k14*k12)-(k2*k16*k11);
    Dz:=(k2*k7*k16)+(k3*k8*k14)+(k4*k6*k15)-(k4*k7*k14)-(k8*k15*k2)-(k3*k6*k16);
    x:=Dx/D;
    y:=Dy/D;
    z:=Dz/D;
    writeln('x=',x:1:2);
    writeln('y=',y:1:2);
    writeln('z=',z:1:2);
   end
else
writeln('net reshenij');
end;
Const
n=3;
m=4;
var a:array[1..n,1..m] of real;
i,j:integer;
Begin
clrscr;
writeln('vvedute znacheniy dlay masuva');
For i:=1 to 3 do
For j:=1 to 4 do
readln(a[i,j]);
for i:=1 to n do
 begin
  for j:=1 to m do
  write(a[i,j]:6:2);
  writeln;
 end;
Kramer;
readln
end.