Automatic segmentation of multiple sclerosis (MS) lesions in brain magnetic resonance imaging (MRI) has been widely investigated in the recent years with the goal of helping MS diagnosis and patient follow-up. In this research work, Gaussian mixture model (GMM) has been used to segment the MS lesions in MRIs, including T1-weighted (T1-w), T2-w, and T2-fluid attenuation inversion recovery. Usually, GMM is optimized by using expectation-maximization (EM) algorithm. The drawbacks of this optimization method are, it does not converge to optimal maximum or minimum and furthermore, there are some voxels, which do not fit the GMM model and have to be rejected. So, GMM is time-consuming and not too much efficient. To overcome these limitations, in this research study, at the first step, GMM was applied to segment only T1-w images by using 100 various starting points when the maximum number of iterations was considered to be 50. Then segmentation results were used to calculate the parameters of the other two images. Furthermore, FAST-trimmed likelihood estimator algorithm was applied to determine which voxels should be rejected. The output result of the segmentation was classified in three classes; White and Gray matters, cerebrospinal fluid, and some rejected voxels which prone to be MS. In the next phase, MS lesions were detected by using some heuristic rules. This new method was applied on the brain MRIs of 25 patients from two hospitals. The automatic segmentation outputs were scored by two specialists and the results show that our method has the capability to segment the MS lesions with dice similarity coefficient score of 0.82. The results showed a better performance for the proposed approach, in comparison to those of previous works with less time-consuming.