x=[1 2 3 4 5];
y=[7.1212e5 -1.2227e5 -1.2227e5 -0.8213e5 -0.8213e5];

plot(x,y,'-o')
xlabel('mode number'),ylabel('\omega^2_p')
title('P=\rho+10log\rho')