a = 1;
b = 1;
h = 0.04;

function forcer = ff(x,y)
  forcer = ((3*x + x.^2) .* y .* (1 - y) ...
            + (3*y + y.^2) .* x .* (1 - x)) .* exp(x + y);
end

function z = lbc(y)
  z = zeros(size(y))(:);
end

function z = rbc(y)
  z = zeros(size(y))(:);
end

function z = bbc(x)
  z = zeros(size(x))(:);
end

function z = tbc(x)
  z = zeros(size(x))(:);
end

[soln, xs, ys] = poissonSolver(a, b, h, @ff, @lbc, @rbc, @bbc, @tbc);
mesh(xs, ys, soln)
%for azim = -20:10:340
% view(azim, 20)
% pause
%end