Q = posets.DivisorLattice(factorial(7))
Q.plot()

W = Q.subposet(flatten(Q.level_sets()[0:3]))
W.plot()

gcd(37,200)

xgcd(37,200)

-27*37+5*200

gcd(4,6,10)

---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-8-58a5f6bfd37d> in <module>()
----> 1 gcd(Integer(4),Integer(6),Integer(10))

TypeError: gcd() takes at most 2 arguments (3 given)

gcd(gcd(4,6),10)

?diophantine

Object `diophantine` not found.

P=Polyhedron(eqns=[(-28,4,6)],ieqs=[(0,1,0),(0,0,1)],base_ring=QQ)

P.integral_points()

var("x,y"); implicit_plot(4*x+6*y==28,(x,0,10),(y,0,10)) + list_plot(P.integral_points(),color='red')