d1 = load('../data/flopPhDat01_0a.m'); % D = 0.01
d2 = load('../data/flopPhDat02_0a.m'); % D = 0.03
d3 = load('../data/flopPhDat03_0a.m'); % D = 0.05

K1 = d1(:,1); y1 = d1(:,2);
K2 = d2(:,1); y2 = d2(:,2);
K3 = d3(:,1); y3 = d3(:,2);

K = 1.76:0.005:2.15;

yIp1 = interp1(K1,y1,K,'cubic');
yIp2 = interp1(K2,y2,K,'cubic');
yIp3 = interp1(K3,y3,K,'cubic');

plot(K,yIp1,'k',K,yIp2,'g',K,yIp3,'b','LineWidth',2);
grid on;
title('y-K_{c} Phase Diagram');
xlabel('K_{c}');
ylabel('y');
legend('D_{b} = 0.01','D_{b} = 0.03','D_{b} = 0.05');