alphaM = 0.209;
nSim = 9;
vect = [1:nSim]';
Ju = [1:nSim];
Jx = [1:nSim];

for i = 1:nSim,
    vect(i) = alphaM/(10-i);
end


for y = 1:9,
    nPassi = 1000; % durata della simulazione
    tempi = [1:nPassi]';

    % Perturbazione sul livello delle scorte
    x = zeros(nPassi,1);

    % Perturbazione sulla domanda che il sito fa al sito a monte
    u = zeros(nPassi,1);

    % Perturbazione sulla domanda che il sito riceve dal sito a valle
    omega = 0.4;
    d = sin(omega*tempi);

    % Ritardo (trasporto+produzione)
    lambda = 7;

    % Costante alfa della politica P: u_k = -alfa*x_k
    alfa = vect(y);

    for k = 1:nPassi-1,

        u(k) = -alfa*x(k);

        if k-lambda > 0,
            x(k+1) = x(k) + u(k-lambda) - d(k);
        else x(k+1) = x(k) - d(k);
        end

    end

    % Prestazioni:
    Jx(y) = sqrt(sum(x.^2))/sqrt(sum(d.^2)) %Jx = rapporto tra norma 2 di x e norma 2 della domanda d dal sito a valle
    Ju(y) = sqrt(sum(u.^2))/sqrt(sum(d.^2)) %Ju = rapporto tra norma 2 della domanda u al sito a monte e norma 2 della domanda d dal sito a valle
end

figure
plot(vect,Jx,'b',vect,Ju,'r')