Faulhaber's formula: Difference between revisions

Content added Content deleted
(add RPL)
(→‎{{header|RPL}}: Improved code using Bernoulli numbers' built-in function)
Line 2,792: Line 2,792:
=={{header|RPL}}==
=={{header|RPL}}==
{{works with|HP|49}}
{{works with|HP|49}}
≪ '''CASE'''
DUP 1 == '''THEN''' 2 / '''END'''
DUP 2 > LASTARG MOD AND '''THEN''' DROP 0 '''END'''
{ }
1 ROT 1 + '''FOR''' m
m INV +
m 1 + → j
≪ '''WHILE''' 'j’ DECR 2 ≥ '''REPEAT'''
j 1 - DUP2 GETI UNROT GET -
OVER * EVAL PUT
'''END'''
'''NEXT'''
HEAD
'''END'''
≫ '<span style="color:blue>BPLUS</span>' STO
≪ → p
≪ → p
≪ 0
≪ 0
p 0 '''FOR''' m
p 0 '''FOR''' m
p 1 + m 1 + COMB
p 1 + m 1 + COMB
p m - <span style="color:blue>BPLUS</span> * EVAL
p m - IBERNOULLI
'n' *
'''IF''' LASTARG 1 == '''THEN''' NEG '''END'''
* EVAL + 'n' *
-1 '''STEP'''
-1 '''STEP'''
p 1 + / EXPAN
p 1 + / EXPAN