Background. Monte Carlo determination of TG-43 brachytherapy dosimetry parameters and dose distribution calculation for (131)CS source model CS-1 are presented in this study. Materials and Methods: The dose distribution was calculated around the (131)CS Model CS-1 located in the center of 30 cm x30 cm x30 cm water, and soft tissue phantoms cube using MCNP code by Monte Carlo method. The percentage depth dose (PDD) variation along the different axis, parallel and perpendicular, the source was calculated. Then, the isodose curves for 100%, 75%, 50% and 25% PDD were constructed. Finally, F(r,theta) and g(r) dosimetry parameters of TG-43 protocol have been determined. Results: Results obtained show that the Monte Carlo method could only calculate dose deposition in high gradient region, near the. source, accurately. The energy cut off was found to be 1 eV and the error in the calculations was less than 2%. Conclusion: The isodose curves of the CS-1 (131)CS source were constructed from dose calculation by MCNP code. The calculated dosimetry parameters for the source were in agreement with previously published results.