// Calcul du PGCD de deux entiers naturels
clear
a=input("a=")
b=input("b=")
//reste devrait fonctionner a la place de modulo
r=modulo(a,b)
  while r<>0
    a=b;
    b=r;
    r=modulo(a,b);
  end
PGCD=b