Función MCD para Octave

[alfonso]

Necesitaba esta función en Octave y no estaba por defecto, así que buscando la he encontrado en un powerpoint de una universidad italiana.
Como tampoco tengo claro en dónde va en el foro, la pongo aquí:

Código:

function [div] = mcd(x,y)

if(x == y)
   div = x;
elseif(x < y)
   for cont = x:-1:1
       if mod(x,cont) == 0 && mod(y,cont) == 0
           div = cont;
           break
       end
   end
else
   for cont = y:-1:1
       if mod(x,cont) == 0 && mod(y,cont) == 0
           div = cont;
           break
       end
   end
end