gridres = 10;
xv = ones(gridres);
yv = zeros(gridres);
vv = zeros(gridres);
nu = 0.8;
rho = 1;

for iterations = 1:10
  assert;
  nxv = xv;
  nyv = yv;
  nvv = vv;
  for i = 2:gridres-1;
    for j = 2:gridres-1;
      nxv(i,j) = ((xv(i-1,j) + xv(i+1,j)) + nu*( xv(i,j-1) + xv(i,j+1)))/4;
      nyv(i,j) = ((yv(i-1,j) + yv(i+1,j)) + nu*( yv(i,j-1) + yv(i,j+1)))/4;
    end
  end
  xv = nxv;
  yv = nyv;
end

imshow(xv,xv,xv)


