x = [0.13:0.001:10.13].';
plik = load('3.dat');
B = plik(:,1);
invB = (1./B);
I = plik(:,6);
U = plik(:,5);
R = U./I;
R_inv_int = interp1(invB, R, x, 'linear');
