require scipy>=1.6 because of scipy.integrate.simpson

require scipy>=1.6 because of scipy.integrate.simpson #4