active elemet impedance

Hi all,
i am trying to calculate the active element impdenace with MATLAB for planar phased array consisting of circular waveguide.
Does anyone have MATLAB code about this topic?
thank you very much!