Расширенный алгоритм Евклида

var
m,n,d,x,y:integer;
begin
readln(m,n);
while (m<>0)and(n(<>0) do
if m>n then m:=m mod n else n:=n mod m;
d:=(m+n);
writeln(d);
end.

while (m<>0)and(n<>0) do

var
  a,b,d,x,y:Longint;
{ðàñøèðåííûé àëãîðèòì Åâêëèäà}
procedure Eq(a,b:longint; var d,x,y:longint);
var
  x1,y1,x2,y2,q,r:Longint;
begin
  if b=0 then
  begin
    d:=a;
    x:=1;
    y:=0
  end
  else
  begin
    x1:=0;
    x2:=1;
    y1:=1;
    y2:=0;
    while b>0 do
    begin
      q:=a div b;
      r:=a-q*b;
      x:=x2-q*x1;
      y:=y2-q*y1;
      a:=b;
      b:=r;
      x2:=x1;
      x1:=x;
      y2:=y1;
      y1:=y
    end;
    d:=a;
    x:=x2;
    y:=y2
  end;
end;
{îñíîâíàÿ ïðîãðàììà}
begin
  readln(a,b);
  Eq(a,b,d,x,y);
  writeln('Ðàñøèðåííûé àëãîðèòì Åâêëèäà(d(ÍÎÄ(a,b))=x*a+y*b):');
  writeln('d=',d,' x=',x,' y=',y)
end.

var
  a,b,d,x,y,x1,y1,x2,y2,q,r:Longint;
begin
  readln(a,b);
  if b=0 then
  begin
    d:=a;
    x:=1;
    y:=0
  end
  else
  begin
    x1:=0;
    x2:=1;
    y1:=1;
    y2:=0;
    while b>0 do
    begin
      q:=a div b;
      r:=a-q*b;
      x:=x2-q*x1;
      y:=y2-q*y1;
      a:=b;
      b:=r;
      x2:=x1;
      x1:=x;
      y2:=y1;
      y1:=y
    end;
    d:=a;
    x:=x2;
    y:=y2
  end;
  writeln('Ðàñøèðåííûé àëãîðèòì Åâêëèäà(d(ÍÎÄ(a,b))=x*a+y*b):');
  writeln('d=',d,' x=',x,' y=',y)
end.

function ExtEvklid(e,phi: integer; var x: integer; var y: integer): integer;
var 
  x1,x2,y1,y2,q,r: integer;
begin
  Result:= 0;
  x1:=0;
  x2:=1;
  y1:=1;
  y2:=0;
  if phi = 0 then 
  begin
    Result:=e;
    x:=1;
    y:=0;
  end
  else
  begin
    while phi > 0 do    
    begin
      q:=e div phi;
      r:=e-q*phi;
      x:=x2-q*x1;
      y:=y2-q*y1;
      e:=phi;
      phi:=r;
      x2:=x1;
      x1:=x;
      y2:=y1;
      y1:=y;
    end;
    Result:=e; // по условиям псевдокода 
    x:=x2;
    y:=y2;
  end;
end;

procedure ExtEvklid(e,phi: integer; var x: integer; var y: integer);
var 
  x1,x2,y1,y2,q,r: integer;
  xx1,yy1: integer;
begin
  x1 := 1; x2 := 0;
  y1 := 0; y2 := 1; // matrix
 
  q := e div phi;
  r := e mod phi;
 
  while(r <> 0) do 
    begin
      xx1 := x1;
      yy1 := y1;
 
      x1 := x2;
      x2 := xx1 - q*x2;
      y1 := y2;
      y2 := yy1 - q*y2;
 
      e := phi; 
      phi := r; 
 
      q := e div phi;
      r := e mod phi;
    end;
 
  x := x2;
  y := y2;
end;

ExtEvklid(e, phi, x, y);

d:= phi * (-x);

if x > 0 then d:= x else 
begin
...
end;