lin_g2.m

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

g = (x1 + x3 - 0.5) / 2;