function f = onda_hartley(fx,fy,p,CD,LD)
np = size(p,1);
UD = CD - LD/2;  % canto inferior do dominio
escala = sqrt(2/np);
for i = 1:np
    x = (p(i,1)-UD(1))/LD(1);  
    y = (p(i,2)-UD(2))/LD(2);
    fase = (fx*x + fy*y + 1/8)*2*pi;
    f(i) = escala*sin(fase);
end