linear regression versus a system of equations

It turns out that the closed form linear regression solution of Y predicted with X with coeffs beta = (X^tX)^-1X^tY is the same as the one that would be obtained by differentiating the sum of squared deviations with each parameter in beta and solving the system of linear equations... courtesy Rakesh Kumar