The problem of pricing European options based on multiple assets with transaction costs is considered. These options include, for example, quality options and options on the minimum of two or more risky assets. The value of these options is the solution of a nonlinear parabolic partial differential equation subject to a final condition given by the payoff function associated with the option. A computationally efficient method to solve this final-value problem is proposed. This method is based on an asymptotic expansion of the required solution with respect to the parameters related to the transaction costs followed by the numerical solution of the linear partial differential equations obtained at each order in perturbation theory. The numerical solution of these linear problems involves an implicit finite-difference scheme for the parabolic equation and the use of the fast Fourier sine transform to solve the resulting elliptic problems. Numerical results obtained on test problems with the method proposed here are shown and discussed.