In this thesis, an approach is presented that is based on the calculation of bifurcations in the spring-mass model. The new approach consists of the transformation of the series of initial value problems on different intervals into a single boundary value problem. Using this technique, discontinuities can be avoided and sophisticated numerical methods for studying parameterized nonlinear boundary value problems can be applied. Thus, appropriate extended systems are used to compute transcritical and period-doubling bifurcation points as well as turning points. We show that the resulting boundary value problems can be solved by the simple shooting method with sufficient accuracy, making the application of the more extensive multiple shooting superfluous. The proposed approach is fast, robust to numerical perturbations and allows to determine highly unstable periodic solutions of the original problem. Asymmetric leg function is often an undesired side-effect in artificial legged systems and may reflect functional deficits or variations in the mechanical construction. It can also be found in legged locomotion in humans and animals, for example after an accident or in specific gait patterns. So far, it is not clear to what extent differences in the leg function of contralateral limbs can be tolerated during walking or running. Here, we address this issue using the spring-mass model for simulating walking and running with compliant legs. With the help of the original realization of the model and the boundary value problem approach, we show that considerable differences between contralateral legs can be tolerated and may even provide advantages to the robustness of the system dynamics. A better understanding of the mechanisms and potential benefits of asymmetric leg operation may help to guide the development of artificial limbs or the design novel therapeutic concepts and rehabilitation strategies.