%%% another 3D vectorfield plot
% rha on vaio 28 feb 2004
% cf cooper p. 221 but no loop
u = inline('1 + 0*x', 'x', 'y', 'z')
v = inline('x + y.^2', 'x', 'y', 'z')
w = inline('z', 'x', 'y', 'z')
x = linspace(-2,3,6);
y = linspace(-1,2,6);
z = linspace(-1,1,4);
[X,Y,Z] = meshgrid(x,y,z);
U = u(X,Y,Z);
V = v(X,Y,Z);
W = w(X,Y,Z);
quiver3(X,Y,Z,U,V,W)
