x=[1 2 3 4 5];
y=[0.1267e5 0.2219e5 0.2771e5 5.6971e5 2.2572e5];

z=[0.2713e5 0.3066e5 0.3836e5 5.4803e5 -1.5354e5];
plot(x,y,'-*',x,z,'--')
xlabel('mode number'),ylabel('\omega^2_p')
title('P=\rho+\kappalog\rho')