This paper discusses a harmonic response estimation method on the L$acute{e}$vy plate with two opposite edges simply supported and the other two edges having free boundary conditions. Since the equation of motion of the plate is not self-adjoint, the modes are not orthogonal to each other on the domain. Noting that the L$acute{e}$vy plate can be expressed using one term sinusoidal function that is orthogonal to other sinusoidal functions, this paper suggested the calculation method that is equivalent to finding a least square error minimization solution of the finite number of algebraic equations. Example problems subjected to a distributed area loading and a distributed line loading are defined and their solutions are provided. The solutions are compared to those of the commercial code, ANSYS. According to the verification results, it is expected that the suggested method will be useful to predict the forced response on the L$acute{e}$vy plate with the distributed area or line loading conditions.