Estimate 
We use the transformation u=1-x, v=1-y,
to transform the singular points to u and v axes, and also the point (1,1); note that the Jacobian is 1.
Thus consider the new integral 
By using the open type quadrature in two dimensions mentioned
earlier, with the uniformly regular matrices, 
we implement the Pascal
program. Partial results are shown below:
