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?

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