#include <iostream>
#include <vector>
#include <gsl/gsl_statistics.h>

int main(void) {
  // using plain array
  double data[5] = {17.2, 18.1, 16.5, 18.3, 12.6};
  double mean, variance;
  mean     = gsl_stats_mean(data, 1, 5);
  variance = gsl_stats_variance(data, 1, 5);
  std::cout<<"     mean = "<<mean<<"\n";
  std::cout<<" variance = "<<variance<<"\n";
  // using std::vector
  std::vector<double> v = {17.2, 18.1, 16.5, 18.3, 12.6};
  mean     = gsl_stats_mean(&v[0], 1, 5);
  variance = gsl_stats_variance(&v[0], 1, 5);
  std::cout<<"     mean = "<<mean<<"\n";
  std::cout<<" variance = "<<variance<<"\n";
  }