SymPy from the ground up: where the pendulum's 1.74 % comes from

Afterwards you can substitute, differentiate, integrate, series-expand, and solve with SymPy, and turn a result into a NumPy function with lambdify.