function z = prob_grooteman_09_parallel_x2z_test(varargin)

%% settings

OPT = struct(...
    'U1', [],...
    'U2', [],...
    'U3', [],...
    'U4', [],...
    'U5', []...
    );

OPT = setproperty(OPT, varargin{:});

%% limit state function

z1  = 2.677-OPT.U1-OPT.U2;
z2  = 2.500-OPT.U2-OPT.U3;
z3  = 2.323-OPT.U3-OPT.U4;
z4  = 2.250-OPT.U4-OPT.U5;
z   = max([z1, z2, z3, z4], [], 2);