Abstract: | Coupling techniques are essential to combining different numerical methods together for the purpose of solving an elliptic boundary value problem. By means of nonconforming constraints, the combinations of various Lagrange finite element methods often cause reduced rates of convergence. In this article, we present a method using penalty plus hybrid technique to match different finite element methods such that the optimal convergence rates in the ‖ · ‖h and zero norms of errors of the solution can always be achieved. Also, such a coupling technique will lead to an optimal asymptotic condition number for the associated coefficient matrix. Moreover, this study can easily be extended for combining the finite difference method with the finite element method to also yield the optimal rate of convergence. |