>>407
Gは逆格子ベクトルそのものではなくて、その長さだから（ルート忘れてるみたいだが）
たんに逆格子点に相当するn1,n2,n3が分かっていればいい。

後は式の通りに書いて、逆格子点について足せばいい。

FORTRANの疑問ではないので、ここで聞く前に周囲の人に聞いた方が早いと思う。


数値微分は単に、格子間距離aを、特定のaとａ＋ｄａに取ったときに計算して
差分をｄａでわればいい。ｄａは小さい数だ。