double medianOfVector(vector<double> &vec) {
    double result = 0;
    //    vector<double> vec2;
    if ((vec.size() % 2) != 0) {
        //        vec2 = vec;
        vector<double>::iterator i = vec.begin();
        vector<double>::size_type m = (vec.size() / 2);

        nth_element(i, i + m, vec.end());
        result = vec.at(m);
    } else {
        vector<double>::iterator i = vec.begin();
        vector<double>::size_type m = (int) (((vec.size() - 1) / 2));
        nth_element(i, i + m, vec.end());
        double min = vec.at(m);
        double max = *min_element(i + m + 1, vec.end());
        result = (min + max) / 2;
    }

    return result;
}