Create gcd.py

Greatest common divisor and factors of n.
by JRG