The function subroutine femodal.m calulates modal parameters for a given structural system. It calculates natural frequency and eigenvector. The eigenvectors are normalized so that the modal mass matrix becomes an indentity matrix.