Найти НОД(а, b). (программа Паскаль)

program z1;
uses crt;
var a, b: integer;
begin
read (a, b);
while a<>b do
if a>b then a:=a-b else b:=b-a;
 writeln (a);
end.