theta_angle=[30 50];
theta_angle;

%% List of the angles in radian
theta=[];
theta = theta_angle * (pi / 180);
theta;

phi=[];
phi = pi * sin(theta);
phi;

%% Filling matrix A which is N x L
A = zeros(N,L);
for l=0:(L-1),
        for k=0:(N-1),
                A(k+1,l+1) = exp(i * k * phi(l+1));
        end;
end;