x=[1 2 3 4 5];
y=[7.2142e5 -1.3958e5 -1.3958e5 -0.9390e5 -0.9390e5];
z=[7.1212e5 -1.2227e5 -1.2227e5 -0.8213e5 -0.8213e5];
t=[5.4513e5 2.6741e5 -0.9480e5 -0.6233e5 -0.6233e5];
u=[7.3981e5 -1.72968e5 -1.7268e5 -1.1650e5 -1.1650e5];
plot(x,y,'-*',x,z,'-o',x,t,'--',x,u,'--')
xlabel('mode number'),ylabel('\omega^2_p')
title('P=\rho+\kappalog\rho')