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: