clc
clear all
%%
%%%%%Valve1%%%%%

load valve_1
z1 = valve_1 (:,1);
P1_ave = valve_1 (:,2);
P1_max = valve_1 (:,3);
P1_min = valve_1 (:,4);

figure (1)
plot (z1,P1_ave,'k','LineWidth',3);
hold on
plot (z1,P1_max,'k','LineWidth',3);
hold on
plot (z1,P1_min,'k','LineWidth',3);
grid on
xlabel ('Z');
ylabel('Inlet Pressure [Kpa]');
title ('Open loop Bifurcation Diagram');

%%
%%%%%valve_2%%%%%

load valve_2
z2 = valve_2 (:,1);
P2_max = valve_2 (:,2);
P2_ave = valve_2 (:,3);
P2_min = valve_2 (:,4);

figure (2)
clf
plot (z2,P2_max,'k','LineWidth',3);
hold on
plot (z2,P2_ave,'k','LineWidth',3);
hold on
plot (z2,P2_min,'k','LineWidth',3);
grid on
xlabel ('Z');
ylabel('Inlet Pressure [Kpa]');
title ('Open loop Bifurcation Diagram');
