function Z = l24f(X,Y)
Z = 10*(Y-X.^2).^2 + (1-X).^2;
% = 10 y^2 -20 yx^2 + 10x^4 + 1 -2x+x^2
