function cr = reflexion_vertical(e,d,angulo_grados)
%Función que calcula el coeficiente de reflexión para polarización
%vertical

ang=angulo_grados*(pi/180);
f=2.4*10^9;
c=3*10^8;

numerador = e*sin(ang) - sqrt(e - cos(ang)*cos(ang));
denominador = e*sin(ang) + sqrt(e - cos(ang)*cos(ang));
rv = numerador/denominador;

numerador = 2*pi*f*sqrt(e - sin(ang));
k=numerador/c;

numerador=1-exp(-i*2*k*d);
denominador=1-(rv^2)*exp(-i*2*k*d);
cr=rv*(numerador/denominador);
end