function Y = uq_Ex08_BermBreakwater_robustness(X, myModel, myInput)
    nSamples = 1000;
    Y = zeros(size(X,1), 1);
    for i = 1:size(X,1)
        B = X(i,1);
        hb = X(i,2);
        Ns = X(i,3);
        hf = X(i,4);
        envSamples = uq_getSample(myInput, nSamples);
        X_full = [repmat([B, hb, Ns, hf], nSamples, 1), envSamples];
        limitStateOutput = uq_evalModel(myModel, X_full);
        meanOutput = mean(limitStateOutput(:,1));
        stdOutput = std(limitStateOutput(:,1));
        CoV = stdOutput / abs(meanOutput + eps); % Add eps to avoid division by zero
        Y(i) = 1 / (CoV + eps); % Robustness is inversely proportional to CoV
    end
end