// Programme d'étude de la décroissance radioactive
// Dominique Lefebvre septembre 2010
// TangenteX.com

// Définition des vecteurs de mesure
t = [0,10,20,30,40,50,60,70,80,90,100,110,120];
N = [6145,5470,4825,4210,3750,3345,2955,2605,2260,2015,1795,1575,1390];

// Tracé de la courbe de mesure N = f(t)
plot2d(t,N, style = 2);

// recherche d'une courbe de régression exponentielle
// complètement hors programme!
n = length(t);
T = [ones(n,1),t'];
LN = log(N);
LNT = LN';
TT = inv(T'*T)*T'*LNT;
a = exp(TT(1,1));
b = TT(2,1);
Y = a*exp(b*t);

// affichage de la courbe de régression en rouge et des paramètres a et b
mprintf('a: %f b: %f\n',a,b);
plot2d(t, Y, style = 5);