To obtain the fundamental solutions for computation of magneto-electro-elastic media by the boundary element method, the general solutions in the case of distinct eigenvalues are derived and expressed in five harmonic functions from the governing equations and the strict differential operator theorem. On the basis of these general solutions, the fundamental solution of infinite magneto-electro-elastic solid are obtained with the method of trial-and-error. Finally, the boundary integral formulation is derived and the corresponding boundary element method program is implemented to perform two numerical calculations(a column under uni-axial tension, uniform electric displacement or uniform magnetic induction, an annular plate simply-supported on outer and inner surfaces under axial loads). The numerical results agree well with the analytical ones.