function H=hessian(f,x,x0)

% H=hessian(f,x) 计算表达式
% H=hessian(f,x,x0) 计算hessian矩阵的值 x0为x的初值

%$copyright by$ LUO sir

switch  nargin
    case 1
    error('please input variables in f(x)')

    case 2
        H=subhessian(f,x);
 
    case 3
        H=subhessian(f,x,x0);
        x=x0;
        H=subs(H);
    otherwise
       
         error('too many arguments or nothing')         
end

function HH=subhessian(f,x,x0)
 
n=length(x);

J=jacobian(f,x);
HH=[];
HH=sym(HH);
for i=1:n
HH(i,:)=jacobian(J(1,i),x);
end