Matlab code for Binomial array antenna

clc;
clear all;
close all hidden ;
f=input('please enter the operating frequency in Hz:\n->');
N=input('please enter the number of elements(N):\n->');
lambda=(3e8)/f%lambda=c/f% the value is in meters
d=input('please enter the distance d in meters between the elemnts:\n->');
k=2*pi/lambda;

theta=0:0.01:2*pi;
U=k*d*cos(theta);%the U here is twice the u used in the block.
A=[0 0 0 0 0;1 0 0 0 0;1 1 0 0 0;3 1 0 0 0;3 4 1 0 0;10 5 1 0 0;10 15 6 1 0;...
   35 21 7 1 0;35 56 28 8 1;126 84 36 9 1];
x=mod(N,2);
if x==0;
    AF=0;
    for n=1:5;
        AF=AF+A(N,n).*cos(((2*n-1)/2).*U);
    end
end
W=abs(AF);
w=W/max(W);
polar(theta,w)
title('normalized eE-field of array antenna inlenera scale');
figure;
afdb=20*log10(W);
afplot=(afdb+abs(afdb))/2;
polar(theta,afplot);
title('non-normalized polar pattern of the array in dB scale');
HPBW=(1.06/(N-1)^0.5)*(180/pi)
Do=1.77*(N)^0.5
Do_db=10*log10(5.5972)