function tvTuning = tvTuning(gp,minTuning)



kc=minTuning(1,:); % loading PO curve PI settings
ti=minTuning(2,:);

tvTuning = [];
for i = 1:length(kc)

gc = controller(kc(i),ti(i));

Gu  = feedback(gc,gp);
Gdu = feedback(-gp*gc,1,1);
sys = [Gu;Gdu];
[u,t]=step(sys);
u = [0 0;u];
tvTuning(1,i) = tv(u(:,1));
tvTuning(2,i) = tv(u(:,2));

end

return