function [gcd] = myGCD(a,b)
  if b == 0
    gcd = a;
    return
  else
    gcd = myGCD(b,mod(a,b));
    return
  end
end