SymPy

Simple example

Nevertheless, let's put here a simple example:

from sympy import Integral, Symbol, pprint

x = Symbol("x")
y = Symbol("y")

e = 1/(x + y)
s = e.series(x, 0, 5)

print(s)
pprint(s)

This example should print the following after the execution:

1/y + x**2*y**(-3) + x**4*y**(-5) - x*y**(-2) - x**3*y**(-4) + O(x**5)
     2    4         3
1   x    x    x    x
─ + ── + ── - ── - ── + O(x**5)
y    3    5    2    4
    y    y    y    y