function [MCD]=mcd(a,b)

if a==b
   MCD=a ;
elseif a>b
   MCD=mcd(a-b,b) ;
else
   MCD=mcd(a,b-a) ;
end

fprintf('... chiamata\n')

end
