The two-sided nonlinear boundary crossing probabilities for one-dimensional Brownian motion and related processes have been studied in Fu and Wu (2010) based on the finite Markov chain imbedding technique. It provides an efficient numerical method to computing the boundary crossing probabilities. In this paper we extend the above results for high-dimensional Brownian motion. In particular, we obtain the rate of convergence for high-dimensional boundary crossing probabilities. Numerical results are also provided to illustrate our results.