// Simulation de décroissance radioactive
// Dominique Lefebvre septembre 2010
// TangenteX.com

// Définition de la fonction à intégrer 
// dN = -lambda*N*dt
function du = Desintegration(t,u)
du = -lambda*u;
endfunction

// Saisie des paramètres de la simulation
N0 = input('Population initiale des noyaux (N0): ');
lambda = input('Constante de desintegration (en s-1): ');
t1 = input('Durée de la mesures (en s): ');

// Définition des conditions initiales
t0 = 0.;
u0 = [N0];

// Autres paramètres de la simulation
Pas = 1000;
dt = (t1 - t0)/Pas;

// Initialisation du vecteur temps
t = t0:dt:t1;

// intégration 
u = ode(u0,t0,t,Desintegration);

// Affichage
mprintf('Demi-vie: %f\n',log(2)/lambda);
clf;
xtitle('Décroissance radioactive', 'Temps (s)', 'Population');
plot2d(t,u,style = 2);