We consider the shape optimization for the clamped and buckling plate eigenvalue problems. We develop a numerical method using Hadamard's shape derivative and the method of fundamental solutions as forward solver which allows one to propose numerical candidates to be minimizers for the first ten eigenvalues of both problems.