|
%I #24 Aug 30 2019 17:48:29
%S 1,2,3,5,7,12,29,41,53,200,253,306,359,665,8286,8951,9616,10281,10946,
%T 11611,12276,12941,13606,14271,14936,15601,31867,79335,111202,190537,
%U 5446238,5636775,5827312,6017849,6208386,6398923,6589460,6779997,6970534,7161071
%N A list of equal temperaments (equal divisions of the octave) whose nearest scale steps are closer and closer approximations to the ratios of two tones of musical harmony: the perfect 4th, 4/3 and its complement the perfect 5th, 3/2.
%C The sequence was found by a computer search of all the equal divisions of the octave from 1 to over 6589460. This is not a perfect recurrent sequence because its self-accumulating nature fails between the 9th and 10th terms, between the 14th and 15th terms, and between the 30th and 31st terms. The examples of recurrence which are present in this sequence are of the same type that is seen in sequences A054540, A060526 and A060527. The numerical value of each term represents a musical scale based on an equal division of the octave. 12, for example, signifies the scale which is formed by dividing the octave into 12 equal parts. - corrected by _K. G. Stier_, Jan 29 2015
%C Also the denominators of increasingly better rational approximations to log(3)/log(2) = 1.5849625... (see A020857). The respective numerators are A254351. The reason why the sequence's "self-accumulating nature fails between the 9th and 10th terms, the 14th and 15th terms and the 30th and 31st terms" (see original comment) is simply that 84/53, 1054/665 and 301994/190537 are very good approximations, thus followed by a jump. (E.g., this phenomenon can also be seen in the numerators and denominators of rational approximations to Pi.). - _K. G. Stier_, Jan 29 2015
%o (Maxima) x:bfloat(log(3)/log(2)),fpprec:100, errold:2,for denominator:1 thru 10000 do (numerator:round(x*denominator), errnew:abs(x-numerator/denominator), if errnew < errold then (errold:errnew, print(denominator))); \\ _K. G. Stier_, Jan 29 2015
%o (PARI) lista(nn) = {d = 2; v = log(3)/log(2); for (den=1, nn, num = round(v*den); newd = abs(v-num/den); if (newd < d, print1(den, ", "); d = newd;););} \\ after Maxima, _Michel Marcus_, Feb 28 2015
%Y Cf. A054540, A060525, A060526, A060527, A254351.
%Y A005664 is a subsequence, A206788 is a supersequence.
%K frac,nonn
%O 1,2
%A Mark William Rankin (MarkRankin95511(AT)Yahoo.com), Apr 12 2001
%E Incorrect term 571611 removed by _K. G. Stier_, Jan 29 2015
%E More terms from _Jon E. Schoenfield_, Feb 06 2015
|
