figure (1)

subplot(2,2,1)
hold on
plot(tout,y(:,2),'g')
ylabel('DP [bar]')
xlabel('Tid [s]')
grid on
hold off

subplot(2,2,2)
hold on
plot(tout,d(:,4),'g')
ylabel('mT_i_n_n [kg/s]')
xlabel('Tid [s]')
grid on
hold off

subplot(2,2,3)
hold on
plot(tout,y(:,1),'g')
ylabel('P_1 [bar]')
xlabel('Tid [s]')
grid on
hold off

subplot(2,2,4)
hold on
plot(tout,z,'g')
ylabel('z [-]')
xlabel('Tid [s]')
grid on
hold off