function drawKYPhase(filename,clrDat,clrInt)

data = load(filename);
Kdat = data(:,1); ydat = data(:,2);
K = 2.35:0.005:3;
y = interp1(Kdat,ydat,K,'cubic');


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