function y = funkce(x)
y =   3*(1-x(1))^2*exp(-(x(1)^2) - (x(2)+1)^2) ...
    - 10*(x(1)/5 - x(1)^3 - x(2)^5)*exp(-x(1)^2-x(2)^2) ...
   - 1/3*exp(-(x(1)+1)^2 - x(2)^2);
