% solve for V for a capacitor

% mesh size
M = 70;
N = 7;

% initialize
V0 = zeros(M,N);
%V0 = ones(M,N)*7;

% bondary conditions
V0(:,1) = -1;
V0(:,end) = 1;
V1 = V0;

% exit tolerance & take one step 
tol = 1e-5;
n = 1
V1 = laplaceStepJacobi(V1,V0)


while ...
        sum(sum(abs(V1(2:end-1,2:end-1)-V0(2:end-1,2:end-1)))) ...
        / ... 
        ((M-2)*(N-2)) > tol
%for k = 1:1005
    V0 = V1;
    n = n+1
    V1 = laplaceStepJacobi(V1,V0)
end





return