function flopKYDraw(filename,clr1,clr2)

data = load(filename);
Kp = data(:,1); yp = data(:,2);
K = 1.76:0.005:2.15;
y = interp1(Kp,yp,K,'cubic');


hold on;
if strcmp(clr1,'np') == 0
    plot(Kp,yp,strcat('o',clr1),'MarkerSize',10);
end
if strcmp(clr2,'np') == 0
    plot(K,y,clr2,'LineWidth',2);
end
grid on;
title('y-K_{c} Phase Diagram');
xlabel('K_{c}');
ylabel('y');
