When you sum the terms in the diagonal you don't get
$$\sum_{i=1}^n |a_i|^2\sum_{j=1}^n |b_j|^2.$$
Instead, what you get is
$$\sum_{i=1}^n |a_ib_i|^2.$$
Let's write $[n] = \{1,\ldots,n\}$. Here it helps to separate the sum in two sums, one where the indexes agree and one where the indexes are different. Notice that
$$
\sum_{i=1}^n |a_i|^2\sum_{j=1}^n |b_j|^2
= \sum_{i,j\in[n]} |a_ib_j|^2
= \sum_{i=j} |a_ib_j|^2 + \sum_{i\neq j} |a_ib_j|^2
= \sum_{i=1}^n |a_ib_i|^2 + \sum_{i\neq j} |a_ib_j|^2.
$$
But if $i\neq j$, of course $j\neq i$, so
$$\sum_{i\neq j} |a_ib_j|^2 = \sum_{1\leq i\lt j\leq n} |a_ib_j|^2 + |a_jb_i|^2.$$
Following your line of thought
\begin{align}
\left|\sum_{i = 1}^na_ib_i\right|^2
&=\left(\sum_{i = 1}^na_ib_i\right)\left(\sum_{j = 1}^n\bar{a}_j\bar{b}_j\right)\\
&= \sum_{i,j\in[n]}a_i\bar{a}_jb_i\bar{b}_j\\
&= \sum_{i=1}^n |a_ib_i|^2 + \sum_{1\leq i\lt j\leq n}a_i\bar{a}_jb_i\bar{b}_j+a_j\bar{a}_ib_j\bar{b}_i \\
&= \sum_{i=1}^n |a_ib_i|^2 + \sum_{i\neq j} |a_ib_j|^2 - \sum_{i\neq j} |a_ib_j|^2 + \sum_{1\leq i\lt j\leq n}a_i\bar{a}_jb_i\bar{b}_j+a_j\bar{a}_ib_j\bar{b}_i \\
&= \sum_{i=1}^n |a_i|^2\sum_{j=1}^n |b_j|^2 -
\left(
\sum_{i\neq j} |a_ib_j|^2 - \sum_{1\leq i\lt j\leq n}a_i\bar{a}_jb_i\bar{b}_j+a_j\bar{a}_ib_j\bar{b}_i
\right) \\
&= \sum_{i=1}^n |a_i|^2\sum_{j=1}^n |b_j|^2 -
\sum_{1\leq i\lt j\leq n} |a_ib_j|^2 + |a_jb_i|^2 - a_i\bar{a}_jb_i\bar{b}_j-a_j\bar{a}_ib_j\bar{b}_i \\
&= \sum_{i=1}^n |a_i|^2\sum_{j=1}^n |b_j|^2-\sum_{1\leq i\lt j\leq n}\lvert a_i\bar{b}_j - a_j\bar{b}_i\rvert^2
\end{align}
You can start with the parametrization
$$ x(s) = a \cos s $$
$$ y(s) = b \sin s $$
and write
$z(s) = x(s) + i y(s) = a (\cos s + i \sin s) + (b-a)(i \sin s)$
$ = a e^{i s} + \frac{b-a}{2}( e^{is}-e^{-is}) $
$ = (a + \frac{b-a}{2}) e^{i s} - (\frac{b-a}{2}) e^{-is}. $
$ = A_+ e^{is} + A_- e^{-is},$
with $A_\pm = \frac{a\pm b}{2} \in \mathbb{R}$.
There are several virtues to this representation. If you want the ellipse to be centered around a point $z_0 = x_0 + i y_0$, just add it:
$z(s) \rightarrow z_0 + A_+ e^{is} + A_- e^{-is}$;
If you want the ellipse to be rotated by an angle $\psi$, just multiply:
$z(s) \rightarrow e^{i \psi}(A_+ e^{is} + A_- e^{-is}) = A_+ e^{i(s+\psi)} + A_- e^{-i(s-\psi)}.$
Application
If you are an engineer like I am, you are probably thinking of these equations in terms of phasors, which are complex numbers with fixed magnitude and linear phase (their phase changes at a constant rate). We are in the business of turning real signals into complex phasors for DSP applications. To do this, one must build a filter that has the effect of adding a $\pi/2$ phase shift. If you are not careful about how you build this filter, you can end up with an additional (backwards propagating) phasor that is identical to the one above. The result is a signal that traces out an ellipse, not a circle, in the complex plane. If you goof up the phase shift and get it wrong by a small amount ($\pi/2-\epsilon$), this equivalent to the above parametrization with $$\frac{A_-}{A_+} = \tan (\epsilon/2).$$
(The ellipse will also be rotated by an angle $\psi = \pi/4$.)
Best Answer
Write $\;z=a+ib\;,\;\;a,b\in\Bbb R\;$ , so:
$$z-5=(a-5)+ib\implies\left|\text{Re}\,(z-5)\right|=|a-5|\le1\iff 4\le a\le 6$$
and you have a pretty precise closed interval on the real axis.
Now,
$$\frac z{\overline z}=\frac{(a+ib)^2}{a^2+b^2}=\frac{a^2-b^2}{a^2+b^2}+i\frac{2ab}{a^2+b^2}\implies \text{Re}\,\frac z{\overline z}=\frac{a^2-b^2}{a^2+b^2}=1\implies\ldots\text{etc.}$$
Solve now.