x=[1 2 3 4 5];
y=[0.1291e5 0.1455e5 0.1812e5 6.1859e5 2.4156e5];


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