The present study deals with an approximate integral equation approach to finite deflection of elastic plates with arbitrary plane form. An integral formulation leads to a system of boundary integral equations involving values of deflection, slope, bending moment and transverse shear force along the edge. The basic principles of the development of boundary element technique are reviewed. A computer program for solving for stresses and deflections in a isotropic, homogeneous, linear and elastic bending plate is developed. The fundamental solution of deflection and moment is employed in this program. The deflections and moments are assumed constant within the quadrilateral element. Numerical solutions for sample problems, obtained by the direct boundary element method, are presented and results are compared with known solutions.