login
Orbit periods of the iteration s(1)=0, s(m+1)=(2*i*Pi*n)/log(3)+log(1/(-1-2^-s(m)))/log(3).
0

%I #38 Feb 06 2016 15:12:58

%S 2,1,4,2,1,3,3,1,2,1,6,2,1,3,3,1,2,4,1,2,1,4,5,1,3,3,1,2,1,5,2,1,3,3,

%T 1,2,5,1,2,1,7,5,1,5,4,1,2,1,4,2,1,3,3,1,2,5,1,2,1,3,3,1,5,4,1,2,1,4,

%U 2,1,3,3,1,2,6,1,2,1,3,3,1,2,4,1,2,1,4,2,1,3,3

%N Orbit periods of the iteration s(1)=0, s(m+1)=(2*i*Pi*n)/log(3)+log(1/(-1-2^-s(m)))/log(3).

%C Using the Reduce[1 + 1/2^a + 1/3^b == 0, {a, b}] command in Mathematica to solve the equation 1+1/2^a+1/3^b=0 gives the following conditional expression: b=(2*i*Pi*n)/log(3)+log(1/(-1-2^-a))/log(3), where n is an integer.

%C By setting n=0 and replacing "a" and "b" with "s", setting s(1)=0 as starting value and then iterating:

%C s(1)=0

%C s(2)=(2*i*Pi*n)/log(3)+log(1/(-1-2^(-s(1))))/log(3)

%C s(3)=(2*i*Pi*n)/log(3)+log(1/(-1-2^(-s(2))))/log(3)

%C s(4)=(2*i*Pi*n)/log(3)+log(1/(-1-2^(-s(3))))/log(3)

%C s(5)=(2*i*Pi*n)/log(3)+log(1/(-1-2^(-s(4))))/log(3)

%C ...

%C s(m+1)=(2*i*Pi*n)/log(3)+log(1/(-1-2^-s(m)))/log(3)

%C The iteration converges to sequences with periods of orbit equal to the entries in this sequence a(n). The entries in the converged period are then zeros to equations as in the examples in the conjecture below:

%C a(0)=2, because after 1000 iterations, we have the periodic sequence of repeating convergents with period of the orbit equal to 2:

%C s(m+1) = -0.60158922138 + 2.05252038097 I

%C s(m+2) = -0.60158922138 - 2.05252038097 I

%C s(m+3) = -0.60158922138 + 2.05252038097 I

%C s(m+4) = -0.60158922138 - 2.05252038097 I

%C ...

%C These two complex numbers are then solutions to:

%C Sum_{k=1..3} 1/n^(s(m+1)) + Sum_{k=1..3} 1/n^(s(m+2)) = 0

%C Conjecture:

%C In general for some suitable index value n to a(n), if:

%C a(n)=1 then:

%C Sum_{k=1..3} 1/n^(s(m+1))=0

%C a(n)=2 then:

%C Sum_{k=1..3} 1/n^(s(m+1))

%C +Sum_{k=1..3} 1/n^(s(m+2))=0

%C a(n)=3 then:

%C Sum_{k=1..3} 1/n^(s(m+1))

%C +Sum_{k=1..3} 1/n^(s(m+2))

%C +Sum_{k=1..3} 1/n^(s(m+3))=0

%C ...

%C a(n)=k then:

%C Sum_{k=1..3} 1/n^(s(m+1))

%C +Sum_{k=1..3} 1/n^(s(m+2))

%C +Sum_{k=1..3} 1/n^(s(m+3))

%C +...

%C +Sum_{k=1..3} 1/n^(s(m+k))=0

%C and so on, by observation.

%C The index n of positions of ones in the sequence a(n) gives a linear looking plot. The sequences of positions of ones starts: 1, 4, 7, 9, 12, 15, 18, 20, 23, 26, 28, 31, 34, 37,... and is equal to the first 34 entries of A184589, but then begins to differ at n=35. The differences between the index positions of ones starts: 3, 3, 2, 3, 3, 3, 2, 3, 3, 2, 3, 3, 3,.. which appears to have initial periodic patterns that fail later on, and each time the pattern fails the periodic patterns get longer and longer. The values of n when a(n)=1, the function in the equation the iteration solves is the partial Riemann zeta function, here 1 + 1/2^s+ 1/3^s = 0.

%C The Mathematica program can treat all partial Riemann zeta functions but caution is needed when changing the parameters of the program in order for the formula to converge.

%t (*Mathematica 8*)

%t Clear[k, m, s, kk, n, aa, b1, t, aaa, kkk, aaa]

%t Print["k can be varied to any integer greater than or equal to 2:"]

%t k = 3;

%t m = 1000;

%t s = 0;

%t Monitor[aaa =

%t Table[Table[

%t Do[s = (2 I \[Pi]*(kk))/Log[k] +

%t N[Round[Log[-1/(Sum[1/n^s, {n, 1, k - 1}])]/Log[k]*10^20]/

%t 10^20, 20], {i, 1, m}];

%t aa = Table[

%t s = (2 I \[Pi]*(kk))/Log[k] +

%t N[Round[Log[-1/(Sum[1/n^s, {n, 1, k - 1}])]/Log[k]*10^20]/

%t 10^20, 20], {i, 1, 100}];, {kk, kkk, kkk}];

%t (*end*)

%t Flatten[

%t Position[Chop[Accumulate[Sum[1/n^(aa), {n, 1, k}]]],

%t 0]][[1]], {kkk, 0, 90}], kkk]

%t b1 = Flatten[Position[aaa, 1]] - 1;

%t Differences[b1];

%K nonn

%O 0,1

%A _Mats Granvik_, Jan 30 2016