clc
clear all
close all

num=5; %case

%% Load data
namePI=strcat('PI',num2str(num));
nameMPC=strcat('MPC',num2str(num),'_newsim');
load(namePI,'-mat')

load(namePI,'-mat')
load(nameMPC,'-mat')


%Calculate integral(abs(delta qout/delta y))
integralMPC=0;
integralPI=0;
deltat=0.1;
steps=300;
for j=1:steps
%     delta=sim.u(j+1,1)-sim.u(j,1);
%     deltaydeltat=abs(delta/deltat);
%     integralMPC=integralMPC+deltaydeltat;  
    delta2=q_out.signals.values(j+1)-q_out.signals.values(j);
    deltaydeltat2=abs(delta2);
    integralPI=integralPI+deltaydeltat2;
end 
for k=1:steps
    deltaMPC(k)=abs(sim.u(k+1)-sim.u(k));
    % deltaydeltat=abs(delta);
   % integralMPC=integralMPC+deltaydeltat;
end 
integralMPC=sum(deltaMPC)/deltat

%integralMPC
integralPI=integralPI/deltat