function tvCombined = tvTuning(gp,minTuning)



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

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

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

Gu  = feedback(gc,gp);
Gdu = feedback(-gp*gc,1,1);
sys = [Gu;Gdu];
[y,t]=step(sys);

tvCombined(i) = tv(y(:,1)) + tv(y(:,2));

end

return