% Izris polja dveh ozkih vzporednih tuljav
mi0=4*pi*1e-7;
I=1; % A
R=20e-2;

z=-20e-2:1e-2:20e-2;
 
z1=z-10e-2;
B1=mi0*I*R^2./(2*(z1.^2+R.^2).^1.5)
 
z2=z+10e-2;
 
B2=mi0*I*R^2./(2*(z2.^2+R.^2).^1.5)
 
B=B1+B2;
 
plot(z,B1, z ,B2,z,B,'*-','Linewidth',2)