d = 0.01;        %distance, meters
C = 0.5*10^-12; %capacitance, farads
Zo = 50;        %Zo, ohms
eps_eff = 5;    %Epsilon eff value


frequency_array = linspace(.3,3000, 9999).*10^6;
omega_array = frequency_array.*2.*pi;
beta = 2.*pi./lambda;

lambda = (2.98*10^8) ./ (sqrt(eps_eff).*frequency_array);


gamma_lower = 1i.*beta;

gamma_upper = 1     ./  (   1 + (  (2 .* Zo)./(1i .*  C .* omega_array) )   );

Z1 = Zo .* (1 + gamma_upper .* exp(-2.*gamma_lower.*d)) ./  (1-gamma_upper .* exp(-2.*gamma_lower.*d));

VoVi = abs (...
        (   (1+gamma_upper)    ./   (1+ gamma_upper .* exp(-2.*gamma_lower.*d)  )   ) ...
        .* (    (1i.*omega_array.*C.*Z1) ./     (1 + 1i.*omega_array.*C.*Z1)    ) ...
        .*  ((1i.*omega_array.*C.*Zo)    ./     (1 + 1i.*omega_array.*C.*Zo)) ...
        .* exp(-2.*gamma_lower.*d)...
        );
VoVidB = 10 .* log10(VoVi);

plot (frequency_array, VoVidB)