#!perl -l

# Looking at function F = pi ^ y
#
#             F = pi ^ y
#    error = dF = y * pi ^ (y-1) * dpi
# =>        dpi = error / ( y * x ^ (y-1) )
#
# So if we want to limit the error in F to 0.1:

use constant ERROR => 0.1;
use constant PI    => 4 * atan2(1,1);

warn "Check: pi is ", PI, "\n";

my @y = (1..9, map(10 * $_, 1..9), map(100*$_,1..10));

foreach my $y (@y) {
    printf "y=%d  dpi=%g\n", $y, dpi($y);
}

sub dpi {
    my $y = $_[0];
    my $denom = $y * PI ** ($y-1);
    return ERROR / $denom;
}

__END__
Check: pi is 3.14159265358979
y=1  dpi=0.1
y=2  dpi=0.0159155
y=3  dpi=0.00337737
y=4  dpi=0.000806288
y=5  dpi=0.00020532
y=6  dpi=5.44627e-05
y=7  dpi=1.48594e-05
y=8  dpi=4.13867e-06
y=9  dpi=1.171e-06
y=10  dpi=3.35468e-07
y=20  dpi=1.79111e-12
y=30  dpi=1.27507e-17
y=40  dpi=1.02116e-22
y=50  dpi=8.72341e-28
y=60  dpi=7.76258e-33
y=70  dpi=7.10495e-38
y=80  dpi=6.6385e-43
y=90  dpi=6.30114e-48
y=100  dpi=6.05568e-53
y=200  dpi=5.8364e-103
y=300  dpi=7.5001e-153
y=400  dpi=1.08428e-202
y=500  dpi=1.67203e-252
y=600  dpi=2.68581e-302
y=700  dpi=0
y=800  dpi=0
y=900  dpi=0
y=1000  dpi=0