辗转相处法求最大公因数及其证明

% 辗转相处法求最大公因数
% 定理：gcd(a,b) = gcd(b,a mod b)
% 证明
% a可以表示成a = kb + r，则r = a mod b
% 假设d是a,b的一个公约数，则有
% d|a, d|b，而r = a - kb，因此d|r
% 因此d是(b,a mod b)的公约数
%
% 假设d 是(b,a mod b)的公约数，则
% d | b , d |r ，但是a = kb +r
% 因此d也是(a,b)的公约数
%
% 因此(a,b)和(b,a mod b)的公约数是一样的，其最大公约数也必然相等，得证

function nResult = my_gcd(m, n)
    % 调整次序
    if m < n, swap(m, n),end
    
    % 第一次
    nResult = mod(m, n);
    last_n = n;

   

    % 不使用递归

    while nResult > 1
        last_n = nResult;
        nResult = mod(n, nResult);
        n = last_n;
    end    

    if nResult == 1

       % 没有最大公因数,1为他们的最大公因数

        nResult = 1;

    else

        % 找到

        nResult = last_n;
    end
end