def eccentricity(jd):
 T = ((jd - jd1950) / 36525.0)
 p = ((-1.26e-07), (-4.193e-05), 0.01673011)
 return np.polyval(p, T)

def eccentricity ( jd ) : T = ( ( jd - jd1950 ) / 36525.0 ) p = ( ( - 1.26e-07 ) , ( - 4.193e-05 ) , 0.01673011 ) return np . polyval ( p , T )
eccentricity of the earths orbit at the requested julian date .