[x, y, z] = meshgrid([-1.5:0.1:+1.5]);
u = (-1).*((y+1)./(x.^2+(y+1).^2)-(y-1)./(x.^2+(y-1).^2)+z*0);%%Bx
v = x.*(1./(x.^2+(y+1).^2)-1./(x.^2+(y-1).^2)+z*0);%%By
w = 0;%%Bz
figure;
hold on
grid on
quiver3(x, y, z*0, u, v,z*0,1);%%campo
title('hilo infinito con corriente');
xlabel('Bx');
ylabel('By');
zlabel('Bz');
t = -2:pi/50:2;
plot3(t*0,t*0+1,t,'r');
plot3(t*0,t*0-1,t,'g');
%%la esfera es por defecto, no tiene abssolutamente nada 
%%que ver con el ejemplo
%%solo para que funcione rapido el grafico
r = 0.001;
[a,b,c] = sphere();
surf(a*r,b*r,c*r);