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

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

u = x - 0.2*y.^3;
v = 0.2*x.^3 + y;

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

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