This project implements a high-performance molecular dynamics simulation engine that compares three numerical integration algorithms—Beeman, Verlet, and Gear Predictor-Corrector—for solving damped and composed oscillator systems. The architecture features a pluggable integrator pattern, allowing seamless switching between integration methods while maintaining sub-millisecond time-step accuracy for systems of up to 1000 coupled particles.

Built for the Simulation Systems course at Buenos Aires Institute of Technology (ITBA).