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