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

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

"Get 15% discount on your first 3 orders with us"
Use the following coupon
"FIRST15"

Order Now