I am using python's package ephem.
I can not locate function to get tilt angle between ecliptic plane and equatorial plane.
Hope some one can help me.
Since 0° right ascension is defined as the point where the equatorial and ecliptic planes cross, they are farthest apart at 90° right ascension, where their distance from each other is the same as the Earth’s tilt on its axis:
from ephem import Equatorial, Ecliptic
p = Equatorial('90', '0', epoch='2000')
print Ecliptic(p).lat
The result is the Earth’s axial tilt at the beginning of the year 2000:
23:26:21.4
Bit of Googling might give you your answer: