lin_g1.m

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

g = (x2 + 1.5) / 2;