Concentration of a diffusing substance in a medium was derived in various cases of uni-dimensional diffusion, including in a semi-infinite medium and a plate-shaped medium. Multi-dimensional diffusion involves boundary conditions in each coordinate direction. The algorithm dealing with uni-dimensional case becomes very complicated in multi-dimensional cases. This study proposes an algorithm, which is called the complementary method, that combines complementary functions of the normalized solution in uni-dimensional diffusion case by multiplication to solve those in various multi-dimensional diffusion cases with dramatically simplified mathematics. Besides, the complementary method is used to solve various kinds of boundary conditions for multi-dimensional diffusion.