#include <valarray> #include <iostream> double standard_dev(const std::valarray<double> &vals) { return sqrt(pow(vals - (vals.sum() / vals.size()), 2).sum() / vals.size()); } int main() { std::cout << standard_dev({2, 4, 4, 4, 5, 5, 7, 8}) << '\n'; }