import sympy as sym
import numpy as np
import math
from IPython.display import display, Math
from sympy.abc import w,x,y,z,a,b,c,d
sym.init_printing()
import math
math.gcd(4,16)
math.gcd(95,100)
# application
a = 16
b = 88
fact = math.gcd(a,b)
display(Math('\\frac{%g}{%g} = \\frac{%g\\times %g}{%g\\times %g}' %(a,b,a/fact,fact,b/fact,fact)))