A method to estimate the torsional critical speed for practical rotors has been developed in this study. First, the rotor with a uniform shaft segment is modeled for undamped torsional motion analysis, while satisfying all the boundary conditions. This eventually generates governing equations for the torsional critical speeds of the system. The set of governing equations has the form of a sparse and banded matrix. The elements of banded matrix can be arranged in partitions, which correspond to the specific boundary of the rotor. This permits an automatic generation of the system matrix using a computer. In order to calculate the determinant generated by the simultaneous equations, which leads to the torsional critical speed, a recurring numerical algorithm for a (3*4) sub-matrix has been used. This numerical algorithm practically examines successive (3*4) sub-matrix, one at a time, instead of treating a huge matrix. The output of the program provides the mode shapes with continuous curves. The method has been implemented to three rotors given as examples : a simple rotor, Prohl's rotor, and Macmillan rotor.