It might be interesting that the recursion
$$a(n+1) = \frac{1}{2} (a(n) + \frac{c}{a(n)})\tag{1}$$
can be solved explicitly.
The solution is
$$a(n) = \sqrt{c}\; \coth {(2^{n-1}\;\text{arccoth}( a(1)/\sqrt{c}))}\tag{2} $$
For large $n$ the argument of the $\text{coth}$ gets large as well so that the $\text{coth} \to 1$ and $a(n) \to \sqrt{c}$, as expected.
Derivation
First we normalize the problem letting
$$d(n) = a(n) \sqrt{c}$$
which leads to the recursion
$$d(n+1) = \frac{1}{2} (d(n) + \frac{1}{d(n)})\tag{3}$$
The sum of $d(n)$ and its inverse suggests the ansatz
$$d(n) = \coth(f(n)) = \frac { e^{f(n)} + e^{- f(n)} } { {e^{ f(n)} - e^{- f(n)} }}\tag{4}$$
Inserting this into the r.h.s of (3) gives
$$\frac{ e^{2 f(n)} + e^{-2 f(n)} } { {e^{2 f(n)} - e^{-2 f(n)} }}$$
which has the same form as $d(n+1)$ if we let $f(n)$ obey the recursion
$$f(n+1) = 2 f(n)\tag{5}$$
But this can be solved immediately
$$f(n) = 2^{n-1} f(1) $$
so that the solution to (3) is given by
$$d(n) = \coth {(2^{n-1}f(1))}$$
The constant $f(1)$ has to be adapted to the initial condition leading to
$$f(1) = \text{arccoth}( d(1)) $$
Hence the solution to (3) becomes
$$d(n) = \coth {(2^{n-1}\text{arccoth}( d(1)))} $$
And, reversing the nomalization, we find the solution (2).