function f = funcao_calombo(x,y)

x1 = x - 0.25;
y1 = y - 0.20;
R = 0.15;
u = (x1^2+y1^2)/R^2;
if u < 1 
    f = (1 - u)^2; 
else 
    f = 0;
end