In this note, singular integration problems of the formwhere Ω = [a1, a2] × [b1, b2],x = (x, y) ∈ Ω and fixed point x0 = (x0, y0) ∈ Ω, is considered. The density function h(x, y) is assumed given, continuous and smooth on the rectangle Ω and belong to the class of functions C 2,α (Ω). Cubature formula for double integrals with algebraic singularity on a rectangle is constructed using the modified spline function SΩ(P ) of type (0, 2). Highly accurate numerical results for the proposed method is given for both tested density function h(x, y) as linear, quadratic and absolute value functions. The results are in line with the theoretical findings.