# -*- coding: utf-8 -*-

import numpy as np
import scipy.integrate as integrate

# START CUT
result = integrate.odeint(
    lambda y, x: np.sin(x),
    0, np.linspace(0, np.pi, 10000)
)
# result est une matrice 2D (1 ligne par instant, 1 colonne par
# dimension du problème).
# result[-1] = dernière ligne.
# result[-1, 0] = 1er élément de la dernière ligne
print "Result =", result[-1, 0]
# END CUT