% Smart antennas 3-D pattern for loop antennas with circumference C/lambda.
% use 100 data points in theta and 100 data points in phi

set(0,'defaultfigurecolor','w')
tend=pi;
%Change C for each required plot
C=3;
fx=inline('besselj(1,3*sin(theta))^2*sin(theta)*cos(phi)');
fy=inline('besselj(1,3*sin(theta))^2*sin(theta)*sin(phi)');
fz=inline('besselj(1,3*sin(theta))^2*cos(theta)');
maxx=besselj(1,3)^2;
figure;
ezmesh(fx,fy,fz,[0,pi,0,2*pi],100)
colormap([0 0 0])
axis([-maxx maxx -maxx maxx -maxx maxx])
axis equal
view(0,24)
grid off
axis off
title .
%zoom(1.1)