diff --git a/gpu/Color.cu b/gpu/Color.cu
index 4c0280da..795ca258 100644
--- a/gpu/Color.cu
+++ b/gpu/Color.cu
@@ -1241,12 +1241,12 @@ __global__  void dvc_MassColorCollideD3Q7(char *ID, double *A_even, double *A_od
 			A_even[n] = 0.3333333333333333*na;
 			B_even[n] = 0.3333333333333333*nb;
 			// Non-Stationary equilibrium distributions
-			feq[0] = 0.1111111111111111*(1+3*ux);
-			feq[1] = 0.1111111111111111*(1-3*ux);
-			feq[2] = 0.1111111111111111*(1+3*uy);
-			feq[3] = 0.1111111111111111*(1-3*uy);
-			feq[4] = 0.1111111111111111*(1+3*uz);
-			feq[5] = 0.1111111111111111*(1-3*uz);
+			feq[0] = 0.1111111111111111*(1+4.5*ux);
+			feq[1] = 0.1111111111111111*(1-4.5*ux);
+			feq[2] = 0.1111111111111111*(1+4.5*uy);
+			feq[3] = 0.1111111111111111*(1-4.5*uy);
+			feq[4] = 0.1111111111111111*(1+4.5*uz);
+			feq[5] = 0.1111111111111111*(1-4.5*uz);
 			// Construction and streaming for the components
 			for (idx=0; idx<3; idx++){
 				//...............................................