function [res] =  R(X,Y)
% This function is used by Finite Element Method
% (Algorithm 12.5) to implement R(x,y).
 res = 0;