すいません。あほなんで基本的には全然わかりません。
状態方程式のＰを求めるためにやってると思うんですけど。
圧力Ｐは結晶エネルギーＥvca(Ω,x)の体積Ωに関係する微係数で与えられるらしいのです。
a=5.225A A=10**-8
Ω=(a**3)/2 rs=(9π/4)**(1/3)/Kf
Ef=h**2*Kf**2/2m c=2(1+0.0155π/Kf)**-1
f(q)=q**2/2(Q**2+cKf**2)
ε(q)=1+(1-f(q))*(4πZe**2/Ωq**2)*(3/2Ef)*(1/2+((4Kf**2-q**2)/8Kfq)*log(2Kf+q/2Kf-q))
CsのＨ近似の場合Rmは4.610 uは-0.9605
KのＨ近似の場合Rmは3.019 uは-0.5682

せめて濡/=0 のＧベクトルが０以外の時を全部足す方法が知りたいです。
お願いします。