function rhs_f = Construct_rhs(n)

%
% Construct grid ponits
%
leng_h = 1 / ( n + 1 );
x      = leng_h:leng_h:n*leng_h;
y      = x;
%
% Compute right hand vector
%
rhs_f  = zeros(n^2,1);
for jj = 1:n
    for ii = 1:n
        rhs_f((jj-1)*n+ii) = 2 * pi^2 * sin(pi * x(ii)) * sin(pi * y(jj));
    end
end

rhs_f = leng_h^2 * rhs_f;