How to code for E-plane and H-plane of rectangular patch antenna in MATLAB?

Have this so far but how do i plot to look like the image?

clear;

clc;

%% Patch Dimensions and Parameters

dieconst=2.2;       % dielectric constant

h=0.1588;           % hieght in cm

f0=10e9;            % resonate in Hz

w=1.186;            % width in cm

effdie=1.972;       % effective dielectric constant

inlength=0.081;     % incremental length in cm

l=0.906;            % length in cm

Le=1.068;           % effective length in cm

E0=1;

v0=h*E0;

k0=(2*pi*f)/(3e8);  % wavenumber

r=1;

phi=linspace(0,pi/2,1e3);

theta=pi/2;

%% Clean up

a=j*((k0*w*v0*e.^-j*k0*r)/(pi*r));

b=((k0*h)/2)*cos(phi);

c=((k0*h)/2)*sin(theta);

d=((k0*Le)/2)*sin(phi);

e=((k0*w)/2)*cos(theta);

%% E-Plane

E=a*[sin(b)/b]*cos(d);

%% H-Plane

H=a*[sin(theta)*(((sin(c)/c))*((sin(e)/e)))];

%% Plot

-E. (E-plane) E, (H-plane) 150 180 270 0 300 -15 E- and H-plane Patterns of Rectangular Microstrip Antenna 120 330 -30 -45 8-90 b=0 dB -45 30 -30 60 -15 60 0 30 90 0