The Fibonacci numbers (see http://en.wikipedia.org/wiki/Fibonacci_number) are based on a recurrence relation:

It is difficult to express this relation directly with NumPy code. However, we can express this relation in a matrix form or use the following golden ratio formula:

with

This will introduce the matrix() and rint() functions. The matrix() function creates matrices and the rint() function rounds numbers to the closest integer, but the result is not an integer.