How to compensate drift in solving numerical differential equation

Hi

I work with experimental data from a force platform for human experiment. This platform records forces in the 3 dimensions, Fx, Fy, Fz. Sampling frequency is 50 Hz or 500 Hz.

Thanks to theses forces I would like to calculate velocities of the center of mass of the body and then the mechanical work of theses forces. To do this I must solve two differentials equations successively.

1st to get velocities V : dV/dt=F/M (M is the body mass)

2nd to get work W : dW/dt=F*V

I used ode113 or a personal function using Adams numerical method.

What's the best way to solve such equations, with F being vectors of experimental data and not mathematical function ?

Drift appears in the solutions, what kind of method should be used to avoid it ?

Thanks