lin_g3.m

function g = lin_g3(x1,x2,x3,x4);
% function used with script jacobi.m

g = (x2 + x4 - 0.5) / 2;