28 lines
906 B
Mathematica
28 lines
906 B
Mathematica
|
%-------------------------------------------------------------------------% %
|
||
|
% Fonction de calcul de la Jacobienne du residu de la fonction de %
|
||
|
% desintegration radioactive du carbone 14 %
|
||
|
%-------------------------------------------------------------------------%
|
||
|
|
||
|
function J_res = J_residu_C14(beta, donnees)
|
||
|
%
|
||
|
% Paramètres en entrés
|
||
|
% --------------------
|
||
|
% beta : vecteur des paramètres
|
||
|
% real(p)
|
||
|
% donnees : Données
|
||
|
% real(n,2)
|
||
|
%
|
||
|
% Paramètres en sortie
|
||
|
% --------------------
|
||
|
% J_res : Matrice jacobienne des résidus
|
||
|
% real(n,p)
|
||
|
%
|
||
|
|
||
|
% res = donnees(:, 2) - beta(1)*exp(-beta(2)*donnees(:, 1));
|
||
|
|
||
|
v1 = -exp(-beta(2)*donnees(:, 1));
|
||
|
v2 = beta(1)*donnees(:, 1).*exp( -beta(2)*donnees(:, 1) );
|
||
|
|
||
|
J_res = [ v1, v2 ];
|
||
|
|
||
|
end
|