#C22.txt, April 8, 2024, Solar Eclipse
Help:=proc(): print(`ProveEclipse()`):
print(`SE(SunCenter,SunRadius,MoonCenter,MoonRadius,EarthCenter,EarthRadius)`):

end:
#km
#ProveEclipse(): outputs the apparent sizes of the Sun and Moon
ProveEclipse:=proc(): local SunRad,SunDis,MoonRad,MoonDis:
SunRad:=696.*10^3;
SunDis:=149.6*10^6;
MoonRad:=1737.;
MoonDis:=384.*10^3;
[SunRad/SunDis, MoonRad/MoonDis]:
end:

with(plots):

SE:=proc(SunCenter,SunRadius,MoonCenter,MoonRadius,EarthCenter,EarthRadius)
local P,t:
P:=plot([SunCenter[1]+SunRadius*cos(t),SunCenter[2]+SunRadius*sin(t),
t=0..2*Pi],axes=none,filled=true,color=yellow):
P:=P,plot([MoonCenter[1]+MoonRadius*cos(t),MoonCenter[2]+MoonRadius*sin(t),
t=0..2*Pi],axes=none,filled=true,color=blue):
display(P):
end: