This thesis contains the development of continuous Kepler orbit- and a discrete numerical integration-based collision detection algorithms in a system of LEO satellites, which in combination with collision algorithm form a simplified space debris evolution model. This model is then used to study the Kessler syndrome. The continuous and discrete algorithms get their names from the solutions of the Two Body Problem (TBP) and the methods for collision detection that they are based on; the analytical and continuous time solution of TBP resulting in the Kepler orbits and the numerical, discrete time Velocity Verlet integration of the TBP. The collision model consists of an algorithm for fragmentation collisions largely based on the NASA Standard Breakup Model and a method for elastic, random scattering collisions. Comparison between the continuous and discrete algorithms shows that on average both predict the same time to the first collision in a system of homogeneously distributed satellites. The algorithms differ in their efficiency depending on the number and the radius of the satellites in and the geometry of the system. For relatively small satellite numbers in large systems, the continuous algorithm is computationally more efficient. However, as more satellites or fragments result from previous collision, the continuous algorithm is outperformed by the discrete algorithm. Consequentially, its time complexity appears to be O(N²). Armed with this knowledge, the continuous algorithm is used to show that an initially small system of satellites is able to evolve into a large population of debris particles within several decades. Similarly, the discrete algorithm is used to show that an ordered collection of satellites in an homogeneously distributed system of debris-like particles exhibits the effect that a collision early on in the simulation can cause a cascade of collisions at a later stage. Hence Both the discrete and continuous algorithms predict a Kessler Syndrome and mimic predictions made by more advanced models from leading space agencies like NASA’s LEGEND, ESA’s DELTA and JAXA’s LEODEEM [Lio+13].Future research could focus on including atmospheric drag and gravitational perturbations to the continuous algorithm, thereby lengthening the time frame during which it can realistically simulate a system of satellites in LEO. To achieve this, it is suggested that one execute the calculations inherent to the algorithm in parallel on a GPU, as these are independent of each other.