double average(double a[], int count)

{

    double remainder = 0;

    double avg = 0;

    double dcount = count;

    for(int i = 0; i < count; ++i)

    {

        double v, vv;

        avg += v = a[i] / dcount;

        vv = v * dcount;

        remainder += a[i] - vv;

    }

    return avg + remainder / dcount;

}