Solution du Tp3

Solution de l'exercice:

function H=hessiantp3(f)

syms x y

a=[x y];

%il suffit de changer a=[x y z t] pour 4 variables et ainsi de suite

for i=1:size(a,1)

for j=1:size(a,2)

H(i,j)=diff(diff(f,1,a(i)),1,a(j));

end

end

pour le cacul de la matrice en un point donné

subs(H,{x,y},{1,1})

A partir de R2011b c'est plus facile de calculer la hessienne avec

Hessian(f)

jacobian(gradient(f))