Added computeNewVelocity in computation.cpp BUT gotta double-check indices (AH)