%%% vfd.m, script for vectorfields in plane
% rha, 26 feb 2004 on vaio
% cf cooper, p. 220
% this exm is (1, x+y^2)
u = inline('0*x +1', 'x', 'y')
v = inline('x + y.^2', 'x', 'y')
x = linspace(-2,3,11);
y = linspace(-1,2,11);
[X,Y] = meshgrid(x,y);
U = u(X,Y); V = v(X,Y);
quiver(X,Y,U,V)
axis image