MATLABを使ってゼータ関数の零点を計算する。\(s\)を\(s=\frac{1}{2}+yi\)とするとゼータ関数のノルムは
\begin{align}
\zeta(s)= \left \|\sum_{i=1}^{\infty} \frac{1}{n^s} \right \|
\end{align}
となる。MATLABの関数を使って描画すれば
以下コード。
Nmin=10;
Nmax=20;
d=0.01;
N=14;
n=50000;
j=sqrt(-1);
result=zeros(1,N);
y1=zeros(1,N);
count=1;
for y=Nmin:d:Nmax
s=1/2+y*j;
result(1,count)=norm(zeta(s));
count=count+1;
end
fig1 = figure('name', 'ゼータ関数の特殊値');
plot(Nmin:d:Nmax,result,'k-')
grid on
xlabel('$y$','Interpreter', 'latex');
ylabel('$\zeta(s)$','Interpreter', 'latex');
コメント