Составить процедуру вычисления среднего значения элементов массива

#include <iostream>
#include <vector>
#include <algorithm>
#include <iterator>
#include <time.h>
#include <cstdlib>
using namespace std;
 
/*объект функции для вычисления среднего арифметического*/
class Average{
    private:
        double average;
        int number;
 
    public:
        Average():average(0), number(0) {}
 
        void operator() (int elem){
            average += elem;
            ++number;
        }
 
        operator double(){
            return static_cast<double>(average) / number;
        }
};
 
/*функция рандомного заполнения*/
void random(double &elem){
    elem = rand() % 21;
}
 
void number_element(vector<double> &coll){
 
    cout << "Array: ";
    copy(coll.begin(), coll.end(), ostream_iterator<double>(cout, " "));
 
    /*считаем среднее значение*/
    double average = for_each(coll.begin(), coll.end(), Average());
 
    cout << "\nAverage: " << average;
 
    /*определяем кол-во элементов больших срденого*/
    int number = count_if(coll.begin(), coll.end(), 
        bind2nd(greater<double>(), average));
 
    cout << "\nThe number of elements greater than the average: "
        << number << "\n\n";
}
 
int main()
{
    srand(time(NULL));
 
    int size(10);
    vector<double> coll1(size), coll2(size);
 
    /*заполняем массивы рандомно*/
    for_each(coll1.begin(), coll1.end(), random);
    for_each(coll2.begin(), coll2.end(), random);
 
    number_element(coll1);
    number_element(coll2);
 
    cout << "\n\n";
    return 0;
}