Bernoulli numbers: Difference between revisions

m
→‎{{header|Sidef}}: minor code simplications
(→‎{{header|Sidef}}: code simplications and added 2 more algorithms)
m (→‎{{header|Sidef}}: minor code simplications)
Line 3,801:
var K = 2*(n! / Num.tau**n)
var d = n.divisors.grep {|k| is_prime(k+1) }.prod {|k| k+1 }
var Nz = ceil((K*d).root(n-1)).primes.prod {|p| 1 - p.float**(-n) }
 
var z = 1
for p in (primes(N)) {
z *= (1 - float(p)**(-n))
}
 
(-1)**(n/2 + 1) * int(ceil(d*K / z)) / d
Line 3,814 ⟶ 3,809:
<lang ruby>func bernoulli_print {
var a = []
for m in (0..60) {
a.append( << 1/(m+1))
for j (flipin (1..m) -> flip) {
(a[j-1] -= a[j]) *= j
}
2,747

edits