The purpose of this thesis is to find effective algorithms to numerically solve certain systems of differential equations that arise from standard Newtonian mechanics. Numerical models of elastica has already been well studied. In this thesis we concentrate on the Kirchhoff problem. The goal is to create an effective and robust numerical method to model the dynamic behavior of springs that have a prescribed natural curvature. In addition to the mathematics, we also provide the implementation details of the numerical method using the computer language Python 3. We also discuss in detail the various difficulties of such a software implementation and how certain auxiliary computations can make the software more effective and robust.