Haga cálculos de tiempo celestial 

Trace la ecuación del tiempo, la diferencia entre la ascensión recta del Sol al mediodía y la hora sideral al mediodía.

In[1]:=

X sunpos = SunPosition[GeoPosition[{0, 0}], DateRange[DateObject[{2014, 1, 1, 12, 0}, TimeZone -> 0], DateObject[{2014, 12, 31, 12, 0}, TimeZone -> 0], 10], CelestialSystem -> "Equatorial"]; stime = SiderealTime[GeoPosition[{0, 0}], DateRange[DateObject[{2014, 1, 1, 12, 0}, TimeZone -> 0], DateObject[{2014, 12, 31, 12, 0}, TimeZone -> 0], 10]]; equationoftime = TimeSeriesThread[ With[{diff = First[#][[1]] - Last[#]}, UnitConvert[ Mod[diff, Quantity[24, "HoursOfRightAscension"], Quantity[-12, "HoursOfRightAscension"]], "MinutesOfRightAscension"]] &, {sunpos, stime}]; DateListPlot[equationoftime, PlotLabel -> "Equation of Time", Axes -> True]

Out[1]=