function [] = plotting(x,y,t)




figure(1)
clf
subplot(3,1,1)
plot(t/60,x(:,1));
title('mass of gas in annulus (x_1)')
subplot(3,1,2)
plot(t/60,x(:,2));
title('mass od gas in riser (x_2)')
subplot(3,1,3)
plot(t/60,x(:,3));
title('mass of liquid in riser (x_3)')



figure(2)
clf
subplot(4,2,1)
plot(t/60,y(:,1),'-r','linewidth',1);
%grid on
title('Injected gas into annulus [kg/s]')

subplot(4,2,2)
plot(t/60,y(:,2),'-r','linewidth',1);
%grid on
title('Pressure at top of annulus [bar]')

subplot(4,2,3)
hold on;
plot(t/60,y(:,3),'-g','linewidth',1);
%grid on
title('Pressure at top of riser [bar]')

subplot(4,2,4)
plot(t/60,y(:,4),'-b','linewidth',1);
%grid on
title('Pressure at bottom of annulus [bar]')

subplot(4,2,5)
plot(t/60,y(:,5),'linewidth',1);
%grid on
title('Pressure at bottom hole [bar]')

subplot(4,2,6)
plot(t/60,y(:,6),'linewidth',1);
%grid on
title('Choke mass flow [kg/s]')

subplot(4,2,7)
plot(t/60,y(:,7),'linewidth',1);
%grid on
title('Mixture density at top of riser [kg/m^3]')

subplot(4,2,8)
plot(t/60,y(:,8),'linewidth',1);
%grid on
title('Liquid volume fraction at top of riser [-]')

end