function [B]=polje(R, rc)
% funkcija izracuna polje krozne zanke polmera rc za oddaljenost r od osi
global k

if (R==rc)
    error('Pri polmeru zanke B ni definiran')
end
if (rc==0)
    error('Polmer zanke ne more biti enak 0')
end



fi=0; N=500; B=0;
dfi=2*pi/N;  
dl=2*pi*rc/N;
dBB=0;
for i=0:N  % INTEGRACIJA DO 2*PI
   
    fi=fi+dfi;
    r2=rc^2+R.^2-2*rc.*R*cos(fi);
    r=sqrt(r2);  % razdalja od tokovnega elementa do tocke
    theta1=acos((rc^2+r2-R.^2)./(2*rc.*r));
    theta=pi/2+theta1; % kot med dl in r
    dB=k*dl*sin(theta)./r2;  % izracun po B-S za dB
   % dBB=[dBB dB]
   B=B+dB;   % superpozicija
end
B=B;