Quick integration of a known function

Using scipy.integrate