说明：模拟二维情况下的晶粒演化，将计算区域离散为正方形网格，采用周期性边界条件。初始相场值为在（-0.001-0.001）之间服从正态分布的随机数，模拟中引入的第二相颗粒近似为圆形。模拟中参数如下：p=36，dx=l，dt=0．1，k=0．5，L=1.0，m=1．0，网格节点数为300x300。

matlab 求解