Matlab求两个数最大公约数

在命令窗口输入两个自然数，输出它们的最大公约数（不要用matlab自带的求最小公倍数和最大公约数的函数。利用求同余的函数判断是否整除）

function greatestcommondivisor()

%求两个数的最大公约数

a=input('请输入正整数a:');

b=input('请输入正整数b:');

%先求把两个数中最小的值赋给tem

if a>b

tem=b;

else tem=a;

end

%进行步长为1的倒序循环，保证求得的公约数是最大的

for n=tem:-1:1

%c和d分别是a、b除以m的余数，c、d同时为零时，m即为最大公约数c=rem(b,n);

d=rem(a,n);

if c==0&d==0

disp(['最大公约数为;',num2str(n)])

break;

end

end