freq = linspace(40e3,9e9,100000);
ckt_DUT = circuit('LC_Wire');
add(ckt_DUT,[1 2],inductor(1.25e-9))
add(ckt_DUT,[2 3],capacitor(1e-12))
add(ckt_DUT,[2 4],inductor(1.25e-9))
setports(ckt_DUT,[1 3],[4 3]);
ckt_C_open = circuit('Copen');
add(ckt_C_open,[1 3],capacitor(100e-15));
add(ckt_C_open,[1 2],resistor(1e9));
add(ckt_C_open,[2 3],capacitor(100e-15));
setports(ckt_C_open,[1 3],[2 3]);
ckt_L_short = circuit('Lshort');
add(ckt_L_short,[1 2],inductor(1e-9));
add(ckt_L_short,[2 3],inductor(1e-9));
add(ckt_L_short,[4 5],resistor(1e-15));
setports(ckt_L_short,[1 4],[3 5]);
S_ckt_DUT = sparameters(ckt_DUT,freq);
Z_ckt_DUT = s2z(S_ckt_DUT);