%% Ukazi za izračun polja prečno na os zanke pri z=0
global k
mi0=4*pi*1e-7;
r0=15e-2;  % polmer kratke tuljave
I=1
k=mi0*I/(4*pi)

N=101;      % število diskretiziranih odsekov
rmax=2*r0; % končni radij
r_znotraj=linspace(0,0.9*r0,N)  % vektor polmerov razdeljen na N linearno razporejenih vrednosti
r_zunaj=linspace(1.1*r0,rmax,N)  % vektor polmerov razdeljen na N linearno razporejenih vrednosti


B_znotraj=polje(r_znotraj,r0)  % izračun polja pri vrednostih r
B_zunaj=polje(r_zunaj,r0)  % izračun polja pri vrednostih r

plot([0 rmax],[0 0],r_znotraj,B_znotraj,r_zunaj,B_zunaj, 'Linewidth',2);

xlabel('Polmer   / m')
ylabel('Gostota magnetnega pretoka  / T')