sub std_dev {
  my ($sum, $avg, $variance, $t);

  $sum += $_ for @_;
  $avg = $sum / @_;
  $variance += ($t = ($avg - $_)) * $t for @_;
  $variance /= @_;

  sqrt $variance;
}