1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
from pylab import *
from scipy import integrate

L = 5                  # constante de décroissance

def derivative(y, t0):
    return -L * y      # dN/dt = -L * N, décroissance radioactive

N = 1000               # valeur à t = 0
t = arange(0, 1, 0.01)                  # intervalle et dates

nt = integrate.odeint(derivative, N, t) # intégration, voir la doc de scipy
print nt

plot(t, nt[:,0])                        # on extrait la première colonne du tableau 2D
show()