function [H_f, res, J_res] = Hess_f_C14(beta, donnees, residu, J_residu)
%
% Paramètres en entrés
% --------------------
% beta     : vecteur des paramètres
%            real(p)
% donnees  : Données
%            real(n,2)
% residu   : fonction qui code les résidus
%            res_beta = residus(beta)
% J_residu : fonction qui code la matrice jacobienne
%            J_res_beta = J_residu(beta);
%
% Paramètres en sortie
% --------------------
% H_f      : Matrice hessienne
%            real(p,p)
% res      : vecteur des résidus
%            real(n)
% J_res    : Matrice jacobienne des résiduis
%            real(n,p)
%

    res = residu(beta);
    J_res = J_residu(beta);
    % cf formule du cours
    S = [ sum(0*donnees(:,1))                             , sum(donnees(:,1).*exp(-beta(2)*donnees(:, 1)));
              sum(donnees(:, 1).*exp(-beta(2)*donnees(:, 1))), sum(-beta(1)*(donnees(:, 1).^2).*exp(-beta(2)*donnees(:, 1))) ];
    H_f = S + J_res.'*J_res;
    
end