The dynamical equations governing pulsation in rotating neutron stars are derived in the framework of general relativity. Stellar models are constructed by using a realistic equation of state for cold neutron matter. Small radial displacement and slow rotation are treated as perturbations on spherically symmetric body. In these models the maximum masses are 1.761 M⊙ at the central density 3.461 × 1015 g cm−3 for a sequence of nonrotating configurations and 2.165 M⊙ for rotating models with the critical angular velocity (GM/R3)½.