 [x,y] = meshgrid(1:2:20);
 quiver(x,y,sin(2*pi*x/20),sin(2*pi*y/20))
 print('vectorfield.jpg','-djpg')

