# -*- coding: utf-8 -*-
import scipy.integrate
import numpy as np
import matplotlib.pyplot as plt

# y' = f(y, x)
def f(y, x):
    return np.cos(10 * y) * np.cos(x)

x = np.linspace(0, 10, 100)
# Résolution numérique de l'équation
y = scipy.integrate.odeint(f, 0, x)

# Tracé du résultat
plt.plot(x, y)
plt.show()