lin_g4.m

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

g = (x3 + 0.5) / 2;