Hi,
For real argument, we could easily interface std::riemann_zeta :
https://en.cppreference.com/w/cpp/numeric/special_functions/riemann_zeta
If you have a compiler (under windows you can install the minGW atoms
module), you can run the following script:
code=[
"#define
Finally got to the end of the problem and replicated the plot of the
Riemann Zeta function on the critical line (s=0.5 + %i*t)
Looks pretty close to that shown on the Wikipedia page for the Riemann Zeta
Function!
function zs1=zeta_0_1(s, n)
// Vectorised versionzs1=0k=linspace(1,n,n);
Hi all,
After a lot of trial and error, I have managed to get a set of functions to
compute the approximations of Riemann's Zeta for negative and positive real
values; values of n > 1e6 seem to give better results:
function zs=zeta_s(z, n)
// Summation loop
zs=1;
if z == 0
zs