In this paper, finite element method is applied to Leland's model for numerical simulation of option pricing with transaction costs. Spatial finite element models based on P1 and/or P2 elements are formulated in combination with a Crank-Nicolson-type temporal scheme. The temporal scheme is implemented using the Rannacher approach. Examples with several sets of parameter values are presented and compared with finite difference results in the literature. Spatial-temporal mesh-size ratios are observed for controlling the stability of our method. Our results compare favorably with the finite difference results in the literature for the model.