function V = uq_OvertoppingFailure(X)
% UQ_SIMPLYSUPPORTEDBEAM computes the overtopping failure
%
%   Model with five input parameters  X= [cot Hs Tp WL]
%         cot:  relative density
%         Hs:   significant wave height
%         Tp:
%         WL:
%
%   Output:  V = 0.001 - 0.026.*sqrt(Lm./Hs./cot).*...
%    exp(-(((2.5.*(hc-WL).*cot.*sqrt(Hs./Lm))./Hs./gamma).^1.3)).*...
%    sqrt(9.806.*(Hs.^3));
%    Tm = Tp/1.1;
%    Lm = 9.81*(Tm.^2)/2/pi();
% 
% Vectorized implementation
gamma = 0.5;
hc = 11.55;
Tm = X(:, 6)./1.25;
Lm = 9.81.*(Tm.^2)/2/pi();
V = 0.001 - 0.026.*sqrt(Lm./X(:, 5)./X(:, 4)).*...
    exp(-(((2.5.*(hc-X(:, 9)).*X(:, 4).*sqrt(X(:, 5)./Lm))./X(:, 5)./gamma).^1.3))...
    .* sqrt(9.806.*(X(:, 5).^3));