diff --git a/doc/source/reference/routines.polynomials.rst b/doc/source/reference/routines.polynomials.rst
index 0763a1cf719a..0da51e57a308 100644
--- a/doc/source/reference/routines.polynomials.rst
+++ b/doc/source/reference/routines.polynomials.rst
@@ -47,23 +47,27 @@ The `~numpy.polynomial.polynomial.Polynomial` class is imported for brevity::
     from numpy.polynomial import Polynomial
 
 
-+------------------------+------------------------------+---------------------------------------+
-|  **How to...**         | Legacy (`numpy.poly1d`)      | `numpy.polynomial`                    |
-+------------------------+------------------------------+---------------------------------------+
-| Create a               | ``p = np.poly1d([1, 2, 3])`` | ``p = Polynomial([3, 2, 1])``         |
-| polynomial object      |                              |                                       |
-| from coefficients [1]_ |                              |                                       |
-+------------------------+------------------------------+---------------------------------------+
-| Create a polynomial    | ``r = np.poly([-1, 1])``     | ``p = Polynomial.fromroots([-1, 1])`` |
-| object from roots      | ``p = np.poly1d(r)``         |                                       |
-+------------------------+------------------------------+---------------------------------------+
-| Fit a polynomial of    |                              |                                       |
-| degree ``deg`` to data | ``np.polyfit(x, y, deg)``    | ``Polynomial.fit(x, y, deg)``         |
-+------------------------+------------------------------+---------------------------------------+
-
++------------------------+----------------------------------------+---------------------------------------+
+|  **How to...**         | Legacy (`numpy.poly1d`)                | `numpy.polynomial`                    |
++------------------------+----------------------------------------+---------------------------------------+
+| Create a               | ``p = np.poly1d([1, 2, 3])``           | ``p = Polynomial([3, 2, 1])``         |
+| polynomial object      |                                        |                                       |
+| from coefficients [1]_ |                                        |                                       |
++------------------------+----------------------------------------+---------------------------------------+
+| Create a polynomial    | ``r = np.poly([-1, 1])``               | ``p = Polynomial.fromroots([-1, 1])`` |
+| object from roots      | ``p = np.poly1d(r)``                   |                                       |
++------------------------+----------------------------------------+---------------------------------------+
+| Fit a polynomial of    |                                        |                                       |
+| degree ``deg`` to data | ``np.polyfit(x, y, deg)``              | ``Polynomial.fit(x, y, deg)``         |
++------------------------+----------------------------------------+---------------------------------------+
+| Evaluate a polynomial  | ``p(2.0)`` or                          | ``p(2.0)`` or ``polyval(2.0, p.coef)``|
+| at a point [2]_        | ``np.polyval([1, 2, 3], 2.0)``         | (use ``p.convert().coef`` after fit)  |
++------------------------+----------------------------------------+---------------------------------------+
 
 .. [1] Note the reversed ordering of the coefficients
 
+.. [2] When evaluating polynomials created with ``fit()``, use ``p(x)`` or ``polyval(x, p.convert().coef)`` to handle domain/window scaling correctly.
+
 Transition Guide
 ~~~~~~~~~~~~~~~~
 
@@ -187,4 +191,4 @@ Documentation for legacy polynomials
 .. toctree::
    :maxdepth: 2
 
-   routines.polynomials.poly1d
+   routines.polynomials.poly1d
\ No newline at end of file
diff --git a/numpy/polynomial/polynomial.py b/numpy/polynomial/polynomial.py
index 6ec0dc58a1de..d9427b199fea 100644
--- a/numpy/polynomial/polynomial.py
+++ b/numpy/polynomial/polynomial.py
@@ -717,6 +717,10 @@ def polyval(x, c, tensor=True):
     -----
     The evaluation uses Horner's method.
 
+   When using coefficients from polynomials created with ``Polynomial.fit()``,
+   use ``p(x)`` or ``polyval(x, p.convert().coef)`` to handle domain/window
+   scaling correctly, not ``polyval(x, p.coef)``.
+
     Examples
     --------
     >>> import numpy as np