function gl=ex321(a,b)
    m=a;n=b;k=1;
    while(k~=0)
        k=mod(m,n);
        m=n;n=k;
    end
    gl=[m,a*b/m];