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

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

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

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

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