We present a Monte-Carlo simulation code, which solves the problem of dust-scattering in interstellar dust clouds with arbitrary light source distribution and dust density structure, and calculate the surface brightness distribution. The method is very flexible and can be applied to radiative transfer problems occurring not only in a single dust cloud, but also in extragalactic dust environment. We compare, for performance test, the result of Monte-Carlo simulation with the well-known analytic approximation for a spherically symmetric homogeneous cloud. We find that the Code approximation gives a very accurate result.