Binomial transform: Difference between revisions
add RPL - basic instruction set
(add RPL) |
(add RPL - basic instruction set) |
||
Line 785:
|
≪ DUP SIZE → seq n
≪ { 1 } 2 n '''FOR''' j
'COMB(j-1,k)' 'k' 0 j 1 - 1 SEQ
seq 1 j SUB
* ∑LIST +
'''NEXT'''
≫ ≫ ‘'''FBITR'''’ STO
≪ DUP SIZE → seq n
≪ { 1 } 2 n '''FOR''' j
'(-1)^(j-k-1)*COMB(j-1,k)' 'k' 0 j 1 - 1 SEQ
seq 1 j SUB * ∑LIST +
'''NEXT'''
≫ ≫ ‘'''RBITR'''’ STO
|
Line 814:
|}
If wanting to stick to the basic HP-28 instruction set and no algebraic expressions:
≪ DUP SIZE → seq n
≪ { } 1 n '''FOR''' j
0 1 j '''FOR''' k
seq k GET j 1 - k 1 - COMB * +
'''NEXT''' + '''NEXT'''
≫ ≫ ‘'''FBITR'''’ STO
≪ DUP SIZE → seq n
≪ { } 1 n '''FOR''' j
0 1 j '''FOR''' k
seq k GET j 1 - k 1 - COMB * -1 j k - ^ * +
'''NEXT''' + '''NEXT'''
≫ ≫ ‘'''RBITR'''’ STO
{{in}}
<pre>
| |||