x=linspace(-1,1,10);
y=linspace(-1,1,10);
z=linspace(-1,1,10);
[X,Y,Z]=meshgrid(x,y,z);
quiver3(X,Y,Z,X.^2,Y.^2+X.^2,Z) 