Programming has been carried out to model the vibration of a mass spring system with no damping and with damping variations using Python programming language coding. This research aims to simulate the simple harmonic vibration of a mass spring (without damping) and the vibration of a damped mass spring. The programming is designed after formulating the equations of motion of damped mass spring vibrations that behave as second-order differential equations and analyzing the numerical formulation of the Feynman-Newton algorithm arrangement. The method used is experimentation using python to simulate the vibration of the spring. Simulation by varying the damping constant (c = 0, c < √4mk, c = √4mk, and c > √4mk). The simulation results show various vibration graphs, namely simple harmonic vibration, under damped vibration, critically damped vibration, and over damped vibration. The shape of the vibration graph is influenced by the mass, spring constant, and damping constant. The greater the damping constant, the less the maximum speed of vibration. This research has succeeded in visualizing the simple harmonic vibration of a massless spring (without damping) and the vibration of a damped massless spring. This modeling can help the physics learning process in high school in understanding the concept of spring vibration to obtain critical reasoning in accordance with physical phenomena. Modeling of damped spring vibrations using Python can be used as a physics learning media on vibration material to train scientific and critical reasoning skills, and as a learning innovation because it is a new thing where the curriculum and high school learning in Indonesia are not used to being delivered.