Continued Fraction of Squareroot

Continued Fraction can be find in http://mathworld.wolfram.com/ContinuedFraction.html

 

Function ContinuedFraction(ByVal n As Long)Dim root As Long, p As Long, q As Long, b() As String, temp As Long, i As Longroot = Int(Sqr(n))p = rootq = n - p * pIf q = 0 Then ContinuedFraction = "{}": Exit FunctionDo While temp <> 2 * roottemp = (root + p) / qi = i + 1ReDim Preserve b(1 To i)b(i) = tempp = temp * q - pq = (n - p * p) / qLoopContinuedFraction = "{" & Join(b, ",") & "}"End FunctionSub main()Debug.Print ContinuedFraction(99999)End Sub

 

It returns:

{4,2,2,1,2,24,1,13,10,1,1,1,5,3,1,13,3,2,2,16,1,2,7,70,7,2,1,16,2,2,3,13,1,3,5,1,1,1,10,13,1,24,2,1,2,2,4,632}