function [u,v] = beanfunc(x,y)

%u = x - y;
%v = x.^3+y;

u = x+(y-y);
v = y+(x-x);

% Jac:
% [ 1    -3y^2 ]
% [ 3x^2     1 ]
%
% det of Jac is 1 + 9x^2y^2.

%u = x - y.^2;
%v = x.^2 + y;
