相场法建模

function [dfdeta] = free_energ_fft_ca_v2（NxNyngrainetasetaigrain）

format long;

A=1.0;
B=1.0;

sum=zeros（NxNy）;

for jgrain=1:ng