The purpose of this paper is to apply the coupling of mixed finite element and boundary integral methods, using two integral equations on the coupling boundary, to study the weak solvability and approximation of certain class of linear and nonlinear elliptic boundary value problems arising in potential theory and elastostatics.