AbstractVariational boundary element formulations are developed. Two functionals based on the dual and complementary energy approach are considered, which provide the upper and lower bound of the solution. Calculated examples of electrostatic fields are demonstrated for their field distributions and capacitances. Validity and the solution accuracy are discussed in comparison with the conventional boundary element solution.