f=10^(3).*[20 30 50 70 90 100];c=10^(-9).*[313.6 264.2 134.2 86.2 58.80 49.3];
r=[12.5 16.1 25.6 38.3 55.2 63.3];p=0.06;c1=10^(-9)*3.6*[1 1 1 1 1 1];
for k=1:length(f)
zsc(k)=r(k)/(1-i*2*pi*r(k)*c(k)*f(k));

x(k)=real(zsc(k));
y(k)=imag(zsc(k));
l(k)=(1/(2*pi*f(k)))*y(k)*(1/p);
r1(k)=1/p*x(k);

ygol(k)=i*2*pi*f(k)*c1(k);
zc(k)=sqrt(zsc(k)/ygol(k));
gamma(k)=(1/p)*sqrt(zsc(k)*ygol(k));
end
%plot(l)
plot(f,r1)