[theta,phi]=meshgrid(0:0.02*pi:2*pi,0:pi/12:1.1*pi);

%E=cos(theta);
the=pi/2-theta;
%phi=ones(size(phi))*pi/4;
% E=(cos(3/2*pi*cos(the)))./sin(the);
% E=cos((pi/2)*cos(the))*sin(the);

% if sin(the)==0
%     E=-i*S.*(1+cos(the))/2
% elseif sin(phi)==0
%     E=-i*S.*(1+cos(the))./2*sin(psi1)./psi1
% elseif cos(phi)==0
%     E=-i*S.*(1+cos(the))./2*sin(psi2)./psi2
% else
Fac=[0,0,0;0,1,0;0,0,0];
%Fac=[1,1,1;1,1,1;1,1,1];
%Fac=[1,2,1;2,4,2;1,2,1];

% end
AE=platecomp(the,phi,Fac);
[x,y,z]=sph2cart(phi,theta,AE);
%z=abs(z);
mesh(x,y,z,'EdgeColor','black');
axis off
axis equal

%polar(theta(1,:),AE(1,:),'-k')
axis off