A numerical method is presented to obtain natural frequencies and mode shapes of uniform elastic beams with static deflections due to dead loads. The differential equation governing the free vibration of beam taken into account the static deflection due to deal loads is derived and solved numerically. The hinged-hinged, clamped-clamped and clamped-hinged end constraints are applied in the numerical examples. As the numerical results, the lowest three nondimensional frequency parameters are reported as functions of nondimensional system parameters; the load parameters, and the slenderness rations. And some typical mode shapes of free vibrations are also presented in figures.