def plot_gaussian_multiply(): xs = np.arange(-5, 10, 0.1) mean1, var1 = 0, 5 mean2, var2 = 5, 1 mean, var = stats.mul(mean1, var1, mean2, var2) ys = [stats.gaussian(x, mean1, var1) for x in xs] plt.plot(xs, ys, label='M1') ys = [stats.gaussian(x, mean2, var2) for x in xs] plt.plot(xs, ys, label='M2') ys = [stats.gaussian(x, mean, var) for x in xs] plt.plot(xs, ys, label='M1 x M2') plt.legend() plt.show()