clc; clear
x = 0.399;
b = 2.84e-6;
l = 50;
% Zc = 2*sqrt(x/(4*b - x*b^2*l^2))
Zc = sqrt(x/b)

Zx = 1j*x*l
Zb = 1/(1j*b*l/2)
par = @(Z1, Z2) Z1*Z2/(Z1+Z2);
Zvlezno = par(par(Zb, Zc) + Zx, Zb);
fprintf('Zvlezno = %.4f + j%.4f\n', real(Zvlezno), imag(Zvlezno));


