An iterative Lagrange method for the regularization of discrete ill-posed inverse problems