Monte Carlo simulations were used to improve the accuracy of light distribution in tissues without increasing computation time. By analytically considering source attenuation separately from the scattering process, the fluence rate was calculated more accurately. Using a non-linear grid helped reduce computation time while maintaining high resolution in important areas. The theory of source attenuation function was discussed, and results showed the validity of the analytic function and the effectiveness of non-linear grids compared to standard techniques.