integration of 1d data

In [1]:

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import cumulative_trapezoid

x = np.arange(30)/3
y = np.sin(x)

integ = cumulative_trapezoid(y, x)

plt.plot(y, label = 'raw')
plt.plot(integ, label = 'integrated')
plt.legend()

Out[1]:

<matplotlib.legend.Legend at 0x7fe64b40caf0>