Beefy Boxes and Bandwidth Generously Provided by pair Networks
laziness, impatience, and hubris
 
PerlMonks  

Re: Exponential Function Programming

by fenLisesi (Priest)
on Nov 22, 2007 at 17:21 UTC ( #652421=note: print w/replies, xml ) Need Help??


in reply to Exponential Function Programming

#!/usr/bin/perl -wTl use strict; use warnings; use bignum 'bexp'; ## get the latest from CPAN use Memoize; use constant MAX_POWER => 30; use constant MAX_DIGITS => 50; my $X = shift || 1; memoize('my_factorial'); my $best = bexp( $X, MAX_DIGITS ); my $estimate = 1; for my $i (1 .. MAX_POWER) { $estimate += ($X ** $i) / my_factorial( $i ); print join "\n", ($i, $best, $estimate, abs( $best - $estimate ), q()); } exit( 0 ); ##--------------------------------------------------+ sub my_factorial { my ($n) = @_; return 1 if $n == 1; $n * my_factorial( $n - 1 ); }

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://652421]
help
Chatterbox?
and all is quiet...

How do I use this? | Other CB clients
Other Users?
Others surveying the Monastery: (8)
As of 2016-12-05 09:26 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    On a regular basis, I'm most likely to spy upon:













    Results (75 votes). Check out past polls.