In this study, wave reflection due to depression of seabed is calculated by using eigenfunction expansion method. The proper numbers of steps and evanescent modes needed for analysis are suggested by applying the eigenfunction expansion method to bottom topography of which slope or curvature varies. While satisfying shallow or intermediate water depth condition, the optimal figure of depression of seabed is obtained by calculating reflection coefficient for various depressions of seabed. The reflection coefficient with distance between the depression of seabeds is then calculated after arraying the optimal geometry in two and three rows.