function[]=maxcd(m,n)
b=0;
c=0;
for a=1:n
    for b=1:m
        if a==b
            c=a;
            if rem(n,a)==0
                if rem(m,a)==0
                f=[c];
                for g=1:f
                    g=f;
                end
            end
        end
    end
end
end
disp('massimo comun divisore')
disp(g)