# File hbody.rb, line 42
  def jerk
    r2 = @pos*@pos
    r3 = r2*sqrt(r2)
    (@vel+@pos*(-3*(@pos*@vel)/r2))*(-@mass/r3)
  end