function f = funcao_hartley(x,y,fx,fy)
u = fx*x + fy*y + 1/8;
f = sin(2*pi*u);
