function [u,v] = azimuthal()
%azimuthal.m A function samples a random direction in 2D space
% 
% FileName      : azimuthal.m
% Creation      : 09/29/2016
% Modifications : September 29, 2016 Initial creation
% Author(s)     : Alex Bielajew
%
% Inputs        : [] i.e. none
%      
% Outputs
%  double
%     u :  x-axis direction cosine
%     v :  y-axis direction cosine
%
% Purpose       : NERS 544 Monte Carlo Methods        
%               : Demonstration, production, part of the code library

phi = 2*pi*rand; % Random azimuthal direction in radians
u = cos(phi);    % x-axis direction cosine
v = sin(phi);    % y-axis direction cosine

return