%%%%%%%%%%%%%%%%%%%%%%
% angular_function.m %
%%%%%%%%%%%%%%%%%%%%%%

clf;

ALPHA = 2.58128140043517;
BETA  = 10;

theta = 0: .01 : pi;
y =  ( (1 + cos (theta) ) / ALPHA ) .^ BETA;

f1 = exp ( - sqrt ( y ) );
f6 = exp ( - sqrt ( 6 * y ) );
f8 = exp ( - sqrt ( 8 * y ) );
f14 = exp ( - sqrt ( 14 * y ) );

plot (theta * 180 / pi, f1, 'c');
axis ([0 120 0.3 1]);
hold on;
plot (theta * 180 / pi, f6, 'm-.');
plot (theta * 180 / pi, f8, 'r--');
plot (theta * 180 / pi, f14, 'g:');
legend ('e^{-{g(\theta)}^{0.5}}', ...
        'e^{-{[6g(\theta)]}^{0.5}}', ...
        'e^{-{[8g(\theta)]}^{0.5}}', ...
        'e^{-{[14g(\theta)]}^{0.5}}', 4);

xlabel('\theta [degree]');
ylabel('screening factor');
title(['Representative screening functions with '...
       'g(\theta) = ((1+cos\theta)/\alpha)^{\beta}' ]);
grid on;

set(gca,'XTick',0:30:180);

print -depsc angular_function.eps;
