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

 

Logo

Annual Appeal: Today, Nov 11 2014, is the 4th anniversary of the launch of the new OEIS web site. 70,000 sequences have been added in these four years, all edited by volunteers. Please make a donation (tax deductible in the US) to help keep the OEIS running.

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 November 20 17:42 EST 2014. Contains 249753 sequences.