Remove commented code

[dataframing]

See:

#(dx, dy, dz) = compute_deltas(x1, x2, y1, y2, z1, z2)
(dx, dy, dz) = (x1-x2, y1-y2, z1-z2)
#e -= compute_energy(m1, m2, dx, dy, dz)

and

#update_rs(r, dt, vx, vy, vz)

[nyuxz]

thanks.