A new method is suggested for the solution of plane elasticity problems. With use of the dislocation model in the crack problems, the basic scheme of this method is to find equilibrium Burgers vectors of dislocations which are distributed along the boundary of the first fundamental boundary value problems. The stress distribution in the region can be found by superposition of the contributions of each dislocation. The method is applied to three cases with known analytical solutions, and to a V-notched specimen under uniaxial tension. The numerical results are compared with other available solutions. This method is effective and simple in its use, compared with other numerical methods. The method also provides very accurate solutions in the region except near the boundary where the discretization error is significant. The extrapolation method is suggested for the stresses in the boundary region. Extensive application are also suggested for a general estimate of the computational efficiency of the method.