If you want $a_0$ to be $0$, then the curve you are trying to fit has equation $y = a_1 x + a_2 x^2$. If you derive the normal equations starting from this equation, you get

$a_1 \sum x_i^2 + a_2 \sum x_i^3 = \sum x_i y_i\\
a_1 \sum x_i^3 + a_2 \sum x_i^3 = \sum x_i^2 y_i$

Or

$\left[\begin{array}{cc|c}
\sum x_i^2 & \sum x_i^3 & \sum x_i y_i\\
\sum x_i^3 & \sum x_i^4 & \sum x_i^2 y_i
\end{array}\right]$

**Derivation of Normal Equations**  
$y = a_1 x + a_2 x^2$

The error $\epsilon_i$ in the $i$<sup>th</sup> term is the difference between the predicted value $y(x_i)$ and the observed value $y_i$:  
$\epsilon_i = y(x_i) - y_i = a_1 x_i + a_2 x_i^2 - y_i$

The total error is the sum of squares of these errors:  
$E = \displaystyle\sum_i \epsilon_i^2 = \displaystyle\sum_i(a_1x_i + a_2x_i^2 - y_i)^2$

To find the values of the parameters $a_1$ and $a_2$ that minimize the total error, let $\dfrac{\partial E}{\partial a_1} = \dfrac{\partial E}{\partial a_2} = 0$:

$
\dfrac{\partial E}{\partial a_1} = \sum 2x_i(a_1 x_i + a_2 x_i^2 - y_i) = 0 \Rightarrow a_1 \sum x_i^2 + a_2 \sum x_i^3 = \sum x_i y_i\\
\dfrac{\partial E}{\partial a_2} = \sum 2x_i^2(a_1 x_i + a_2 x_i^2 - y_i) = 0 \Rightarrow a_1 \sum x_i^3 + a_2 \sum x_i^4 = \sum x_i^2 y_i\\
$