How to realize Riemann Siegel Theta Function with MATLAB?
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
yifei wang
el 1 de Jun. de 2020
Comentada: 兰花 西
el 24 de Ag. de 2020
How to realize Riemann Siegel Theta Function with MATLAB?
0 comentarios
Respuesta aceptada
John D'Errico
el 1 de Jun. de 2020
Editada: John D'Errico
el 2 de Jun. de 2020
You start by writing the formula for it. That can be found on virtually the first line in this link.
However in MATLAB, the gamma function seems to be defined only for a real variable, unless you use the symbolic toolbox version of gamma. So this should work:
t = (1:5)';
[t,angle(double(gamma(sym(1/4 + i*t/2)))) - log(pi)/2*t]
ans =
1 -1.76754795281229
2 -2.52591091881613
3 -2.99456469601083
4 -3.29063503121648
5 -3.45962037536346
That is valid for real arguments t. I'd need to do some reading to know if this works for complex arguments, though I assume it does not, since the wiki link I gave indicates it is valid for real arguments. Hopefully you care only about real arguments, and are uninterested in an extension to the complex plane.
4 comentarios
兰花 西
el 24 de Ag. de 2020
Excuse me, when I set the parameter to a relatively large imaginary number(-1.93097838020833 + 38142819340462.3 i) in the above formula, I cannot get a valid result. Do you know how to solve it? Thank you !
Más respuestas (0)
Ver también
Categorías
Más información sobre Startup and Shutdown en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!