figure(1)

subplot(3,1,1)
hold on
grid on
plot(tout,y(:,1),'b')
ylabel('P_1 [bar]')
hold off

subplot(3,1,2)
hold on
grid on
axis([0 1000*20 0 4])
plot(tout,y(:,2),'b')
ylabel('DP [bar]')
hold off

subplot(3,1,3)
hold on
grid on
plot(tout,z,'b')
ylabel('z[-]')
hold off
