% Smart antennas 3-D pattern for loop antennas with radius a
% use 100 data points in theta and 100 data points in phi

set(0,'defaultfigurecolor','w')
tend=pi;
%a is inserted below before the pi as 2*a*pi;
% this case is for a = .2
fx=inline('((cos(pi*cos(theta))-cos(pi))/sin(theta))^2*sin(theta)*cos(phi)');
fy=inline('((cos(pi*cos(theta))-cos(pi))/sin(theta))^2*sin(theta)*sin(phi)');
fz=inline('((cos(pi*cos(theta))-cos(pi))/sin(theta))^2*cos(theta)');
figure;
%ezsurf(fx,fy,fz,[pi/2,3*pi/2,0,tend],100)
ezmesh(fx,fy,fz,[pi/2,2*pi,0,pi],100)
%shading interp
colormap([0 0 0])
%brighten(.5)
axis square
axis equal
%axis([-.5 .5 -.5 .5 0 1])
%set(gca,'XDir','reverse','YDir','reverse')
view(121,24)
%lightangle(-63,24)
%camlight right
%camlight(20,-20,'infinite')
%material dull
%lighting phong
grid off
axis off
title .
zoom(1.5)