function BP=beampattern(c,dlam,phivec)

%%% Vector method
phivec=phivec(:)';			% make sure it's a row vector
M=length(c);
mvec=(0:M-1)';
uvec=dlam*sin(pi*phivec/180);
V=1/sqrt(M)*exp(-j*2*pi*mvec*uvec);
BP=c'*V;

%%% FFT method
% Nfft=1024;
% BP=fftshift(fft(conj(c),Nfft));
% w=2*pi*(0:Nfft-1)/Nfft;		
% mid=ceil(Nfft/2)+1;
% w(mid:Nfft)=w(mid:Nfft)-2*pi;
% u=fftshift(w)/(2*pi);
% phi=asin(1/dlam*u);