x=[1 2 3 4 5];
y=[7.2142e5 -1.3958e5 -1.3958e5 -0.9390e5 -0.9390e5];

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