In this paper, the governing equations for a vibratory beam with moderately large deflection are derived using the first order shear deformation theory. These equations which are a system of nonlinear partial differential equations with constant coefficients are solved analytically with the perturbation technique and the natural frequencies and the buckling load of the system are determined. A parametric study is performed and the effects of the geometrical and material properties on the natural frequency and buckling load are investigated and the effect of normal transverse strain and axial load on natural frequency are examined. Some results based on the first order shear deformation theory are consistent with classic theories of beams and some yield different results. Formulation presented to calculate the transverse frequency, determines the axial frequency too. Also, the natural frequencies and buckling load are calculated with the finite elements method by applying one and three-dimensional elements and the results are compared with the analytical solution.