% Buler matrix array

N=input('What is the number of elements N? (2,4,8,...)');
d=.5;
th1=-pi/2:.01:pi/2;
ang=th1*180/pi;
figure;
for j=-(N-1):2:N-1
  AF=sin(N*pi*d*sin(th1)-j*pi/2)./(N*pi*d*sin(th1)-j*pi/2);
AF=abs(AF);
polarpattern(th1,abs(AF),'k')
hold on
end