for i = 2:ny
  Vinl = [];
  for j = 2:nx
    vals = [V(i - 1,j-1:j), V(i, j-1:j)];
    datMat = [vals; xDeriv(V, i, j); yDeriv(V, i , j)];
    if(j == 2)
      Vinl = [Vinl, interpolationFunction(datMat, len)];
    else
      Vinl = [Vinl, interpolationFunction(datMat, len)(:, 2:end)];
    end
  end
  if(i == 2)
    Vinter = [Vinter; Vinl];
  else
    Vinter = [Vinter; Vinl(2:end, :)];
  end
end