
f1=[1000 4000 10000 20000 40000 100000 150000 200000 400000 500000 1000000];

f=4000:4000:1060000;d=4000;
r1=[268 268 268 269  271 282 295 312 390 425 603];
l1=10^(-6).*[680 680 678 675 669 650 642 635 619 608 574];
%c5=polyfit(f1,r1,3);

%c15=polyfit(f1,l1,3);
%r=polyval(c5,f);
%l=polyval(c15,f);
r=spline(f1,r1,f);
l=spline(f1,l1,f);
c=10^(-9)*45.5.*ones(size(f));
g=zeros(size(f));z1=zeros(size(f));y=zeros(size(f));d1=3*.3048; d5=6*.3048;d7=1.5*.03048; %ultima portiune;
for k=1:length(f);
z1(k)=r(k)+2*pi*f(k)*l(k)*i;
y(k)=g(k)+2*pi*f(k)*c(k)*i;
zc(k)=sqrt(z1(k)/y(k));
gamma(k)=sqrt(z1(k)*y(k));
alfa(k)=real(gamma(k));
beta(k)=imag(gamma(k));
