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

N=101;      % stevilo diskretiziranih odsekov vzdolz osi
rmax=2*r0; % koncni 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 za vrednosti r<r0
B_zunaj=polje(r_zunaj,r0);  % izračun polja za vrednosti r>r0

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

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

% Primerjam vrednosti polja pri r=0: na osi
disp('B v osi numericno:'); B_znotraj(1)  
disp('B v osi analiticno:');B0=mi0*I/(2*r0)

hold on
plot(0,B0,'*r') % prikazem z zvezdico