We present the results of FUV dust scattering simulation, which is based on the Monte-Carlo method. In this simulation, we focus on the multiple scattering in the Ophiuchus complex region because the single scattering case in the region already reported by Lee et al. 2008. We compare the simulation result to the FUV intensity with FIMS and the single scattering result. We also discuss the parameters related to the results of this simulation, such as asymmetry factor, albedo and other different setting-ups.