#Nuray Kutlu
#April 14
#OK to post

#since M/V is the density of an object, if we show that (Mass of sun / Volume of sun) 
#< (Mass of Moon / Volume of Moon), we would have shown that the density of the moon is 
#greater than the mass of the sun. 
# We can say that M_moon/(D_moon )^3 is greater than M_sun/(D_sun)^3, due to that fact 
#that tidal force is proportional to M/D_object_from_earth^3, and the tidal force 
#of the moon is far greater than that of the sun. 
#We showed that Sun_distance/Sun_radius = Moon_distance/Moon_radius in class. 
#This also means Sun_distance^3/Sun_radius^3 = Moon_distance^3/Moon_Radius^3.
#Since the two values are equal, we can multiply each one on both sides
#and maintaining the inequality:
#( M_moon/(D_moon)^3 )*(Moon_distance^3/Moon_Radius^3)>( M_sun/(D_sun)^3)*(Sun_distance^3/Sun_radius^3)
#Simplifying this, we can get M_moon/Moon_Radius^3 > M_sun/Sun_radius^3
#Since the  volume of a sphere is directly proportionate to it's radius cubed,
#we can say that M_moon/V_moon > M_sum/V_sun and so, D_moon>D_sun


Help:=proc(): print(`SE()`):
end:

SE:=proc()
local P,t,SunCenter,SunRadius,MoonCenter,MoonRadius,EarthCenter,EarthRadius:
SunCenter:=[0,0]:
SunRadius:=3:
MoonCenter:=[5.9,5.9]:
MoonRadius:=.3:
EarthCenter:=[7,7]:
EarthRadius:=1:
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):
P:=P,plot([EarthCenter[1]+EarthRadius*cos(t),EarthCenter[2]+EarthRadius*sin(t),
t=0..2*Pi],axes=none,filled=true,color=purple):
P:=P,plot([EarthCenter[1]+.3+EarthRadius*5*cos(t),EarthCenter[2]+.3+EarthRadius*1.5*sin(t),
t=0..2*Pi],axes=none,filled=false,color=green):
P:=P,plot([EarthCenter[1]-.18+EarthRadius*.6*cos(t),EarthCenter[2]-.18+EarthRadius*.6*sin(t),
t=0..2*Pi],axes=none,filled=true,color=lavender):
P:=P,plot([SunCenter[1]+.75+SunRadius*3*cos(t),SunCenter[2]+.75+SunRadius*3*sin(t),
t=0..2*Pi],axes=none,filled=false,color=purple):
display(P):
end: