function f=fibonaccir ( m )
# prova con m=2000
f = [1 1 ]; % Il vettore f contiene i primi 2 elementi della successione
for n  =  2 : m
    f(n+1)=f(n)+sign(rand-0.5)*f(n-1);
end
semilogy(abs(f),'b-;RFib;')
c=1.13198824;
hold on
semilogy(c.^[1:m],'r-;Exp;')
hold off
endfunction