sub standard_deviation {
  my ($mean, $std_dev) = (0,0);

  $mean += $_ for @_;
  $mean /= @_;

  $std_dev += ($_ - $mean)**2 for @_;
  return $std_dev / @_;
}