% Plotter resultatene som kommer fra simuleringen

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


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

% % Plotter konsentrasjonene ved målepunktene
% figure(3)
% subplot(2,1,1)
% plot(t,y(:,10))
% ylabel('Molfraksjon C øvre målepkt')
% subplot(2,1,2)
% plot(t,y(:,11))
% ylabel('Molfraksjon A nedre målepkt')
