function f = funcao_calombos(x,y)
N = 15;
f = 0;
R0 = 0.086;
X0 = -0.05;
Y0 = 0;
M = 6.5;
A = log(2)/M;
B = 2*pi/M;
for i = 1:N
    %calcula centro (x0,y0), largura L e altura h do calombo
    r = R0*exp((i-1)*A);
    t = (i-1)*B;
    x0 = r*cos(t) + X0;
    y0 = r*sin(t) + Y0;
    L = 0.40*r;
    h = 1.0/(2 - i/N);
    d = hypot(x-x0,y-y0);
    if d < L
        %acumula o valor do calombo no ponto (x,y)
        u = d^2/L^2;
            fi = h*(1 - u)^2;
            f = f + fi;
    end
end
    
    



