function [f,gf,Hf]=testfunc03(x)

f=-10*x(1)^2*x(2)+5*x(1)^2+4*x(2)^2+x(1)^4+2*x(2)^4;

gf=[-20*x(1)*x(2)+10*x(1)+4*x(1)^3
    -10*x(1)^2+8*x(2)+8*x(2)^3];

Hf=[-20*x(2)+10+12*x(1)^2   -20*x(1);
    -20*x(1)   8+24*x(2)^2];