aux=open('solOPT.mat');
sol=aux.solOPT;



% solMH.t=Timevec;
% solMH.qhf=qhf;
% solMH.qhr=qhr;
% solMH.mr=mr;
% solMH.qair=qair;
% solMH.Tf=Tf;
% solMH.Tr=Tr;
% solMH.co2=co2;
% solMH.J=EnergyCost;
% solMH.price=price;
% solMH.To=To;
% solOPT=solMH;





qinsim.time= sol.t;
qinsim.signals.values=sol.qair;
qinsim.signals.dimensions=[1];


Tosim.time= sol.t;
Tosim.signals.values=sol.To;
Tosim.signals.dimensions=[1];


Pricesim.time= sol.t;
Pricesim.signals.values=sol.price;
Pricesim.signals.dimensions=[1];



[init,par]=getData
sim('testFeedback1.mdl')
Linewidth=1.5;
% figure(4)
% plot(timefeedback/3600,costfeedback,'g','LineWidth',Linewidth)
