Update GrainGrowth.m

Correction in the derivative of the free energy functional.

Author: abhinavroy1999

GrainGrowth.m

@@ -131,25 +131,25 @@ function GrainGrowth(Nx, Ny, dx, dy, end_time, time_step)
 for temp = start_time : end_time
   for i = 1 : Nx
       for j = 1 : Ny
-            geta1(i,j) = -eta1(i,j) + eta1(i,j)*eta1(i,j)*eta1(i,j)*eta1(i,j) ... 
+            geta1(i,j) = -eta1(i,j) + eta1(i,j)*eta1(i,j)*eta1(i,j) ... 
                         + 2*eta1(i,j)*(eta2(i,j) + eta3(i,j) + eta4(i,j) + eta5(i,j) + eta6(i,j) + eta7(i,j) + eta8(i,j) + eta9(i,j) + eta10(i,j));
-            geta2(i,j) = -eta2(i,j) + eta2(i,j)*eta2(i,j)*eta2(i,j)*eta2(i,j) ... 
+            geta2(i,j) = -eta2(i,j) + eta2(i,j)*eta2(i,j)*eta2(i,j) ... 
                         + 2*eta2(i,j)*(eta1(i,j) + eta3(i,j) + eta4(i,j) + eta5(i,j) + eta6(i,j) + eta7(i,j) + eta8(i,j) + eta9(i,j) + eta10(i,j));
-            geta3(i,j) = -eta3(i,j) + eta3(i,j)*eta3(i,j)*eta3(i,j)*eta3(i,j) ... 
+            geta3(i,j) = -eta3(i,j) + eta3(i,j)*eta3(i,j)*eta3(i,j) ... 
                         + 2*eta3(i,j)*(eta2(i,j) + eta1(i,j) + eta4(i,j) + eta5(i,j) + eta6(i,j) + eta7(i,j) + eta8(i,j) + eta9(i,j) + eta10(i,j));
-            geta4(i,j) = -eta4(i,j) + eta4(i,j)*eta4(i,j)*eta4(i,j)*eta4(i,j) ... 
+            geta4(i,j) = -eta4(i,j) + eta4(i,j)*eta4(i,j)*eta4(i,j) ... 
                         + 2*eta4(i,j)*(eta2(i,j) + eta3(i,j) + eta1(i,j) + eta5(i,j) + eta6(i,j) + eta7(i,j) + eta8(i,j) + eta9(i,j) + eta10(i,j));
-            geta5(i,j) = -eta5(i,j) + eta5(i,j)*eta5(i,j)*eta5(i,j)*eta5(i,j) ... 
+            geta5(i,j) = -eta5(i,j) + eta5(i,j)*eta5(i,j)*eta5(i,j) ... 
                         + 2*eta5(i,j)*(eta2(i,j) + eta3(i,j) + eta4(i,j) + eta1(i,j) + eta6(i,j) + eta7(i,j) + eta8(i,j) + eta9(i,j) + eta10(i,j));
-            geta6(i,j) = -eta6(i,j) + eta6(i,j)*eta6(i,j)*eta6(i,j)*eta6(i,j) ... 
+            geta6(i,j) = -eta6(i,j) + eta6(i,j)*eta6(i,j)*eta6(i,j) ... 
                         + 2*eta6(i,j)*(eta2(i,j) + eta3(i,j) + eta4(i,j) + eta5(i,j) + eta1(i,j) + eta7(i,j) + eta8(i,j) + eta9(i,j) + eta10(i,j));
-            geta7(i,j) = -eta7(i,j) + eta7(i,j)*eta7(i,j)*eta7(i,j)*eta7(i,j) ... 
+            geta7(i,j) = -eta7(i,j) + eta7(i,j)*eta7(i,j)*eta7(i,j) ... 
                         + 2*eta7(i,j)*(eta2(i,j) + eta3(i,j) + eta4(i,j) + eta5(i,j) + eta6(i,j) + eta1(i,j) + eta8(i,j) + eta9(i,j) + eta10(i,j));
-            geta8(i,j) = -eta8(i,j) + eta8(i,j)*eta8(i,j)*eta8(i,j)*eta8(i,j) ... 
+            geta8(i,j) = -eta8(i,j) + eta8(i,j)*eta8(i,j)*eta8(i,j) ... 
                         + 2*eta8(i,j)*(eta2(i,j) + eta3(i,j) + eta4(i,j) + eta5(i,j) + eta6(i,j) + eta7(i,j) + eta1(i,j) + eta9(i,j) + eta10(i,j));
-            geta9(i,j) = -eta9(i,j) + eta9(i,j)*eta9(i,j)*eta9(i,j)*eta9(i,j) ... 
+            geta9(i,j) = -eta9(i,j) + eta9(i,j)*eta9(i,j)*eta9(i,j) ... 
                         + 2*eta9(i,j)*(eta2(i,j) + eta3(i,j) + eta4(i,j) + eta5(i,j) + eta6(i,j) + eta7(i,j) + eta8(i,j) + eta1(i,j) + eta10(i,j));
-            geta10(i,j) = -eta10(i,j) + eta10(i,j)*eta10(i,j)*eta10(i,j)*eta10(i,j) ... 
+            geta10(i,j) = -eta10(i,j) + eta10(i,j)*eta10(i,j)*eta10(i,j) ... 
                         + 2*eta10(i,j)*(eta2(i,j) + eta3(i,j) + eta4(i,j) + eta5(i,j) + eta6(i,j) + eta7(i,j) + eta8(i,j) + eta9(i,j) + eta1(i,j));        
       end
   end