function v = funcao_cratera(x,y)
S = 0.7;  %fator de escala geral 
R0 = 0.40*S;
h = 0.04*S;
R1 = 0.15*S;
FX = 0.3*pi/R1;
FY = 0.4*pi/R1;

r2 = x^2 + y^2;
%pico central
%v1 = exp(-5*r2/(R1*R1))*cos(FX*x + FY*y);
v1 = exp(-5*r2/(R1*R1));
%parede circular 
z = (r2 - (R0*R0))/(2*R0);
v0 = exp(-0.5*(z/h)^2); 
v = 0.5*v0 + v1;




% %original
% function v = funcao_cratera(x,y)
% S = 0.7;  %fator de escala geral 
% R0 = 0.35*S;
% h = 0.03*S;
% R1 = 0.11*S;
% FX = 0.3*pi/R1;
% FY = 0.4*pi/R1;
% 
% r2 = x^2 + y^2;
% z = (r2 - (R0*R0))/(2*R0);
% v0 = exp(-0.5*(z/h)^2);
% v1 = exp(-0.5*r2/(R1*R1))*cos(FX*x + FY*y);
% v = 0.5*v0 + v1;