
disp(' ')
disp('     Computing the flux through the boundary   ')
disp('  We use the vector field ' )
disp('      F(x,y)  = (u(x,y),v(x,y)) = ( y, 2xy) ')
disp('      div(F) = 2x  ')

corners = [-2 2 -2 2];
u = inline(' y', 'x', 'y');
v = inline(' 2*x.*y', 'x', 'y');

flux2(u,v,corners) 
pause 
close
