An Euler-Bernoulli beam and a concentrated mass on this beam are considered as a beam-mass system. The beam is supported by immovable end conditions, thus leading to stretching during the vibrations. This stretching produces cubic non-linearities in the equations. Forcing and damping terms are added into the equations. The dimensionless equations are solved for five different set of boundary conditions. Approximate solutions of the equations are obtained by using the method of multiple scales, a perturbation technique. The first terms of the perturbation series lead to the linear problem. Natural frequencies and mode shapes for the linear problem are calculated exactly for different end conditions. Second order non-linear terms of the pertu...