% Plotter resultatene som kommer fra simuleringen

kons=y(:,1:3);
temp=y(:,4);
n=length(t);
% Plotter utgangene
figure(1)
subplot(3,1,1)
plot(t,kons,t,ynom(1),'--')
grid
legend('xA','xC','xB')
ylabel('Molfraksjon')
xlabel('Tid')
title('Utganger')
axis([0 t(n) 0.9 1]);
subplot(3,1,2)
plot(t,temp,t,ynom(4),'--')
grid
ylabel('Temperatur [K]')
xlabel('Tid')
axis([0 t(n) -22 -10]);
subplot(3,1,3)
plot(t,y(:,5),t,y(:,6))
grid
ylabel('Holpup')
xlabel('Tid')
legend('Topp','Reboiler')
% axis([0 t(n) -22 -10]);


% Plotter pådragene
figure(2)
subplot(3,2,1)
plot(t,L)
legend('L')
grid
xlabel('Tid')
ylabel('Molstrøm')
axis([0 t(n) 0 1.2]);
subplot(3,2,2)
plot(t,V)
legend('V')
grid
xlabel('Tid')
ylabel('Molstrøm')
axis([0 t(n) 0 1.2]);
subplot(3,2,3)
plot(t,S)
legend('S')
grid
xlabel('Tid')
ylabel('Molstrøm')
axis([0 t(n) 0 1.2]);
subplot(3,2,4)
plot(t,Rl)
legend('Rl')
grid
xlabel('Tid')
ylabel('Molstrøm')
axis([0 t(n) 0 1.2]);
subplot(3,2,5)
plot(t,D)
legend('D')
grid
xlabel('Tid')
ylabel('Molstrøm')
axis([0 t(n) 0 1.2]);
subplot(3,2,6)
plot(t,B)
legend('B')
grid
xlabel('Tid')
ylabel('Molstrøm')
axis([0 t(n) 0 1.2]);