mean1 = data[:,0].mean () mean2 = data[:,1].mean () std1 = data[:,0].std () std2 = data[:,1].std () params = np.array ( [(mean1 - 0.2, mean2 - 1, std1, std2, 0), (mean1 + 0.2, mean2 + 1, std1, std2, 0)] ) weights = np.array ( [ 0.5, 0.5 ] )