How to realize Riemann Siegel Theta Function with MATLAB?
5 ビュー (過去 30 日間)
古いコメントを表示
How to realize Riemann Siegel Theta Function with MATLAB?
0 件のコメント
採用された回答
John D'Errico
2020 年 6 月 1 日
編集済み: John D'Errico
2020 年 6 月 2 日
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 件のコメント
兰花 西
2020 年 8 月 24 日
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 !
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Startup and Shutdown についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!