In this paper, we study the stationary Black-Scholes model arising in finance with transaction costs. This model becomes interesting when the time does not play a role such as, for instance, in perpetual options. The equation describing this model is a nonlinear second-order boundary value problem and there is no analytic solutions in closed form for such a nonlinear equation. After discretization via the centered finite difference formula we have to solve a nonlinear algebraic system which would be a serious problem when we use a small discretization mesh. We solve this nonlinear system by the residual-based Broyden's method, which is an efficient quasi-Newton method and is convenient to implement by a desk computer. We give a convergence analysis of the Broyden's method by assuming a lower and upper bound of the converged solution of the Black-Scholes model. Numerical results are given to show that the convergence rate of the method is robust with respect to the discretization mesh and the problem parameters.