This article will combine the finite element method, the interpolated coefficient finite element method, the eigenfunction expansion method, and the search-extension method to obtain the multiple solutions for semilinear elliptic equations. This strategy not only grently reduces the expensive computation, but also is successfully implemented to obtain multiple solutions for a class of semilinear elliptic boundary value problems with non-odd nonlinearity on some convex or nonconvex domains. Numerical solutions illustrated by their graphics for visualization will show the efficiency of the approach.