function V = uq_OvertoppingFailure2(X)

gamma = 0.5;
SWL = 5;
X1 = 3.3606e+00; 
X2 = 1.1791e+00;
X3 = 1.0221e+01;
Rc = X1 + X2 + X3 - SWL - X(:, 3);
Tm = X(:, 2)./1.25;
Lm = 9.81.*(Tm.^2)/2/pi();
Hs = X(:, 1);
cot = 2.00515;
V = 0.001 - 0.026.*sqrt(Lm./Hs./cot).*...
    exp(-((2.5.*Rc.*cot./sqrt(Hs.*Lm)./gamma).^1.3))...
    .* sqrt(9.81.*(Hs.^3));