A=1e-4;
l=0.24;
I=1;
N=1000;
H1000=N*I/l

mi0=4*pi*1e-7;
H=0:10:5500;



plot(H,BB(H)); xlabel('H  / A/m'); ylabel('B   / T')

figure; plot(H,B./(mi0*H)); xlabel('H  / A/m'); ylabel('mi r')
break
% lineariziram pri H=1000
B1000=BB(1000);
mir1000=B1000/(mi0*1000)
Fluks=B1000*A
I1000=1000*l/N;
L=N*Fluks/I1000

plot(H,ones(length(H))*L)
hold on
% racunam z magnetilno krivuljo
B=BB(H);
Fluks=B*A;
II=H*l/N;
L=N*Fluks./II
plot(H,L); xlabel('H  / A/m'); ylabel('L   / H')


% Z zracno rezo
lz=0.5e-3;
Bz1000=mir1000*mi0*1000;
Fluks=Bz1000*A
I1000=(1000*l+Bz1000/mi0*lz)/N;
L=N*Fluks/I1000
plot(H,ones(length(H))*L,'--')

% z rezo in upostevanjem magnetilne krivulje
% racunam z magnetilno krivuljo
B=BB(H);
Fluks=B*A;
II=H*l/N+B*lz/mi0/N;
L=N*Fluks./II
plot(H,L,'--'); xlabel('H  / A/m'); ylabel('L   / H')