Investigated herein is the buckling of .nanobeams based on a nonlocal Timoshenko beam model by the method of initial values within the framework of nonlocal elasticity. Since the nonlocal Timoshenko beam theory is of higher order than the nonlocal Euler-Bernoulli beam theory, it is known to be superior in predicting the small-scale effect. The buckling determinants and critical loads for bars with various kinds of supports are presented. The Carry-Over matrix (Transverse Matrix) is presented and the priorities of the method of initial values are depicted. To the best of the researchers' knowledge, this is the first work that investigates the buckling of nonlocal Timoshenko beam with the method of initial values.