We consider eigenvalue problems of the biharmonic operator $\triangle^2$ in a square region $\Omega=(0, 1) \times (0, 1)$ whose boundary condition is that of “free edges”. A finite difference scheme is introduced, and we make numerical computations and compare the results with the physical experiments by Marry D. Waller (1939).