In this paper, we develop a direct formula for determining the coefficients in canonical basis of best polynomial degree $M$ that approximates $N>M$ on symmetric interval $\mathcal{L}^2$-norm. We also formally prove using is more computationally efficient than classical matrix multiplication approach and provide an example to illustrate it numerically stable approach.
