login
This site is supported by donations to The OEIS Foundation.

 

Logo

Many excellent designs for a new banner were submitted. We will use the best of them in rotation.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A072895 Least k for the Theodorus spiral to complete n revolutions. 6
17, 54, 110, 186, 281, 396, 532, 686, 861, 1055, 1269, 1503, 1757, 2030, 2323, 2636, 2968, 3320, 3692, 4084, 4495, 4927, 5377, 5848, 6338, 6849, 7378, 7928, 8497, 9087, 9695, 10324, 10972, 11640, 12328, 13036, 13763, 14510, 15277, 16063, 16869 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

"For n = 16 (which gives sqrt(17)) this sum is 351.15 degree, while for n = 17 the sum is 364.78 degree. That is, perhaps Theodorus stopped at sqrt(17) simply because for n > 16 his spiral started to overlap itself and the drawing became 'messy.'" Nahin p. 34.

There exists a constant c = 1.07889149832972311... such that b(n) = IntegerPart[(Pi*n + c)^2 - 1/6] differs at most by 1 from a(n) for all n>=1. At least for n<=4000 we indeed have a(n)=b(n). - Herbert Kociemba, Sep 12 2005

The preceding constant and function b(n) = a(n) for all n < 21001. - Robert G. Wilson v, Mar 07 2013; Update: b(n) = a(n) for all n < 10^9. - Herbert Kociemba, Jul 15 2013

The preceding constant, c, is actually is -K/2, where K is the Hlawka's Schneckenkonstante (A105459). - Robert G. Wilson v, Jul 10 2013

REFERENCES

Julian Havil, The Irrationals, A Story of the Numbers You Can't Count On, Princeton University Press, Princeton and Oxford, 2012, page 35.

Paul J. Nahin, "An Imaginary Tale, The Story of [Sqrt(-1)]," Princeton University Press, Princeton, NJ. 1998, pgs 33-34.

David Wells, "The Penguin Dictionary of Curious and Interesting Numbers," Revised Edition, London, England, 1997, page 76.

LINKS

Robert G. Wilson v, Table of n, a(n) for n = 1..21000

Herbert Kociemba  The Spiral of Theodorus

FORMULA

a(n) = k means k is the least integer such that Sum_( atan( 1/sqrt(i)), i=1..k) > 2n*Pi.

a(n) = A137515(n) + 1. [Robert G. Wilson v, Feb 27 2013]

MATHEMATICA

s = 0; k = 1; lst = {}; Do[ While[s < (2Pi)n, (* change the value in the parentheses to change the angle *) s = N[s + ArcTan[1/Sqrt@k], 32]; k++]; AppendTo[lst, k - 1], {n, 50}]; lst (* Robert G. Wilson v, Oct 14 2012 *)

K = -2.15778299665944622; f[n_] := Floor[(n*Pi - K/2)^2 - 1/6]; Array[f, 41] (* Robert G. Wilson v, Jul 10 2013 *)

K = -2.1577829966594462209291427868295777235; a[n_] := Module[{a = -(K/2) + n Pi, b}, b = a^2 - 1/6; If[Floor[b] == Floor[b + 1/(144 a^2)], Floor[b], Undefined]] (* defined at least for all n < 10^9, Herbert Kociemba, Jul 15 2013 *)

CROSSREFS

Cf. A002194, A105459.

Sequence in context: A146405 A228244 A158968 * A097059 A117390 A141841

Adjacent sequences:  A072892 A072893 A072894 * A072896 A072897 A072898

KEYWORD

nonn

AUTHOR

Robert G. Wilson v, Jul 29 2002

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified April 19 19:16 EDT 2014. Contains 240777 sequences.