[f,p,ang,temp] = textread('kertoimet.txt');

c=4.95e-12;
l=1.59e-9;
zl  = (p+1)./(1-p);
z_c = 1./(1i*2*pi.*f*l);
z_l = 1i*2*pi.*f*l;

z=(zl.^-1+z_c.^-1).^-1+z_l;

p_new = (z-1)./(z+1);