Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #17 Aug 27 2019 17:49:52
%S 1,1,1,2,2,2,2,3,2,2,4,3,2,3,5,3,3,4,5,3,3,5,5,3,4,6,5,4,5,7,5,4,6,7,
%T 4,5,7,7,4,5,8,6,5,6,8,6,5,7,9,6,5,8,9,6,6,9,8,5,6,10,8,6,7,10,8,6,8,
%U 10,7,6,9,10,7,6,10,10,7,7,11,10,7,7,11
%N a(n) = nearest integer to f(n), where f(0) = f(1) = 1, f(n) = (n-f(n-1))/f(n-2).
%C Comments from _Jon Maiga_, Aug 25 2019: (Start)
%C The function seems to follow the square root of n.
%C The first differences center around the x-axis with a similar look.
%C Changing the initial values to for example f(0)=1 and f(1)=3 creates a very different graph.
%C See the image for fractional versions with f(0)=f(1)=1 and f(0)=1, f(1)=3.
%C (End)
%H Robert Israel, <a href="/A308622/b308622.txt">Table of n, a(n) for n = 0..10000</a>
%H J. Maiga, <a href="http://jonkagstrom.com/look/index.html">A few good looking recursive functions</a>, 2019.
%H Jon Maiga, <a href="/A308622/a308622.png">Plots of f(n) for n up to 1000</a>
%p L[0]:= 1: R[0]:= 1: A[0]:= 1:
%p L[1]:= 1: R[1]:= 1: A[1]:= 1:
%p for n from 2 to 100 do
%p x:= (n - R[n-1])/R[n-2];
%p y:= (n - L[n-1])/L[n-2];
%p L[n]:= 10^(-100)*floor(x*10^100);
%p R[n]:= 10^(-100)*ceil(y*10^100);
%p if round(L[n]) <> round(R[n]) then printf("Oops: %d \n",n); break fi;
%p A[n]:= round(L[n])
%p od:
%p seq(A[n],n=0..100); # _Robert Israel_, Aug 24 2019
%t f[0] = f[1] = 1;
%t f[n_] := f[n] = N[(n - f[n - 1])/f[n - 2]]
%t Round[Array[f, 83, 0]]
%K nonn,look
%O 0,4
%A _Jon Maiga_, Jun 11 2019