sub almost_normal {
    my ($mean, $variance) = @_;
    return sqrt($variance) * atan2( rand(2) - 1, 1 ) + $mean;
}

##</code><code>##

#!/usr/bin/perl -slw

my $PI = 4*atan2(1,1);
sub almost_normal {
    my ($mean, $variance) = @_;
    my $x = rand(2*$PI);
    return sqrt($variance) * sin($x)/cos($x) + $mean;
}

our $PRECISION ||= 1;
our $ITERS = 1000 * $PRECISION;

my %plot;

$plot{ int( $PRECISION * almost_normal( 100, 5 ) ) / $PRECISION }++
    for 1 .. $ITERS;

print "Rand 100 +-5";
printf "%7.2f : %-3d : %s\n", $_, 
    $plot{ $_ }, '#' x ( $plot{ $_ } / $ITERS * 100 * $PRECISION )
    for sort{ $a <=> $b } keys %plot;