function H = hess_ex5(x,iex)
% Define gradients for exercise 5

if nargin==1
    iex=1;
end

if iex==1
    H=zeros(2,2);
    %g(1)=4*(x(1)-1)+1;
    %g(2)=2*(x(2)-2)+1;
    H(1,1)=4; H(1,2)=0;
    H(2,1)=0; H(2,2)=2;
elseif iex==2
    H=zeros(2,2);
    %g(1)=-400*(x(2)-x(1)^2)*x(1)-2*(1-x(1));
    %g(2)=200*(x(2)-x(1)^2);
    H(1,1)=-400*(x(2)-3*x(1)^2)+2; H(1,2)=-400*x(1);
    H(2,1)=-400*x(1); H(2,2)=200;
end