


load press_dist
figure
hold on
plot(press_dist(:,1),press_dist(:,2),'r','linewidth',2)
hold all
plot(press_dist(:,1),press_dist(:,3),'b','linewidth',2)
legend('Pressure(nominal)','Pressure(10% disturbance)')
xlabel 'Valve Opening'
ylabel 'Pressure [bar] '
title 'Pressure sensitivity at 10% Disturbance in Reservoir Gas Fraction'
%axis([20 110 0 3000])
hold off

