Fibonacci numbers (Forth)

From LiteratePrograms
Jump to: navigation, search
Other implementations: bc | C | C Plus Plus templates | dc | E | Erlang | Forth | FORTRAN | Haskell | Icon | Java | JavaScript | Lisp | Lua | occam | Oz | PIR | Prolog | Python | Sed | sh | sh, iterative | T-SQL | Visual Basic .NET

The Fibonacci numbers are the integer sequence 0, 1, 1, 2, 3, 5, 8, 13, 21, ..., in which each item is formed by adding the previous two. The sequence can be defined recursively by


 F(n) =
 \begin{cases}
   0             & n = 0 \\
   1             & n = 1 \\
   F(n-1)+F(n-2) & n > 1 \\
  \end{cases} .

Fibonacci number programs that implement this definition directly are often used as introductory examples of recursion. However, many other algorithms for calculating (or making use of) Fibonacci numbers also exist.


Contents

[edit] Implementation

The Fibonacci numbers in Forth:

<<fib.fs>>=
fib-iter
fib-rec
main

Note that even these implementations are only suitable for small values of n, since the Fibonacci function grows exponentially and even on a 64-bit machine unsigned Forth cells can only hold the first 92 Fibonacci numbers.

[edit] Iteration

This is the easiest way to calculate fibonacci numbers. The values of fibi-2 and fibi-1 are stored on the stack, to avoid all the recalculations in the recursive implementation.

( n -- f ) is a stack comment documenting the inputs (n) and outputs (fibn) of the function. It is only a comment, but this style is a de facto standard.

0 1 rot ( 0 1 n ) sets up the initial terms of the sequence and rotates the loop count to the top of the stack.

0 ?do ... loop ( n -- ) takes the count from the stack leaving the two previous Fibonacci terms on the stack for the body of the loop. ?do does an initial check against the limit compared to do.

over ( a b -- a b a ) copies the item under the top of the stack to the top. dup ( a -- a a ) duplicates the top item on the stack.

+ ( a b -- a+b ) adds the top two items on the stack, in this case fibn-1 and fibn-2, yielding fibn.

swap ( a b -- b a ) swaps the top two items on the stack.

drop ( a -- ) gets rid of the extra term before returning from the function.

<<fib-iter>>=
: fib-iter ( n -- f )
  0 1 rot 0 ?do over + swap loop drop ;

[edit] Recursive

This is a recursive implementation. This will become slow on big numbers, because the numbers are recalculated for each recursion.

dup 2 u< if exit then returns n if it is 0 or 1. u< is an unsigned comparison compared to <.

1- ( n -- n-1 ) is the decrement operator.

recurse causes the current function to recurse.

<<fib-rec>>=
: fib-rec ( n -- f )
  dup 2 u< if exit then
  1- dup recurse  swap 1- recurse  + ;

[edit] Test driver

This test driver treats the two Fibonacci functions anonymously. ' fib-iter returns the execution token or xt of the word, and execute performs its action. An xt is similar to a C language function pointer.

<<main>>=
: test ( xt count -- )
  0 do
    i over execute .
  loop drop ;

cr .( fib-iter 0..9: )  ' fib-iter 10 test
cr .( fib-rec  0..9: )  ' fib-rec  10 test
Download code
hijacker
hijacker
hijacker
hijacker