I can't get beam pattern plot of the code. The plot should similar to image attached

%d=lamda/2
clear all; clc;
A11=2*pi;
f1 = @(x) exp(-1i*pi/2*(sqrt(3)*sin(x)+cos(x))) ; A12 = integral(f1,0,2*pi);
f2 = @(x) exp(-1i*pi/2*(sqrt(3)*sin(x)-cos(x))) ; A13 = integral(f2,0,2*pi);
f3 = @(x) exp(1i*pi/2*(sqrt(3)*sin(x)+cos(x))) ; A21 = integral(f3,0,2*pi);
A22 = 2*pi;
f4 = @(x) exp(1i*pi*cos(x)) ; A23 = integral(f4,0,2*pi);
f5 = @(x) exp(1i*pi/2*(sqrt(3)*sin(x)-cos(x))) ; A31 = integral(f5,0,2*pi);
f6 = @(x) exp(-1i*pi*cos(x)) ; A32 = integral(f6,0,2*pi);
A33 = 2*pi;
f7 = @(x) (exp(-1i*pi/2*sqrt(3)*sin(x))) .* sinc(x-(pi/3)) ; b1 = integral(f7,0,2*pi);
f8 = @(x) (exp(1i*pi/2*cos(x))).* sinc(x-(pi/3)) ; b2 = integral(f8,0,2*pi);
f9 = @(x) (exp(-1i*pi/2*cos(x))).* sinc(x-(pi/3)) ; b3 = integral(f9,0,2*pi);
A = [A11 A12 A13; A21 A22 A23; A31 A32 A33];
B = [b1; b2; b3];
W = inv(A)*B;
W1 = W(1) W2 = W(2) W3 = W(3)
%Beam Pattern 
theta = 0:0.01:2*pi; phi=0:0.01:2*pi; y1=exp(1i*pi/2*sqrt(3)*sin(theta)); y2=exp(-1i*pi/2*cos(theta)); y3=exp(1i*pi/2*cos(theta));
E = W1*y1+W2*y2+ W3*y3; F = (cos((pi/2)*cos(phi)))/sin(phi); Radiation= E*F; modeval= abs(Radiation); power = modeval.^2;

[Shoaib Salim]

[Dyuman Joshi]

I don't see a command for plotting anything in your code.

What is the reference for the image shown in the comment above?