function b=rhs(k)

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