The classical model of the universe

> with(plots):with(plottools):

> equator := spacecurve([cos(t),sin(t),0,t=0..2*Pi],
thickness=2,color=blue):

> ecliptic := rotate(spacecurve([cos(t),sin(t),0 ,t=0..2*Pi],color=red) ,23*Pi/180,0,0 ):

> zodiac := rotate( display({seq(hexahedron([cos(T/12*2*Pi),
sin(T/12*2*Pi),0], .02,color=black),T=0..11)}),23*Pi/180,0,0):

> celeste := sphere([0,0,0],1,style=contour):

> earth := sphere([0,0,0],.1,style=point):

> sun := day -> rotate(sphere([cos(-day*2*Pi/360),sin(-day*2*Pi/360),0],.05,
style=point,color=yellow),23.5*Pi/180,0,0):

> universe := display([equator,ecliptic,zodiac,celeste,earth]):

> current := (day,hour)-> display([rotate(sun(day),0,0,hour*2*Pi/24),
rotate(universe,0,0,
hour*2*Pi/24)],orientation=[0,81]):

> current(130,14):

> display([seq(current(i*15,0),i=0..23)],insequence=true);