plotting dipole antenna far-field pattern, with multiple slices

i'm attempting to plot the far field pattern for multiple dipoles in a line, which make up one large one,

however, i need an input for current for each slice, so they can all be plotted on top of eachother so it appears as one dipole.

however, i am running into bugs in ths program, and i need it to ot only calculate the entire thing up at one point, but to do it in a circle so it can plot the field pattern

deltaL=.003366 %meters

beta=12.5664 %

eta=377

i=1.5 %current

betasq=157.914 %beta squared

addpath c:\antennas\ewa

l=.68

a=.001

ker='a'

basis='d'

R=1 %distance, one meter from dipole

E=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,100/.001,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

pfield(l,a,E,ker,basis)

%this spits out a large complex vector of many rows of the current

for k=0:1:49

i=abs(ans(k))%this will extract a different rom/value from the answer for each count

E=-x*(377*i*betasq*deltaL)/(2*pi)*e^(-sqrt(-1)*beta*R)*cos(theta)(1/(sqrt(-1)*beta*r)^2+1/(sqrt(-1)*beta*r)^3)-y*(377*i*betasq*deltaL)/(4*pi)*e^(-sqrt(-1)*beta)*R*sin(theta)*(1/(sqrt(-1)*beta*r)+1/(sqrt(-1)*beta*r)^2+1/(sqrt(-1)*beta*r)^3)

E=E*2 %second half of dipole

plot(x,y)

hold on

end