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