Evca(Ω,x)を求めるプログラムを作りたいのですが教えてください。
Evca(Ω,x)=Ei(Ω,x)+E0(Ω,x)+E1(Ω,x)+E2(Ω,x)
Ei(Ω,x)=-1.79186*Z**(5/3)/r
E0(Ω,x)=Z(2.21/r**2-0.916/r+Ecorr)
Ecorr=-(0.115-0.031logr)
E1(Ω,x)=2πZe**2((1-X)Z**Cs*(Rm**Cs)**2(1+2/3*u**Cs)+xZ**K(Rm**k)**2(1+2/3u**K))
E2(Ω,x)=1/2(4πe**2/ΩG**2)**-1(V(G))**2*1/(1-f(G))*(1/ε(G)-1)
狽ﾍG/=0
Gは逆格子ベクトルG(n1,n2,n3)
V(q)=(1-x)V**Cs+xV**k(q)
Vi(q)=-4πZe**2/Ωq**2((1+u)cos(qRm**i)-usin(qRm**i)/qRm) (i=Cs,K)
u=qRmcos(qRm)/sin(qRm)-qRmcos(qRm)
わからない値は後から自分で入れるので教えてください
状態方程式の計算をするらしいです