23/02/06 07:27:05.79 BB/Q99JB.net
>>517
kx=1e4
ky=3e3
xx=seq(0,2*pi,length=kx)
yy=seq(0,2,length=ky)
re=0
for(i in xx){
for(j in yy){
re=re+f(i,j,FALSE)
}
}
re/(kx*ky)*2*pi*2
> re/(kx*ky)*2*pi*2
[1] 9.422767
理論値3πに近似した値が得られた。
23/02/06 07:27:05.79 BB/Q99JB.net
>>517
kx=1e4
ky=3e3
xx=seq(0,2*pi,length=kx)
yy=seq(0,2,length=ky)
re=0
for(i in xx){
for(j in yy){
re=re+f(i,j,FALSE)
}
}
re/(kx*ky)*2*pi*2
> re/(kx*ky)*2*pi*2
[1] 9.422767
理論値3πに近似した値が得られた。