Optimization of b-values in DWI-IVIM is necessary for precision and reproducibility of parameter estimation. The IVIM biexponential model with Total variation correction (BE+TV) method was used here, which reduces sudden changes in parameter values and removes non-physiological heterogeneity in maps. Various b-value combinations were utilized in prostate: 4, 8 and 13 b-values. 8 and 13 b-values showed similar trend in coefficient of variation (CV), which was confirmed using Bland-Altman plots, including comparable qualitative maps; in contrast, 4 b-values showed variable CV and noisy maps. In prostate, 8 b-values can be used with Total Variation method for faster acquisition protocol and low variability in parameter maps.