zero_crossings = InterpolatingUnivariateSpline_function.derivative().roots() extrema_values = InterpolatingUnivariateSpline_function(zero_crossings) zero_values = np.zeros(len(zero_crossings)) # generate an array of zeros