\documentclass[svgnames,a4paper]{standalone} \usepackage[T1]{fontenc} \usepackage[frenchb]{babel} \usepackage[garamond]{mathdesign} \usepackage{pstricks-add} \usepackage{pst-caelum} \usepackage{xint-meeus} \usepackage{readarray} \input{xint-sun.tex} \readdef{lune01.txt}\lignes \readarray{\lignes}\lune[120,2] \JulianDay(2021,6,10,11,0,58) \newcommand\calcTfor[1]{\xintdeffloatvar T := (JD+#1-2451545.0)/36525;} \calcTfor{0}% \def\RPSun{\FloatToFix{(90-xmDS(T))/20}} \def\APSun{\FloatToFix{-90-xmADS(T)}} \begin{document} \begin{pspicture}(-7,-7)(7,7) \caelum{Graticule}% \caelum[linewidth=0.2pt,linecolor=DimGray,linestyle=dotted]{Boundaries} \caelum[linecolor=Crimson]{Ecliptic}% \psdots[linecolor=Crimson,dotsize=4pt](\RPSun;\APSun)% \multido{\i=1+6}{20}{% \def\dec{\lune[\i,2]} \def\asc{\lune[\i,1]} \def\RPMoon{\FloatToFix{(90-\dec)/20}} \def\APMoon{\FloatToFix{-90-\asc}} \psdots[linecolor=DimGray,dotsize=4pt](\RPMoon;\APMoon) } \end{pspicture} \end{document}