In classical higher-order discontinuous boundary element formulation for two-dimensional elastostatics, interpolation functions for different boundary variables (i.e., boundary displacements and tractions) are assumed to be the same. However, there is a derivational relationship between these variables. This paper presents a boundary element formulation, called Mixed Boundary Element Formulation, for two dimensional elastostatic problems in which above mentioned relationship is taking into account. The formulations are performed by using discontinuous first and second-order mixed boundary elements. Based on the formulations presented in this study, two computer softwares are developed and verified through some example problems. The results show that the present formulation is credible.