In [4]:
import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as stats
In [59]:
a = 1
L = 5
N = np.sqrt(2 * np.sqrt(8 * a ** 3 / np.pi))
x = np.linspace(-L, L, 1e5)
y = N ** 2 * x * x * np.exp(-x ** 2 * a)
In [62]:
plt.plot(x, y)
plt.show()
In [64]:
np.var(x)
Out[64]:
In [ ]: