// 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);