python/sympy_cheat/logarithm_integral.py
#!/usr/bin/env python3
from sympy import *
x = symbols('x')
myli = integrate(sympify(1)/ln(x), x)
# It recognizes our definition as its own li! Beauty.
assert myli.equals(li(x))
for r in range(-2, 2):
for i in range(-2, 2):
print(f'{r} {i} {li(r + i*I).evalf()}')