Vaja2

mu0=4*pi*1e-7;
B1=[0.1:0.1:1.3];
B2=B1/4;
% kA/m pomnozimo s 1000, da dobimo A/m
H=[0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4,5,6,7,8,10,12,14,16,18,20,22,24,26]*1000;
B=[0, 0.45, 0.635, 0.75, 0.84, 0.89, 0.93, 0.96, 0.98, 1.02, 1.05, 1.07, 1.1, 1.14, 1.17, 1.2, 1.23, 1.25, 1.27, 1.29, 1.31, 1.32];
figure(1);
plot(H,B,'.-', 'linewidth',2)

H1=[];
for Bx=B1
  %Bx=B1(7)
  i=sum(B<Bx);
  %B(i), B(i+1), H(i), H(i+1)
  Hx=H(i) + ( H(i+1) - H(i) ) / ( B(i+1) - B(i) ) *(Bx-B(i));
  H1=[H1, Hx];
end

H2=[];
for Bx=B2
  %Bx=B1(7)
  i=sum(B<Bx);
  %B(i), B(i+1), H(i), H(i+1)
  Hx=H(i) + ( H(i+1) - H(i) ) / ( B(i+1) - B(i) ) *(Bx-B(i));
  H2=[H2, Hx];
end
figure(2);
plot(H,B,'.-', H1, B1,'.-', H2, B2, '.-')

%Amperov zakon!
B2;
I=(H1*0.355+H2*0.03+B2/mu0*1.4e-3)/3000;
figure(3);
plot(I,B1, 'linewidth',2, I, B2, 'linewidth',2)