%typeset_mode False import scipy as sp from scipy import integrate from sage.symbolic.integration.integral import indefinite_integral import numpy as np import math