function [res] = GG1(X,Y,A,B,C,I,I1,I2)
% This program is used by Finite Element Method
% (Algorithm 12.5).
 res = G1(X,Y)*(A(I,I1)+B(I,I1)*X+C(I,I1)*Y)*(A(I,I2)+B(I,I2)*X+C(I,I2)*Y);