The answer is much simply presented in the book by John D Anderson - Modern Compressible Flow with historical perspectives. Refer Chapter 5.
\begin{equation}
\rho = \mbox{density}
\end{equation}
a = speed of sound, Ma = Mach number, u = speed of fluid, A = area
Area velocity relationship can be deduced from the continuity relationship -
\begin{equation}
\frac{d\rho}{\rho} + \frac{du}{u} + \frac{dA}{A} = 0
\end{equation}
Watch now ...
\begin{equation}
\frac{dA}{A} = - \{\frac{d\rho}{\rho} + \frac{du}{u}\}
\end{equation}
From momentum conservation relationships between two points, we can write
\begin{equation}
dp = -\rho udu \\
\implies \frac{dp}{\rho} = \frac{dp}{d\rho}\frac{d\rho}{\rho} = -udu
\end{equation}
and
\begin{equation}
\frac{dp}{d\rho} = a^2
\end{equation}
Therefore,
\begin{equation}
\frac{d\rho}{\rho} = -\frac{udu}{a^2} = -Ma^2 \frac{du}{u}
\end{equation}
From these relations, we can write
\begin{equation}
\frac{dA}{A} = (Ma^2 - 1)\frac{du}{u}
\end{equation}
Now check for Ma = 1, you will get
\begin{equation}
\frac{dA}{A} = 0
\end{equation}
This means the area function is either at maxima or minima. The practical solution is that this is minimum. This minimum area is called throat of the nozzle.