The modified quadrature formulas with singularities of f over quadrilateral elements are discussed and modified integration points and weights are derived for these elements. The accuracy of this modified quadrature scheme is checked over a flat square element by evaluating the integrals of the equations of the boundary element methods and results are compared with those using Gauss-quadrature rule over this element.