function z = prob_grooteman_14_series_multiple_DP_x2z_test(varargin)
%% settings

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

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

%% limit state function

z1  = 0.1*(OPT.U1-OPT.U2).^2-(OPT.U1+OPT.U2)/sqrt(2)+3;
z2  = 0.1*(OPT.U1-OPT.U2).^2+(OPT.U1+OPT.U2)/sqrt(2)+3;
z3  = OPT.U1-OPT.U2+3*sqrt(2); 
z4  = OPT.U2-OPT.U1+3*sqrt(2);
z   = min([z1, z2, z3, z4], [], 2);