Trajeto diário do Sol através do ano 

Trace a posição do Sol para a sua localização entre o nascer e o pôr do sol para os solstícios e o equinócio.

mostre o input completo de Wolfram Languageoculte o input

In[1]:=

X summersolstice = DateObject[{2013, 6, 20}]; autumnalequinox = DateObject[{2013, 9, 22}]; wintersolstice = DateObject[{2013, 12, 21}]; coordsToCartesian = Evaluate[CoordinateTransformData["Spherical" -> "Cartesian", "Mapping", {1, \[Pi]/2 - (#2 Degree), 2 Pi - (#1 Degree)}]] &; getValues = Function[{series}, Map[QuantityMagnitude, series["Values"], {2}]]; getSunriseToSunset = DateRange[Sunrise[#], Sunset[#], {30, "Minute"}] &; latitude = QuantityMagnitude[Latitude[$GeoLocation], "Radians"]; Graphics3D[{ Black, Table[ Text[{"W", "S", "E", "N"}[[j]], 1.1 {Cos[Pi/2 j], Sin[Pi/2 j], 0}], {j, 4}], Rotate[{ Arrow[{{0, 0, 0}, {0, 0, 1.4}}], Arrowheads[.02], Arrow@Table[{.1 Cos[t], .1 Sin[t], 1}, {t, 0, -3 Pi/2, -.2}]}, latitude, {0, 1, 0}], , Thickness[.007], Transpose[{{RGBColor[.49, 0, 0], RGBColor[.88, .63, .23], RGBColor[.27, .2, .75]}, (Line@(coordsToCartesian @@@ (getValues@ SunPosition[ getSunriseToSunset[#]]))) & /@ {summersolstice, autumnalequinox, wintersolstice}}], RGBColor[.6, .73, .36], Cylinder[{{0, 0, -.005}, {0, 0, 0}}], Opacity[.1], Sphere[]}, PlotRange -> {0, 1.4}, Boxed -> False, SphericalRegion -> True, ViewAngle -> \[Pi]/12, ViewCenter -> {.5, .5, .2}]

Out[1]=