Beefy Boxes and Bandwidth Generously Provided by pair Networks
Welcome to the Monastery
 
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 the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others musing on the Monastery: (8)
As of 2019-07-23 13:48 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    If you were the first to set foot on the Moon, what would be your epigram?






    Results (27 votes). Check out past polls.

    Notices?