Abstract: Based on the fundamental equations of the plane magnetoelectroelastic solids and the basic idea of virtual boundary element method for elasticity, a virtual boundary element-least square collocation method (VBEM) for plane magnetoelectroelastic solids is presented. Besides all the advantages of the conventional boundary element method (BEM) over domain discretization methods, the method avoids the computation of singular integral on the boundary by introducing the virtual boundary. This method, merely using collocations technology on the real and virtual boundaries, is mesh-free and integration-free. In the end, several numerical examples are presented to demonstrate the performance of the proposed method. The results show that they agree well with the exact solutions.