Using cumtrapz() instead of the equations of motion is the correct method. This way, it can use the change over time, not just an incrementing figure. James Tursa has an excellent display of why this is the correct method. Question answered.