function H=htest(x)

H =[
[12*(exp(x(1))-x(2)).^2.*exp(x(1)).^2+4*(exp(x(1))-x(2)).^3.*exp(x(1))+56*x(1).^6,        -12*(exp(x(1))-x(2)).^2.*exp(x(1)),                                                                                                             0,                                                                                                               0]
[                                              -12*(exp(x(1))-x(2)).^2.*exp(x(1)),12*(exp(x(1))-x(2)).^2+3000*(x(2)-x(3)).^4,                                                                                          -3000*(x(2)-x(3)).^4,                                                                                                               0]
[                                                                               0,                      -3000*(x(2)-x(3)).^4, 3000*(x(2)-x(3)).^4+12*tan(x(3)-x(4)).^2.*(1+tan(x(3)-x(4)).^2).^2+8*tan(x(3)-x(4)).^4.*(1+tan(x(3)-x(4)).^2), 12*tan(x(3)-x(4)).^2.*(1+tan(x(3)-x(4)).^2).*(-1-tan(x(3)-x(4)).^2)+8*tan(x(3)-x(4)).^4.*(-1-tan(x(3)-x(4)).^2)]
[                                                                               0,                                         0,12*tan(x(3)-x(4)).^2.*(1+tan(x(3)-x(4)).^2).*(-1-tan(x(3)-x(4)).^2)-8*tan(x(3)-x(4)).^4.*(1+tan(x(3)-x(4)).^2),                   12*tan(x(3)-x(4)).^2.*(-1-tan(x(3)-x(4)).^2).^2-8*tan(x(3)-x(4)).^4.*(-1-tan(x(3)-x(4)).^2)+2]];
