x=[1 2 3 4 5];
y=[0.2713e5 0.3066e5 0.3836e5 5.4803e5 -1.5354e5];

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