Scipy

from scipy.stats import probplot

probplot(residuals)


from scipy.interpolate import UnivariateSpline