%%
% vector field
% E(x,y,z)=[y z x]'
%
[x,y,z]=meshgrid(-4:1:4);
u=y;v=z;w=x;
quiver3(x,y,z,u,v,w)

xlabel('{\it x}','FontName','Times','FontSize',16)
ylabel('{\it y}','FontName','Times','FontSize',16)
zlabel('{\it z}','FontName','Times','FontSize',16)
text(-2.5,-5,4.5,'div {\it E} = 0, rot {\it E} = [-1 -1 -1]^T',...
    'FontName','Times','FontSize',14)