Simple and accurate approximations to the Riemann zeta function

We develop approximations for the Riemann zeta function that enable high-precision computation within the critical strip and other vertical strips. These approximations combine the main sum of the Riemann-Siegel formula with a simple approximation of the remainder term, which involves only elementary functions and certain precomputed coefficients obtained via Gaussian quadrature. Additionally, we provide approximations for the derivative of the Riemann zeta function and present extensive numerical evidence demonstrating the accuracy of these approximations.