function b=rhs(n)

% a simple rhs for testing the iterative solvers
b = ones((n-2)^2,1);
