https://stats.stackexchange.com/q/592473/69508

592473_69508.py

import matplotlib.pyplot as plt
import numpy as np
from scipy.stats import describe

if __name__ == '__main__':
    x = np.random.normal(2, 3, size=10**4)
    z = x + x
    q = 2 * x
    plt.hist(x, bins=100, alpha=0.1)
    plt.hist(z, bins=100, alpha=0.1)
    plt.hist(q, bins=100, alpha=0.1)
    plt.show()
    print('x', describe(x))
    print('z', describe(z))
    print('q', describe(q))