function [B]=polje(R, rc)
% funkcija izracuna polje krozne zanke polmera rc pri radiju R
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




%R=0.1
fi=0; N=200; B=0;
dfi=2*pi/N;
dBB=0;
for i=1:N
   
    fi=fi+dfi;
    r2=rc^2+R.^2-2*rc.*R*cos(fi);
    r=sqrt(r2);
    theta1=acos((rc^2+r2-R.^2)./(2*rc.*r));
    theta=pi/2+theta1;
    dB=k*sin(theta)./r2;
   % dBB=[dBB dB]
   B=B+dB;
end

%plot(dBB(2:N))

% B0=mi0*I/(2*rc);