function relaxVPt()

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

v(2:n) = (1-omega)*v(2:n) + omega*(v(1:n-1)/(h^2)+v(3:n+1)/(h^2))/((2/h^2) + sigma);
