integrate (2*x+y,(x,1,3),(y,2,4))

integrate(2*x + y, (x, 1, 3), (y, 2, 4))

len(str(28))

factorint(28, limit=100)

factorint(28, limit=256)

Experiment with SymPy at SymPy Live.