Realize cálculos com tempo astronômico 

Faça o gráfico da equação de tempo—a diferença entre a ascensão reta do Sol ao meio-dia e a hora sideral ao meio-dia.

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]=