Abstract: A direct method of using polynomial terms as stress function and combining with Matlab program design is introduced for solving problems of elasticity. First the present method has very strong rule for choice stress function, on the other hand, it simplifies the process of determining unknown coefficients. The method improves disadvantages of traditional inverse method or semi-inverse method that is necessary that choice stress function through repeated tasting and modifies for obtaining satisfactory results. Meanwhile we denote that some problems doesn't possess polynomial solution. The present approach is successful conjunction example of between modern computer technology and traditional elasticity analytical method. It pump new activity into elasticity analytical method.