Algorithmic construction of R-optimal designs for second-order response surface models