gold:= proc(n::integer,epsilon ::realcons)
local t,tau:
tau:= (sqrt(5)+1)/2:
if abs (evalf(n/tau – round(n/tau)))<=epsilonthen t:= round(evalf(n/tau))
 else t:=NULL
fi:
[n,t]
end:

number:= proc(n::integer)#epsilon = 0.1
map(i->if nops(gold(i, 0.1)) = 2  then gold(i, 0.1) fi,[seq(k,k=1..n]);
end

Ταπρώτα γινόμεναμήκους 2 με epsilon = 0.1
[5,3], [8,5],[13,8], [21,13], [26,16], [29,18], [34,21], [42,26], [47,29], [50,31]
[55,34], [60,37], [63,39],[68,42], [76,47], [81,50], [84,52], [89,55], [94,58], [97,60]

Αριθμοίτου Fibonacci
with(combinat):
seq(fibonacci(i),i = 1..50);

number f:= proc (n::integer)
map(i->if nops(gold(fibonacci(i),0.1))>=2then gold(fibonacci(i),0.1)fi,
[seq(k,k=1..n)])
end

Λήμμα: Μετο epsilon = 0.1, οι αριθμοίτου Fibonacci εκτόςαπό: 1, 2, 3
             παράγουνγινόμενομήκους 2.

Παρατήρηση:Υπάρχουνπολλαπλάσιατων αριθμώντου Fibonacci πουπαράγουν
                       γινόμενομήκους 2 αλλάόχι όλα:
                       π.χ.  Ναι  26 από το 13,      Όχι  39 από το 13
Επιπλέοντο δεύτερομέρος τουγινομένου ενόςπολλαπλασίουδεν είναιαπαραίτητα τοπολλάπλασιο τουαρχικού:
Ναι:[13,8] και [26,16]   ή   [21,13] και [42,26]
Όχι: [5,3]και [50,31]

Παρατήρηση: Υπάρχουν καιάλλου τύπουαριθμοί πουπαράγουνγινόμενο μήκους2
                        π.χ. 29,47