Transverse electromagnetic transmission cells are analyzed with high accuracy using a mixed discretization based on both of constant and linear elements in a boundary element method, in which we can take account into discontinuous values of a normal derivative of the potential at corners in the cells.