def METHOD_NAME(jd):
    """
    Eccentricity of the Earth's orbit at the requested Julian Date.
    Parameters
    ----------
    jd : scalar or array-like
        Julian date at which to compute the eccentricity
    Returns
    -------
    eccentricity : scalar or array
        The eccentricity (or array of eccentricities)
    References
    ----------
    * Explanatory Supplement to the Astronomical Almanac: P. Kenneth
      Seidelmann (ed), University Science Books (1992).
    """
    T = (jd - jd1950) / 36525.0
    p = (-0.000000126, -0.00004193, 0.01673011)
    return np.polyval(p, T)