Copy here speed from JPL Horizons ephemeris: dRA*cosD d(DEC)/dt OR dRA*cosD, d(DEC)/dt line