Untitled

function[affidabilita, indici_aff] = calcolaAffidabili(energia, sottoprodotto)

    indici_aff = energia < 50 &  ~(sottoprodotto == 43 | sottoprodotto == 61 | sottoprodotto > 82);

    affidabilita = (sum(indici_aff) / length(energia)) * 100;
end

function [distr, indici_distr] = calcolaDistruttivi(energia, sottoprodotto)

    [~, indici_aff] = calcolaAffidabili(energia, sottoprodotto);

    indici_distr = ~indici_aff;

    distr = (sum(indici_distr) / length(energia)) * 100;
end

energia = input('Inserisci il vettore energia: ');
sottoprodotto = input('Inserisci il vettore sottoprodotti: ');

while length(energia) ~= length(sottoprodotto) || any(energia < 0) || any(sottoprodotto<0) ||
any(sottoprodotto ~= round(sottoprodotto))

    energia = input('Inserisci il vettore energia: ');
    sottoprodotto = input('Inserisci il vettore sottoprodotti: ');
else

    [affidabilita, indici_aff] = calcolaAffidabili(energia, sottoprodotto);

fprintf('Percentuale esperimenti affidabili: %f', affidabilita);


figure;
plot(find(indici_aff), energia(indici_aff), 'go');
hold on;
plot(find(~indici_aff), energia(~indici_aff), 'rx');

end