import numpy as np
import matplotlib.pyplot as plt
from dtk.process import derivative

x = np.linspace(-10.0, 10.0, num=201)
y = x**2

fig, axes = plt.subplots(2, 1, layout='constrained')

axes[0].plot(x, y)
axes[1].plot(x, derivative(x, y, method='combination'))