A comparison between closed form stability functions and finite element stability functions shows that the former has numerical problems under small axial forces and the later is inaccurate under large a~ial forces. The least square fitting technique is applied to improve the finite element solution to get better accuracy under large axial forces and to solve numerical prebleIT's under small axial forces.