function [total, specular, volume] = exponential_bscatter(x,diel,s,A)
%Input parameters:
%   diel = dielectric constant
%   rms = rms height (meters)
%   L = correlation length (meters)
%   theta = incidence angle (radians)

if x >= pi/2
    total = NaN;
    return
end

ref = ((1 - sqrt(diel))/(1 + sqrt(diel)))^2; %Fresnel reflection coefficient
specular = ( (3.*(ref)) ./ ((s.^2).* (cos(x).^4) )) .* exp( (-(sqrt(6) .* tan(x)) ) ./ ((s) ));
volume = volume_scatter(A, x);
total = specular + volume;

end
