sub std_dev {
  my ($sum, $sqr_sum);
  for (@_) {
    $sum += $_;
    $sqr_sum += $_ * $_;
  }
  sqrt(($sqr_sum - $sum * $sum / @_)/@_);
}