sim_t = 10;

k = 2;
T = 3;

Kp = 2;
Ti = 0.5;

figure;
color_ptr = 0;

x = 0.1:0.1:1;
plots = numel(x);
for i = 0.1:0.1:1
    Ti = T;
    Kp = T / (i * T * k);
    
    sim ('pid_model1.slx');
    plot (tout, y, 'Color', [1 - color_ptr * 1 / plots, 0, 0]);
    hold on;
    color_ptr = color_ptr + 1;
end

sim ('model1.slx');

plot (tout, y);