clear;

s = 0.01;
l = 0.1;

theta = 0:1:80;
np = length(theta);

eps = (9 - 1i*0.9);

for n = 1:np
    [e_v(n) e_h(n)] = SemiEmp_EmissivityModel_RoughSurf(eps, s,l,theta(n));
end

figure(1)

plot(theta, e_v, '-r', theta, e_h,'-b')
xlabel('\theta (deg)')
ylabel('Emissivity')
legend('e_v', 'e_h', 3)
grid
