Вводятся натуральные числа X, Y, m и n. Вычислить сумму Xm+Yn (для определения степени числа написать процедуру)

//Алгоритм Евклида. Стоит отметить, что он также работает только с натуральными числами
//дробные и отрицательные числа не являются натуральными
//чаще и 0 не считается натуральным, но в этом решении считается
//сам алгоритм можно посмотреть в ВИКИ. Он простой 
function algorithEvklid(a, b : integer) : integer;
begin 
   while (b <> 0) and (a <> 0) do
   begin
      if a > b then a := a mod b
      else b := b mod a;
   end; 
   algorithEvklid := a + b;
end;
 
//Функция возведения в степень. a и pow должны быть натуральными числами
//дробные и отрицательные числа не являются натуральными
//чаще и 0 не считается натуральным, но в этом решении считается
function pow(a, pow: integer) : integer;
begin
   if (a = 0) and (pow = 0) then result := 0
   else result := 1;
   for var i := 1 to pow do result := result * a;
   pow := result;
end;
 
//Читаем число из консоли. Если оно удовлетворяет диапазону (startRange..endRange), то 
//возвращаем число. Если не удовлетворяет, то просим ввод заново.
function readNumber(text : String; startRange, endRange : integer) : integer;
begin
   repeat
      write(text, ' [', startRange, ' <= число <= ', endRange, ']: ');
      readln(result);
   until (result >= startRange) and (result <= endRange);
   readNumber := result;
end;
 
var
   a, b : integer;
   x, y : integer;
   m, n : integer;
begin
   writeln('Задание №1: Вычислить сумму X^m+Y^n');
   x := readNumber('Введите число Х', 0, 1000);
   y := readNumber('Введите число Y', 0, 1000);
   m := readNumber('Введите степень m', 0, 1000);
   n := readNumber('Введите степень n', 0, 1000);
   var sum := pow(x, m) + pow(y, n);
   writeln('Ответ: ', x, '^', m, ' + ', y, '^', n, ' = ', sum, ';');
   
   writeln('==============================================');
   writeln('Задание №2: Найти НОД 2-х чисел по алгоритму Евклида');
   a := readNumber('Введите число a', 0, 1000);
   b := readNumber('Введите число b', 0, 1000);;  
   writeln('Ответ: НОД(', a, ', ', b, ') = ', algorithEvklid(a, b), ';');
end.

//Алгоритм Евклида. Стоит отметить, что он также работает только с натуральными числами
//дробные и отрицательные числа не являются натуральными
//чаще и 0 не считается натуральным, но в этом решении считается
//сам алгоритм можно посмотреть в ВИКИ. Он простой 
procedure algorithEvklid(a, b : integer);
begin 
   write('НОД(', a, ',', b, ') = ');
   while (b <> 0) and (a <> 0) do
   begin
      if a > b then a := a mod b
      else b := b mod a;
   end; 
   writeln(a + b, ';');
end;
 
//Функция возведения в степень. a и pow должны быть натуральными числами
//дробные и отрицательные числа не являются натуральными
//чаще и 0 не считается натуральным, но в этом решении считается
procedure pow(var a : integer; pow : integer);
begin
   if pow >= 2 then for var i := 2 to pow do a := a * a
   else if (pow = 0) and (a = 0) then a := 0
   else if (pow = 0) and (a <> 0) then a := 1;
end;
 
//Читаем число из консоли. Если оно удовлетворяет диапазону (startRange..endRange), то 
//возвращаем число. Если не удовлетворяет, то просим ввод заново.
function readNumber(text : String; startRange, endRange : integer) : integer;
begin
   repeat
      write(text, ' [', startRange, ' <= число <= ', endRange, ']: ');
      readln(result);
   until (result >= startRange) and (result <= endRange);
   readNumber := result;
end;
 
var
   a, b : integer;
   x, y : integer;
   m, n : integer;
begin
   writeln('Задание №1: Вычислить сумму X^m+Y^n');
   x := readNumber('Введите число Х', 0, 1000);
   y := readNumber('Введите число Y', 0, 1000);
   m := readNumber('Введите степень m', 0, 1000);
   n := readNumber('Введите степень n', 0, 1000);
   
   write('Ответ: ', x, '^', m, ' + ', y, '^', n, ' = ');
   pow(x, m);
   pow(y, n);
   writeln(x + y, ';');
   
   writeln('==============================================');
   writeln('Задание №2: Найти НОД 2-х чисел по алгоритму Евклида');
   write('Введите число a: ');
   readln(a);   
   write('Введите число b: ');
   readln(b);  
   algorithEvklid(a, b);
end.