Algorithm and code are presented which solve the dispersion equation for cylindrical poroelastic structures. The algorithm is based on the spectral method, which discretizes the underlying wave equations with the help of spectral differentiation matrices and solves the corresponding equations as a generalized eigenvalue problem. The results are illustrated for the case of a fluid-saturated free cylinder with open- and closed-pore boundary conditions on its surface. The computed dispersion curves are in good agreement with analytical results, which confirms the accuracy of the method.