Based on hierarchic shape functions and an effective convergence procedure, the p-version and h-p adaptive analysis capabilities were incorporated into a finite element software system, called COSMOS/M. The range of the polynomial orders can be varied from 1 to 10 for two dimensional linear elastic analysis. In the h-p adaptive analysis process, a refined mesh are first achieved via adaptive h-refinement. The p-refinement is then added on to the h-version designed mesh by uniformly increasing the degree of the polynomials. Some numerical results computed by COSMOS/M are presented to illustrate the performance of these p and h-p analysis capabilities.