function Divergenzfeld()
  [x,y] = meshgrid( -6.5:2:6.5 );
  vx = x.^2/2;
  vy = x.*y;
  hold on
  quiver( x,y,vx,vy );         % Darstellung des Vektorfeldes
  div = divergence( x, y, vx, vx );  % Berechnung der Divergenz
  contour( x,y,div )           % Höhenliniendarstellung
  hold off
end