implemented the Riemann-Siegel theta function (for real series)