Bernoulli numbers: Difference between revisions

Line 886:
 
def bernoulli2(limit)
m, ar = 0, [] of BigRational
(0..limit).each do |m|
while m < limit
ar << BigRational.new(1, m+1)
m.downto(1) { |j| ar[j-1] = j * (ar[j-1] - ar[j]) }
yield ar[0] # use value in required block
m += 1
end
end
Anonymous user