DOC: Update randn() to use rng.normal()

Updates all references to `randn` to use `normal` from Generator.
I left the matlib versions alone.
I removed all references to `seed=123` and left  `# may vary`.
Longer outputs were reformatted to use Jupyter style output.

[skip actions] [skip azp] [skip cirrus]

Author: bmwoodruff

numpy/polynomial/hermite.py

@@ -1461,7 +1461,7 @@ def hermfit(x, y, deg, rcond=None, full=False, w=None):
     >>> from numpy.polynomial.hermite import hermfit, hermval
     >>> x = np.linspace(-10, 10)
     >>> rng = np.random.default_rng()
-    >>> err = rng.normal(scale=1./10,size=len(x))
+    >>> err = rng.normal(scale=1./10, size=len(x))
     >>> y = hermval(x, [1, 2, 3]) + err
     >>> hermfit(x, y, 2)
     array([1.02294967, 2.00016403, 2.99994614]) # may vary

numpy/polynomial/hermite_e.py

@@ -1388,7 +1388,7 @@ def hermefit(x, y, deg, rcond=None, full=False, w=None):
     >>> from numpy.polynomial.hermite_e import hermefit, hermeval
     >>> x = np.linspace(-10, 10)
     >>> rng = np.random.default_rng()
-    >>> err = rng.normal(scale=1./10,size=len(x))
+    >>> err = rng.normal(scale=1./10, size=len(x))
     >>> y = hermeval(x, [1, 2, 3]) + err
     >>> hermefit(x, y, 2)
     array([1.02284196, 2.00032805, 2.99978457]) # may vary

numpy/polynomial/laguerre.py

@@ -1442,7 +1442,7 @@ def lagfit(x, y, deg, rcond=None, full=False, w=None):
     >>> from numpy.polynomial.laguerre import lagfit, lagval
     >>> x = np.linspace(0, 10)
     >>> rng = np.random.default_rng()
-    >>> err = rng.normal(scale=1./10,size=len(x))
+    >>> err = rng.normal(scale=1./10, size=len(x))
     >>> y = lagval(x, [1, 2, 3]) + err
     >>> lagfit(x, y, 2)
     array([1.00578369, 1.99417356, 2.99827656]) # may vary