def calc_velocity(vec)
      v = vec.inject(0.0) { |s, i| s + i * i }
      return Math.sqrt(v)
    rescue => e
      raise
    end
def calc_velocity(vec)
      v = vec.inject(0.0) { |s, i| s + i * i }
      return Math.sqrt(v)
    rescue => e
      raise
    end
=========================================================================
 天体の速度ベクトルから実際の速度を計算

 @param:   vec  (ベクトル)
 @return:  v    (速度)
=========================================================================