x=[1 2 3 4 5];
y=[0.0704e5 0.0792e5 0.0984e5 0.0765e5 0.0302e5];
plot(x,y,'-*')
title('p=\rho+log\rho')
xlabel('mode number'),ylabel('\omega^2')