function A = SetUp(f,xVals,yVals)
% f is a function of the form f(x,y)
% xVals is an array of length Nx
% yVals is an array of length Ny
% A is an Ny-by-Nx matrix with A(i,j) = f(xVals(j),yVals(i))
Nx = length(xVals);
Ny = length(yVals);
A = zeros(Ny,Nx);
for i=1:Ny
for j=1:Nx
A(i,j) = f(xVals(j),yVals(i));
end
end