This thesis discusses the structure of the spatially beam as the difference conditions between curvature and torsion in the cross-section of the beam ask the change relation of the natural vibration. Every point which follows the centerline of the beam has three translation displacements and three rotational angles of the cross-section area. The Eigenvalue problem is solved the natural frequencies of varying pre-twisted angles of the rectangular cross-section area, slenderness ratio, E/G ratio, all kinds of end conditions. Using the method of Runge-Kutta integration gets the difference equations. Then, adjoint operator system operates them. Finally, the equations of the difference end conditions write the fortran program by controlling seven factors. The datum of the program draw the charts by using Tecplot. The datum discuss the relational change.