Two numerical methods are used to solve the time independent Schrödinger equation for the torsional potential functions for HONO and HSSH. High accuracy energy levels are obtained with a relatively low computational cost. Theoretical results are compared with the experimental data and different aspects of the dynamic of internal rotation are discussed.