function hesse

%Anlegen der symbolischen Variablen
syms 'x1';
syms 'x2';
syms 'x3';
var=[x1;x2;x3];

%Eingabe der Funktion
f='x1 * x3 * exp(-x1^2 - x2^2)';

%Berechnen des Gradienten
abl=diff(f,var(1));
abl2=diff(f,var(2));
abl3=diff(f,var(3));
nabla=[abl;abl2;abl3]
abl=char(abl);
abl2=char(abl2);
abl3=char(abl3);

%Berechnen der Hessematrix
for i=1:3
    for j=1:3
        hess(i,j)=diff(nabla(i),var(j));
        %hess(i,j)=char(hess(i,j));
        %hess(i,j)=inline(hess(i,j));
    end;
end;
hess

%Einsetzen des Punktes in die Hessematrix
x1=2;
x2=6;
x3=4;
for i=1:3
    for j=1:3
        hess(i,j)=feval(char(hess(i,j)),x1,x2);
    end;
end;
hess