radiation pattern for square loop antenna

how i can

Write a Matlab code that analyzes a square loop antenna with half-wavelength side and fed at the mid-point of one of the sides. Sinusoidal current distribution is assumed.
Code Output:
Radiation pattern in the loop plane (the XY plane).
Radiation pattern in any plane phi=constant
The radiation resistance

i need it too