function V1 = laplaceStepJacobi(V1,V0)
%V1 = laplaceStepJacobi(V0);
% assumes that B.C. specify V on the exterior edges in 2D

for ii = 2:size(V0,1)-1
    for jj = 2:size(V0,2)-1
        V1(ii,jj) = ...
            (V0(ii+1,jj) + V0(ii-1,jj) + V0(ii,jj+1) + V0(ii,jj-1))/4;
    end
end

return