In the discrete formulation of the bubble stabilized Legendre Galerkin methods, the system of equations includes the artificial viscosity term as the parameter. We investigate the estimation of this parameter to get the least-squares solution which minimizes the sum of the squares of errors at each node points. Some numerical results are reported.