from __future__ import division
from scipy import *

a = 44551
b = 455791

print "This program calculates the greatest common divisor of a and b."
print " a =", a
print " b =", b

# Make sure that b is smaller than a
if a < b:
    # Swap a and b
    temp = a
    a = b
    b = temp

while b > 0:
    a = a % b
    # Swap a and b
    temp = a
    a = b
    b = temp

print " gcd(a,b) =", a

