% Smart antennas figure 3.7. 3-D pattern for cos(theta)^4 pattern
% use 100 data points in theta and 100 data points in phi

tend=pi/2;
set(0,'defaultfigurecolor','w')
fx=inline('abs(sin(3*pi*sin(theta))/(3*pi*sin(theta)))*sin(theta)*cos(phi)');
fy=inline('abs(sin(3*pi*sin(theta))/(3*pi*sin(theta)))*sin(theta)*sin(phi)');
fz=inline('abs(sin(3*pi*sin(theta))/(3*pi*sin(theta)))*cos(theta)');
figure;
ezsurf(fx,fy,fz,[pi,-pi,0,tend],100)
shading interp
colormap(gray)
brighten(.5)
axis square
axis equal
axis([-.5 .5 -.5 .5 0 1])
set(gca,'XDir','reverse','YDir','reverse')
view(-63,24)
camlight(20,-20,'infinite')
material dull
lighting phong
grid off
axis off
title .