In the nonlinearly numerical simulation of finite element for materials, components and structures, during the simulation it is in general impossible for the load increment control method to simulate the peak point and the softening decent segment of special equilibrium path for load-displacement relation. In this paper, a new theory of displacement control proposed by Du Xiuli et al. was used to overcome this problem. To make a new stiffness matrix to be symmetric and banded, the main diagonal elements of stiffness matrix in the conventional finite element incremental equation were multiplied by a huge number, and the boundary conditions were also introduced to the finite element equation. Therefore, only a little modification is needed for the program of traditional load control method which can simplify the iterative process of calculation. Finally, it was programmed by MATLAB. As applications, some practically numerical examples were given. It was shown by the analysis that the method presented in the paper is not only effective, feasible, but also simple.