あとはE2(Ω,x)=1/2(4πe**2/oG**2)**-1(V(G))**2*1/(1-f(G))*(1/ε(G)-1)
狽ﾍG/=0 を数値微分して全部足すと出るんですけどここがわからないです。
G=2π/a(n1*n1+n2*n2+n3*n3) 収束すると思うのである程度までの結果を出していきたいです。
n1=n2=n3/=0の時以外を計算させる方法がわかりません。数値微分もわかりません。
V(G)=-4*Z*2/oG**2((1+uk)cos(G*Rmk)-uk*sin(G*Rm)/G*Rm)
ε(G)=1+(1-f(G))*(4*3.14*z*2/o*G**2)*(3/2*Ef)*(1/2+(4*Kf**2-G**2)/(8*Kf*G)*log((2*Kf+G)/2*Kf-G)
f(G)=G**2/4*(G**2+c*Kf**2)+G**2/4*c*Kf**2
です。