
%�³“d‰×‚Ç‚¤‚µ

%-----------
% design parameter
Sc=2; % axis scale
K=1;
%-----------
x = -2:0.5:2;
y = -2:0.5:2;
z = -2:0.5:2;

x=Sc*x;
y=Sc*y;
z=Sc*z;

[X,Y,Z] = meshgrid(x, y, z);


U = K*((X+2).^2+Y.^2+Z.^2).^(-3/2).*(X+2) + K*((X-2).^2+Y.^2+Z.^2).^(-3/2).*(X-2);
V = K*((X+2).^2+Y.^2+Z.^2).^(-3/2).*Y + K*((X-2).^2+Y.^2+Z.^2).^(-3/2).*Y;
W = K*((X+2).^2+Y.^2+Z.^2).^(-3/2).*Z + K*((X-2).^2+Y.^2+Z.^2).^(-3/2).*Z;
quiver3(X,Y,Z,U,V,W)
xlabel('x')
ylabel('y')
zlabel('z')
