The calculation of torsional vibration for marine diesel engine propulsion shafting is normally carried out by equalizing exciting energy to damping energy, or using the dynamic magnifier. But, with these methods, the vibration amplitudes are calculated only for resonance points and vibration amplitudes of other running speeds of engine are determined by the estimation. Recently, many energy-saving ships have been built and on these ships, two-stroke, supercharged, super-long stroke diesel engines which have a small number of cylinders are usually installed. In these cases, the first order critical-torsional vibrations of these engine shaftings appear ordinarily near the MCR speed and the stress amplitudes of their vibration skirts exceed the limit stress defined by the rules of classification society. To predict the above condition in the design stage, the synthesized vibration amplitudes of all orders which are summed up according to their phase angles must be calculated from the drawings of propulsion shaft systems. In this study, a theoretical method to fulfill the above calculation is derived and a computer program is developed according to the derived method. And a shafting system of two-stroke, super-long stroke diesel engine which was installed in a bulk carrier is analyzed with this method. The measured values of this engine shafting are compared with those of calculated results and they show a fairly good agreement.