By request, I am providing information on how I derived this formula here.
If you list the values of square roots of numbers from $1$ to whatever, the whole numbers place in this list can be seen to follow the sequence $2n+1$, where $n$ is the whole numbers place value of $\sqrt{x}.$ Meaning there are $3$ ones, $5$ two's, $7$ three's,..ect. So, if you wanted to know how many square roots of whole numbers begin with the number $3$, the equation would yield $2*3+1=7$, which is correct. Further, because each first new whole number is also the square root of the perfect square in the beginning of each interval of numbers, $n=\lfloor\sqrt{x}\rfloor$, which is the perfect square before any number $x.$ So, $n^2$ is this perfect square.
It is also intuitive to note that in each interval of values i.e.;(4,8), the values of the decimal places nearly have a uniform distribution with respect to the interval. Using this assumption, its possible to estimate the decimal place values with a fraction. We can represent the location of any number $x$ in an interval using the perfect square before it as, $x-n^{2}.$ Knowing the location of $x$, and the distance in the interval, this fraction can take the form of location over distance as follows,
$\frac{x-n^{2}}{2n+1}.$ This is a beginning approximation to the decimal expansion of $\sqrt{x}.$
Consequently, in knowing that the whole number value will be $n,$ a beginning approximation for $\sqrt{x}$ can be formulated as follows,
$\sqrt{x}\approx\frac{x-n^{2}}{2n+1}+n=\frac{x+n+n^{2}}{2n+1}=\frac{x+n(n+1)}{2n+1}.$
The graph of the error in this, which I will refer to as the "error function", is
$\sqrt{x}-\frac{x+n(n+1)}{2n+1}.$
You can find that the maximum error in the "error function" occurs at about $n(n+1).$ Simply make $x$ equal to $n(n+1).$ Now we have,
$\sqrt{n(n+1)}\approx\frac{2n(n+1)}{2n+1}.$
Now, if you divide the the distance $2n+1$ in each interval by the maximum error given by the error function, you will get a sequence of values that are approaching integers at a rapid rate. This sequence truncated is, ${37,101,197,325,...}$ This is a known sequence whose formula is $4(2n+1)^{2}+1.$ So, if the error divided by $2n+1$ gives almost perfectly $4(2n+1)^{2}+1,$ then of course the error is almost perfectly $\frac{2n+1}{4(2n+1)^{2}+1}.$
Thus, $\sqrt{n(n+1)}\approx\frac{2n(n+1)}{2n+1}+\frac{2n+1}{4(2n+1)^{2}+1}.$
The square root in my post and many other forms I've derived follow from the equation above using simple manipulation. For example, if you evaluate the above equation at $n=n^{2}$ and then divide the whole thing by $n,$ you get the formula in my original post.
That is a lot I know, but this was a long time to derive, research and make perfect. I'm trying to explain it all in one little post. In short, the equation starts from an approximation of the decimal expansion in a square root itself.
Note:In every form above that contained $x$ and $n$ in it together, $n=\lfloor\sqrt{x}\rfloor$ and for the ones containing only $n,$ its just n.