用de405星历在AU中的地日距离

Question

在一个项目中，我必须用de405星历来计算地日距离。我使用了pyephem，jplephem，skyfield也安装了de405包，但无法实现并编写了以下代码：

from skyfield.jpllib import Ephemeris
import de405

def esundist(d_year, d_month, d_day):
    eph = Ephemeris(de405)
    astro = eph.earth(utc=(d_year, d_month, d_year)).observe(eph.sun)
    ra, dec, distance = astro.radec
    print distance.AU
esundist(2005, 10, 21)

但这段代码给出了如下错误：

^{pr2}$

我只想计算地球到太阳的距离。如果可能，在pyephem中使用de405？在