z1=80;
z2=80;
z3=60;
for f=1:3
d1=pi/4*2./f;
d2=pi/4*2./f;
d3=pi/3*2./f;
Z0=50;
Z=-j*z3*cot(d3);
Y=1./Z;
Y0=1/z1;
Y1=1/z2;
......................