分散を求める | Programming Place Plus　Ｃ言語編　逆引き

#include <assert.h>
#include <stdio.h>

#define SIZE_OF_ARRAY(array) (sizeof(array)/sizeof(array[0]))

double sum(const double* array, size_t size)
{
    double result = 0.0;

    for (size_t i = 0; i < size; ++i) {
        result += array[i];
    }
    return result;
}

double average(const double* array, size_t size)
{
    assert(size > 0);
    return sum(array, size) / size;
}

/* 偏差平方和を求める */
double sum_of_squares(const double* array, size_t size)
{
    double result = 0.0;

    // 平均値を求める
    double avg = average(array, size);

    // 各要素の値と平均値との差の２乗を合計する
    for (size_t i = 0; i < size; ++i) {
        result += (array[i] - avg) * (array[i] - avg);
    }

    return result;
}

/* 分散を求める */
double variance(const double* array, size_t size)
{
    assert(size > 0);
    return sum_of_squares(array, size) / size;
}

int main(void)
{
    const double array[] = {5, 5, 2, 14, 4};

    printf("%lf\n", variance(array, SIZE_OF_ARRAY(array)));
}