“…For simulation of the dynamics of a constrained membrane, it is necessary to accurately compute many frequencies and mode shapes; this is difficult to compute using existing methods due to singularity. For example, Nagaya [14], Singh and Kothari [16], and Singh et al [17] proposed to combine the two equations for symmetric modes in Eq. (13) into the following single equation…”