The purpose of this basic C integrator is to find the numerical position-momentum solution of a Newtonian particle given arbitrary initial conditions. It implements 4 basic techniques (forwards Euler, backwards Euler, trapezoid, and leapfrog).