function f = funcao_doispicos(x,y)

x1 = x - 0.25;
y1 = y - 0.20;
R1 = 0.19;

u1 = (x1^2+y1^2)/R1^2;
if u1 < 1 
    f1 = (1 - u1)^2; 
else 
    f1 = 0;
end


x2 = x + 0.25;
y2 = y + 0.20;
R2 = 0.03;
u2 = (x2^2+y2^2)/R2^2;
if u2 < 1 
    f2 = (1 - u2)^2; 
else 
    f2 = 0;
end
f = f1 + f2;