steve
/
factor
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

math.extras: fix bernoulli.

db4
John Benediktsson 2013-04-10 09:29:21 -07:00
parent b29d3a0f4a
commit 7d2ad13830
2 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
extra/math/extras/extras-tests.factor
View File

@ -6,6 +6,9 @@ tools.test ;
IN: math.extras.test
{ { 1 -1/2 1/6 0 -1/30 0 1/42 0 -1/30 0 } }
[ 10 iota [ bernoulli ] map ] unit-test
{ -1 } [ -1 7 jacobi ] unit-test
{ 0 } [ 3 3 jacobi ] unit-test
{ -1 } [ 127 703 jacobi ] unit-test

4
extra/math/extras/extras.factor
View File

@ -28,10 +28,10 @@ MEMO: stirling ( n k -- x )
:: ramanujan ( x -- y )
pi sqrt x e / x ^ * x 8 * 4 + x * 1 + x * 1/30 + 1/6 ^ * ;
<PRIVATE
DEFER: bernoulli
<PRIVATE
: (bernoulli) ( p -- n )
[ iota ] [ 1 + ] bi [
0 [ [ nCk ] [ bernoulli * ] bi + ] with reduce