function relaxPt()

global n
global h
global v
global w
global omega
global sigma

for i=2:n
    w(i) = (1-omega)*v(i) + omega*(v(i-1)/(h^2)+v(i+1)/(h^2))/((2/h^2) + sigma);
end
    v = w;
%    s1(i) = norm(v2,inf);
