function dx = cell_equation(t,X,P)


%reshape the 1xN input for the size of the image
x=reshape(X,P.SizeX,P.SizeY);

dx=zeros(P.SizeX,P.SizeY);


%reshape back from Nx1
dx=reshape(dx,P.SizeX*P.SizeY,1);