function?[dfdeta]?=?free_energ_fft_ca_v2（NxNyngrainetasetaigrain）

format?long;

A=1.0;
B=1.0;

sum=zeros（NxNy）;

for?jgrain=1:ng