% Temperature profile on stages in the prefractionator and the main column
stage=1:48;
stage1=stage(13:36);
Temp_pre=temp_nom(1:24);
plot(stage1,Temp_pre,'b')
hold on
Temp_main=temp_nom(25:72);
plot(stage,Temp_main,'r')

%Plot Yopt (change in temp profile because of disturbances)
%stage=1:48;
%stage1=stage(13:36); 
%delYopt1 = z2(1:24);
%plot(stage1,delYopt1,'b')
%hold on
%delYopt2=z2(25:72);
%plot(stage,delYopt2,'r')




%hold on
%axis([0 2*pi -1 1])
%X1=1:100;
%Y1=sqrt(X1);
%X2=1:100;
%Y2=cos(X2);
%figure;
%h1=subplot(2,1,1);
%h2=subplot(2,1,2);

%plot(h1,X1,Y1);
%plot(h2,X2,Y2);
%http://www.matlabtips.com/advanced-plotting-in-matlab/
%linkaxes([h1 h2],'x');
%http://jitkomut.lecturer.eng.chula.ac.th/matlab/graphic.html