Abstract
Algorithms for calculating sequences of upper and lower bounds for an estimate of the optimal backward error for linear least squares problems are developed, based on the Gauss quadrature theory. Numerical results show that the bounds converge quickly and are therefore useful in practice.