The function $f(x)=x^2$ is continuous everywhere.
I wrote all this for good measure, but my question is that in most proofs the calculations below are avoided and instead $\delta$ is chosen to be less than $1$ or chosen to be the $\mathrm{min}(1,\frac{\varepsilon}{2|a|+1})$. For example, as done in this question Show Continuity Using Epsilon Delta Definiton . How do we know to choose $\delta < 1$ or $\mathrm{min}(1,\frac{\varepsilon}{2|a|+1})$ as done in the answer of the linked question? Shouldn't the calculations below be necessary in any conclusive proof of the question.
Using the definition of continuity as
$$\forall \epsilon >0\ \forall x_0 >0 \in S \ \exists \delta > 0\ \forall x \in S [\ | x-x_0| < \delta \Rightarrow \ |f(x)-f(x_0)| <\epsilon ]$$, I realize that we have to define our $\delta$ both in terms of $x_0$ and $\epsilon$. So $$| x-x_0| < \delta$$ $$ -\delta +x_0 < x < \delta+x_0$$ $$ -\delta +2x_0 < x + x_0 < \delta+ 2x_0$$ $$-\delta +2|x_0| \le-\delta +2x_0 < x + x_0 < \delta+ 2x_0 \le \delta+ |2x|_0$$
$$-\delta +2|x_0| < x + x_0 < \delta+ |2x_0|$$
$$ |x + x_0| < \delta+ |2x_0|$$
Now since $| x-x_0| < \delta$, it follows $$| x-x_0||x + x_0| < \delta[\delta+ |2x_0|]$$
$$|f(x)-f(x_0)| < \delta^{2}+ 2|x_0|\delta \ .$$
So we let $$\delta^{2}+ 2|x_0|\delta =\epsilon$$ $$\delta^{2}+ 2|x_0|\delta - \epsilon=0$$ and solve using the quadratic formula, that is,
$$\delta= \frac{ (-2|x_0|)\pm \sqrt{ {(2|x_0|)^2}+4\epsilon}}{2} = \frac{-|x_0|\pm \sqrt{ (|x_0|^2)+ \epsilon}}{1}.$$
Since $\delta > 0$ we let $\delta = -|x_0| $ + $ \sqrt{ (|x_0|^2)+ \epsilon}$.
I always feel it is strongly advisable to work all the inequalities out as "rough working" before you attempt to write the formal proof. To show that $$\lim_{x\to a}x^2=a^2$$ you could do the following: $$\eqalign{|x^2-a^2| &=|x-a|\,|x+a|\cr &=|x-a|\,|(x-a)+2a|\cr &\le|x-a|\bigl(|x-a|+2|a|\bigr)\cr}$$ - the reason for these steps is that we want to write the expression in terms of $|x-a|$. Now IMHO the crucial point is that
So $(*)$ I am going to assume $|x-a|<1$, and my inequality implies $$|x^2-a^2|\le|x-a|(1+2|a|)\ .$$ My aim is that the RHS should be less than $\varepsilon$, so I am also going to assume $$|x-a|<\frac{\varepsilon}{1+2|a|}\ .$$ Therefore my choice for $\delta$ will be $$\delta=\min\Bigl(1,\,\frac{\varepsilon}{1+2|a|}\Bigr)\ ,$$ because this guarantees that both the inequalities I need will be true.
Notice, however, that there is really nothing special about the choice $|x-a|<1$. Starting again from $(*)$, I could have assumed $|x-a|<12345$; then through similar working (please try it for yourself) I would have ended up choosing $$\delta=\min\Bigl(12345,\,\frac{\varepsilon}{12345+2|a|}\Bigr)\ .$$ I could even have assumed $|x-a|<\varepsilon$, leading to $$\delta =\min\Bigl(\varepsilon,\,\frac{\varepsilon}{\varepsilon+2|a|}\Bigr)\ .$$
Exercise. Here's one which looks even simpler, but it doesn't work - can you explain why? Assume $|x-a|<|a|$ and hence choose $$\delta=\min\Bigl(|a|,\,\frac{\varepsilon}{3|a|}\Bigr)\ .$$
Note that the above will be rough working, and the actual proof will be as follows.
I'll leave you to fill in the dots. Hope this helps.