Move Earth

Move Sun

81.6527

Sun-Earth

Soleil1=Cercle[S, Rs]

inclinaison

Dplc1=((theta, 0))

Aphélie

Périhélie

ArcTrajectoire1_2=Arc[Trajectoire_T, T1, T4]

176.0576

Date9=2455927.50028 + t

Date2=floor((floor(Date9 + 0.5) - 1867216.25) / 36524.25)

Date1=floor(Date9 + 0.5) + 1 + Date2 - floor(Date2 / 4)

Date3=Date1 + 1524

Date4=floor((Date3 - 122.1) / 365.25)

Date5=floor(365.25Date4)

Date7=(Date9 + 0.5 - floor(Date9 + 0.5)) 24

Date8=floor(Date7)

Date13=floor((Date7 - Date8) (60))

Date10=Reste[floor(Date9 - 0.5) + 0.5 + 1.5, 7]

Date6=floor((Date3 - Date5) / 30.6001)

Date11=Date3 - Date5 - floor(30.6001Date6) + Reste[Date9 + 0.5, 1]

Date12=Date6 - Si[Date6 < 13.5, 1, 13]

Date14=Date4 - Si[Date12 > 2.5, 4716, 4715]

June solstice

September equinox

December solstice

March equinox

Laurent.Costa@ac-lyon.fr

excentricité

annual period

Périhélie

y

ua

Aphélie

Paramètre

x

ArcTrajectoire2_1=Arcd'éllipse[Trajectoire_T, T4, T1]

ArcTrajectoire1_2=Arc[Trajectoire_T, T1, T4]

27.16-27.32158

x

y

angle/écliptique

z

Earth-Moon

R_E

-0.0783 t + 76.923

Lune1=Cercle[L, Rl]

P_{NORD1}=Point3D[Dplc, -1.3 Rt cos(66.5°), 0, 1.3 Rt sin(66.5°)] + T

P_{SUD1}=Point3D[Dplc, 1.3 Rt cos(66.5°), 0, -1.3 Rt sin(66.5°)] + T

Equateur2=Conique[Point3D[Dplc, 0, -Rt cos(23.5°), 0Rt sin(23.5°)] + T, Point3D[Dplc, Rt cos(23.5°), 0, Rt sin(23.5°)] + T, Point3D[Dplc, Rt cos(°) cos(23.5°), Rt sin(°) cos(23.5°), Rt sin(23.5°)] + T, Point3D[Dplc, 0, Rt cos(23.5°), 0Rt sin(23.5°)] + T, Point3D[Dplc, -Rt cos(23.5°), 0, -Rt sin(23.5°)] + T]

axeN2=Segment[P_{NORD}, P_{NORD1}]

AxeS2=Segment[P_{SUD1}, P_{SUD}]

P_{Cancer2}=Point3D[Dplc, Rt cos(47°), 0, Rt sin(47°)] + T

Cancer2=Conique[Point3D[Dplc, 0, -cos(23.5°) Rt cos(23.5°), 0] + EqC1, Point3D[Dplc, cos(23.5°) Rt cos(23.5°), 0, cos(23.5°) Rt sin(23.5°)] + EqC1, Point3D[Dplc, cos(23.5°) Rt cos(°) cos(23.5°), cos(23.5°) Rt sin(°) cos(23.5°), cos(23.5°) Rt sin(23.5°)] + EqC1, Point3D[Dplc, 0, cos(23.5°) Rt cos(23.5°), 0Rt sin(23.5°)] + EqC1, Point3D[Dplc, -cos(23.5°) Rt cos(23.5°), 0, -cos(23.5°) Rt sin(23.5°)] + EqC1]

phi ≤ 3.14159 / 2

N

S

N

S

Lune3=Cercle[L, Rl]

phi > 3.14159 / 2

EstEntier[floor(10(t + 180) / 180) / 2]

phi ≤ 3.14159 / 2

phi ≤ 3.14159 / 2

Programmation : Laurent COSTA 
Lycée professionnel Emile Bejuit 69500 Bron
Sous la direction de Charles-Henri Eyraud

F18 ≟ 0° ∧ (phi ≥ 3.14159 / 2) ∧ (phi < 3.14159)

F18 ≟ 0.5° ∧ (phi ≥ 3.14159 / 2) ∧ (phi < 3.14159)

PS_{NORD1}=Point3D[Dplc, 0, 0, Rs] + S

axesoleil22=Segment[PS_{SUD1}, PS_{SUD}]

A_3=Point3D[Dplc, Rt cos((20 (t + 10))°) cos(23.5°), Rt sin((20 (t + 10))°) cos(23.5°), Rt cos((20 (t + 10))°) sin(23.5°)] + T

A_{10}=Point3D[Dplc, Rt cos((10 (t + 90))°) cos(23.5°), Rt sin((10 (t + 90))°) cos(23.5°), Rt cos((10 (t + 90))°) sin(23.5°)] + T

charles-henri.eyraud@ens-lyon.fr

<--Nouvelle Lune  ** ajustement_L-->

<--Pleine lune   ** ajustement_L-->

D + (3, -2) - ((x(F) - x(D))⁵, (y(F) - y(D))⁵)

Sud

Nord

Segment[Rotation[Eq_{fixe3}, 23.5°, Tfixe], Symétrie[Rotation[Eq_{fixe3}, 23.5°, Tfixe], l_2]]

N

N

S

S

Rayon2=Tfixe + (-Rt2 - 1, Rt2/4)

Eq_{fixe1}+abs(sin(t)) (Eq_{fixe2}-Eq_{fixe1})

lune1=Point[arc_{CacheL}]

Cliquer sur        pour relancer l'animation

tmax =

G = (-4.9654811232033795, -2.386409588617408)
E = (-6.476880240789697, -3.5129182476879506)
F = (-8.2, 0.58006)
D = (-5.969951344207951, -3.5129182476879506)
u = 1.3 (slider from 0.6 to 5, step 0.1)
R69 = 0
R68 = 0
R67 = 0
R66 = 0
R65 = 0
R64 = 0
R63 = 0
R62 = 0
R61 = 0
R60 = 0
R59 = 0
R58 = 0
R57 = 0
R56 = 0
R55 = 0
R54 = 0
R53 = 0
R52 = 0
R51 = 0
R1 = 76.5
R50 = 22
R49 = 19
R48 = 27
R47 = 27
R46 = 25
R45 = 22
R44 = 39
R43 = 30
R42 = 34
R41 = 28
R40 = 34
R39 = 31
R38 = 34
R37 = 35
R36 = 35
R35 = 34
R34 = 33
R33 = 43
R32 = 35
R31 = 45
R30 = 38
R29 = 49
R28 = 42
R27 = 50
R26 = 47
R25 = 53
R24 = 53
R23 = 53
R22 = 58
R21 = 55
R20 = 62.5
R19 = 55
R18 = 65
R17 = 57
R16 = 66.5
R15 = 58
R14 = 67
R13 = 61
R12 = 66
R11 = 64.5
R10 = 67
R9 = 69
R8 = 67
R7 = 73.5
R6 = 68
R5 = 78
R4 = 70.5
R3 = 81.5
R2 = 73.5
P31 = 435.871
P49 = 702.057
P48 = 686.678
P47 = 672.057
P46 = 657.5505
P45 = 643.107
P44 = 627.068
P43 = 613.107
P42 = 598.158
P41 = 583.994
P40 = 568.845
P39 = 554.3862
P38 = 539.565
P37 = 524.7495
P36 = 510.105
P35 = 495.485
P34 = 480.91666
P33 = 465.485
P32 = 451.438
P30 = 421.895
P29 = 406.349
P28 = 392.2
P27 = 376.865
P26 = 362.4739
P25 = 347.4045
P18 = 243.664
P14 = 184.869
P10 = 126.234
P24 = 332.658
P23 = 317.964
P22 = 302.8685
P21 = 288.585
P20 = 273.221
P19 = 259.174
P17 = 229.746
P16 = 214.2277
P15 = 200.267
P13 = 170.711
P12 = 155.55
P11 = 141.075
P9 = 111.389
P8 = 96.89
P7 = 81.6527
P6 = 67.4458
P5 = 51.9846
P4 = 37.9561
P3 = 22.362
P2 = 8.3557
P1 = -6.548
P50 = 716.437
P51 = 731.5
P52 = 0
P53 = 0
P54 = 0
P55 = 0
P56 = 0
P57 = 0
P58 = 0
P59 = 0
P60 = 0
P61 = 0
P62 = 0
P63 = 0
P64 = 0
P65 = 0
P66 = 0
P67 = 0
P68 = 0
t = 407.486 (slider from 0 to Q100, step 1.0E-7)
S68 = If(t ≤ P68, If[t > P67, (R68 - R67) * (t - P67) / (P68 - P67) + R67, 0], 0)
S67 = If(t ≤ P67, If[t > P66, (R67 - R66) * (t - P66) / (P67 - P66) + R66, 0], 0)
S66 = If(t ≤ P66, If[t > P65, (R66 - R65) * (t - P65) / (P66 - P65) + R65, 0], 0)
S65 = If(t ≤ P65, If[t > P64, (R65 - R64) * (t - P64) / (P65 - P64) + R64, 0], 0)
S64 = If(t ≤ P64, If[t > P63, (R64 - R63) * (t - P63) / (P64 - P63) + R63, 0], 0)
S63 = If(t ≤ P63, If[t > P62, (R63 - R62) * (t - P62) / (P63 - P62) + R62, 0], 0)
S62 = If(t ≤ P62, If[t > P61, (R62 - R61) * (t - P61) / (P62 - P61) + R61, 0], 0)
S61 = If(t ≤ P61, If[t > P60, (R61 - R60) * (t - P60) / (P61 - P60) + R60, 0], 0)
S60 = If(t ≤ P60, If[t > P59, (R60 - R59) * (t - P59) / (P60 - P59) + R59, 0], 0)
S59 = If(t ≤ P59, If[t > P58, (R59 - R58) * (t - P58) / (P59 - P58) + R58, 0], 0)
S58 = If(t ≤ P58, If[t > P57, (R58 - R57) * (t - P57) / (P58 - P57) + R57, 0], 0)
S57 = If(t ≤ P57, If[t > P56, (R57 - R56) * (t - P56) / (P57 - P56) + R56, 0], 0)
S56 = If(t ≤ P56, If[t > P55, (R56 - R55) * (t - P55) / (P56 - P55) + R55, 0], 0)
S55 = If(t ≤ P55, If[t > P54, (R55 - R54) * (t - P54) / (P55 - P54) + R54, 0], 0)
S54 = If(t ≤ P54, If[t > P53, (R54 - R53) * (t - P53) / (P54 - P53) + R53, 0], 0)
S53 = If(t ≤ P53, If[t > P52, (R53 - R52) * (t - P52) / (P53 - P52) + R52, 0], 0)
S52 = If(t ≤ P52, If[t > P51, (R52 - R51) * (t - P51) / (P52 - P51) + R51, 0], 0)
S51 = If(t ≤ P51, If[t > P50, (R51 - R50) * (t - P50) / (P51 - P50) + R50, 0], 0)
S50 = If(t ≤ P50, If[t > P49, (R50 - R49) * (t - P49) / (P50 - P49) + R49, 0], 0)
S49 = If(t ≤ P49, If[t > P48, (R49 - R48) * (t - P48) / (P49 - P48) + R48, 0], 0)
S48 = If(t ≤ P48, If[t > P47, (R48 - R47) * (t - P47) / (P48 - P47) + R47, 0], 0)
S47 = If(t ≤ P47, If[t > P46, (R47 - R46) * (t - P46) / (P47 - P46) + R46, 0], 0)
S46 = If(t ≤ P46, If[t > P45, (R46 - R45) * (t - P45) / (P46 - P45) + R45, 0], 0)
S45 = If(t ≤ P45, If[t > P44, (R45 - R44) * (t - P44) / (P45 - P44) + R44, 0], 0)
S44 = If(t ≤ P44, If[t > P43, (R44 - R43) * (t - P43) / (P44 - P43) + R43, 0], 0)
S43 = If(t ≤ P43, If[t > P42, (R43 - R42) * (t - P42) / (P43 - P42) + R42, 0], 0)
S42 = If(t ≤ P42, If[t > P41, (R42 - R41) * (t - P41) / (P42 - P41) + R41, 0], 0)
S41 = If(t ≤ P41, If[t > P40, (R41 - R40) * (t - P40) / (P41 - P40) + R40, 0], 0)
S40 = If(t ≤ P40, If[t > P39, (R40 - R39) * (t - P39) / (P40 - P39) + R39, 0], 0)
S39 = If(t ≤ P39, If[t > P38, (R39 - R38) * (t - P38) / (P39 - P38) + R38, 0], 0)
S38 = If(t ≤ P38, If[t > P37, (R38 - R37) * (t - P37) / (P38 - P37) + R37, 0], 0)
S37 = If(t ≤ P37, If[t > P36, (R37 - R36) * (t - P36) / (P37 - P36) + R36, 0], 0)
S36 = If(t ≤ P36, If[t > P35, (R36 - R35) * (t - P35) / (P36 - P35) + R35, 0], 0)
S35 = If(t ≤ P35, If[t > P34, (R35 - R34) * (t - P34) / (P35 - P34) + R34, 0], 0)
S34 = If(t ≤ P34, If[t > P33, (R34 - R33) * (t - P33) / (P34 - P33) + R33, 0], 0)
S33 = If(t ≤ P33, If[t > P32, (R33 - R32) * (t - P32) / (P33 - P32) + R32, 0], 0)
S32 = If(t ≤ P32, If[t > P31, (R32 - R31) * (t - P31) / (P32 - P31) + R31, 0], 0)
S31 = If(t ≤ P31, If[t > P30, (R31 - R30) * (t - P30) / (P31 - P30) + R30, 0], 0)
S30 = If(t ≤ P30, If[t > P29, (R30 - R29) * (t - P29) / (P30 - P29) + R29, 0], 0)
S29 = If(t ≤ P29, If[t > P28, (R29 - R28) * (t - P28) / (P29 - P28) + R28, 0], 0)
S28 = If(t ≤ P28, If[t > P27, (R28 - R27) * (t - P27) / (P28 - P27) + R27, 0], 0)
S27 = If(t ≤ P27, If[t > P26, (R27 - R26) * (t - P26) / (P27 - P26) + R26, 0], 0)
S26 = If(t ≤ P26, If[t > P25, (R26 - R25) * (t - P25) / (P26 - P25) + R25, 0], 0)
S25 = If(t ≤ P25, If[t > P24, (R25 - R24) * (t - P24) / (P25 - P24) + R24, 0], 0)
S24 = If(t ≤ P24, If[t > P23, (R24 - R23) * (t - P23) / (P24 - P23) + R23, 0], 0)
S23 = If(t ≤ P23, If[t > P22, (R23 - R22) * (t - P22) / (P23 - P22) + R22, 0], 0)
S22 = If(t ≤ P22, If[t > P21, (R22 - R21) * (t - P21) / (P22 - P21) + R21, 0], 0)
S21 = If(t ≤ P21, If[t > P20, (R21 - R20) * (t - P20) / (P21 - P20) + R20, 0], 0)
S20 = If(t ≤ P20, If[t > P19, (R20 - R19) * (t - P19) / (P20 - P19) + R19, 0], 0)
S19 = If(t ≤ P19, If[t > P18, (R19 - R18) * (t - P18) / (P19 - P18) + R18, 0], 0)
S18 = If(t ≤ P18, If[t > P17, (R18 - R17) * (t - P17) / (P18 - P17) + R17, 0], 0)
S17 = If(t ≤ P17, If[t > P16, (R17 - R16) * (t - P16) / (P17 - P16) + R16, 0], 0)
S16 = If(t ≤ P16, If[t > P15, (R16 - R15) * (t - P15) / (P16 - P15) + R15, 0], 0)
S15 = If(t ≤ P15, If[t > P14, (R15 - R14) * (t - P14) / (P15 - P14) + R14, 0], 0)
S14 = If(t ≤ P14, If[t > P13, (R14 - R13) * (t - P13) / (P14 - P13) + R13, 0], 0)
S13 = If(t ≤ P13, If[t > P12, (R13 - R12) * (t - P12) / (P13 - P12) + R12, 0], 0)
S12 = If(t ≤ P12, If[t > P11, (R12 - R11) * (t - P11) / (P12 - P11) + R11, 0], 0)
S11 = If(t ≤ P11, If[t > P10, (R11 - R10) * (t - P10) / (P11 - P10) + R10, 0], 0)
S10 = If(t ≤ P10, If[t > P9, (R10 - R9) * (t - P9) / (P10 - P9) + R9, 0], 0)
S9 = If(t ≤ P9, If[t > P8, (R9 - R8) * (t - P8) / (P9 - P8) + R8, 0], 0)
S8 = If(t ≤ P8, If[t > P7, (R8 - R7) * (t - P7) / (P8 - P7) + R7, 0], 0)
S7 = If(t ≤ P7, If[t > P6, (R7 - R6) * (t - P6) / (P7 - P6) + R6, 0], 0)
S6 = If(t ≤ P6, If[t > P5, (R6 - R5) * (t - P5) / (P6 - P5) + R5, 0], 0)
S5 = If(t ≤ P5, If[t > P4, (R5 - R4) * (t - P4) / (P5 - P4) + R4, 0], 0)
S4 = If(t ≤ P4, If[t > P3, (R4 - R3) * (t - P3) / (P4 - P3) + R3, 0], 0)
S3 = If(t ≤ P3, If[t > P2, (R3 - R2) * (t - P2) / (P3 - P2) + R2, 0], 0)
S2 = If(t ≤ P2, If[t > P1, (R2 - R1) * (t - P1) / (P2 - P1) + R1, 0], 0)
P69 = 0
S69 = If(t ≤ P69, If[t > P68, (R69 - R68) * (t - P68) / (P69 - P68) + R68, 0], 0)
S100 = Sum(S2:S69)
O3 = If(t < P3, If[t > P2, (P3 - t) / (P3 - P2), 0], 0)
O4 = If(t < P4, If[t > P3, (P4 - t) / (P4 - P3), 0], 0)
O5 = If(t < P5, If[t > P4, (P5 - t) / (P5 - P4), 0], 0)
O6 = If(t < P6, If[t > P5, (P6 - t) / (P6 - P5), 0], 0)
O7 = If(t < P7, If[t > P6, (P7 - t) / (P7 - P6), 0], 0)
O8 = If(t < P8, If[t > P7, (P8 - t) / (P8 - P7), 0], 0)
O9 = If(t < P9, If[t > P8, (P9 - t) / (P9 - P8), 0], 0)
O12 = If(t < P12, If[t > P11, (P12 - t) / (P12 - P11), 0], 0)
O13 = If(t < P13, If[t > P12, (P13 - t) / (P13 - P12), 0], 0)
O16 = If(t < P16, If[t > P15, (P16 - t) / (P16 - P15), 0], 0)
O17 = If(t < P17, If[t > P16, (P17 - t) / (P17 - P16), 0], 0)
O20 = If(t < P20, If[t > P19, (P20 - t) / (P20 - P19), 0], 0)
O21 = If(t < P21, If[t > P20, (P21 - t) / (P21 - P20), 0], 0)
O22 = If(t < P22, If[t > P21, (P22 - t) / (P22 - P21), 0], 0)
O23 = If(t < P23, If[t > P22, (P23 - t) / (P23 - P22), 0], 0)
O24 = If(t < P24, If[t > P23, (P24 - t) / (P24 - P23), 0], 0)
O10 = If(t < P10, If[t > P9, (P10 - t) / (P10 - P9), 0], 0)
O11 = If(t < P11, If[t > P10, (P11 - t) / (P11 - P10), 0], 0)
O14 = If(t < P14, If[t > P13, (P14 - t) / (P14 - P13), 0], 0)
O15 = If(t < P15, If[t > P14, (P15 - t) / (P15 - P14), 0], 0)
O18 = If(t < P18, If[t > P17, (P18 - t) / (P18 - P17), 0], 0)
O19 = If(t < P19, If[t > P18, (P19 - t) / (P19 - P18), 0], 0)
O25 = If(t < P25, If[t > P24, (P25 - t) / (P25 - P24), 0], 0)
O26 = If(t < P26, If[t > P25, (P26 - t) / (P26 - P25), 0], 0)
O27 = If(t < P27, If[t > P26, (P27 - t) / (P27 - P26), 0], 0)
O28 = If(t < P28, If[t > P27, (P28 - t) / (P28 - P27), 0], 0)
O29 = If(t < P29, If[t > P28, (P29 - t) / (P29 - P28), 0], 0)
O30 = If(t < P30, If[t > P29, (P30 - t) / (P30 - P29), 0], 0)
O33 = If(t < P33, If[t > P32, (P33 - t) / (P33 - P32), 0], 0)
O34 = If(t < P34, If[t > P33, (P34 - t) / (P34 - P33), 0], 0)
O35 = If(t < P35, If[t > P34, (P35 - t) / (P35 - P34), 0], 0)
O36 = If(t < P36, If[t > P35, (P36 - t) / (P36 - P35), 0], 0)
O37 = If(t < P37, If[t > P36, (P37 - t) / (P37 - P36), 0], 0)
O38 = If(t < P38, If[t > P37, (P38 - t) / (P38 - P37), 0], 0)
O39 = If(t < P39, If[t > P38, (P39 - t) / (P39 - P38), 0], 0)
O40 = If(t < P40, If[t > P39, (P40 - t) / (P40 - P39), 0], 0)
O41 = If(t < P41, If[t > P40, (P41 - t) / (P41 - P40), 0], 0)
O42 = If(t < P42, If[t > P41, (P42 - t) / (P42 - P41), 0], 0)
O43 = If(t < P43, If[t > P42, (P43 - t) / (P43 - P42), 0], 0)
O44 = If(t < P44, If[t > P43, (P44 - t) / (P44 - P43), 0], 0)
O45 = If(t < P45, If[t > P44, (P45 - t) / (P45 - P44), 0], 0)
O46 = If(t < P46, If[t > P45, (P46 - t) / (P46 - P45), 0], 0)
O47 = If(t < P47, If[t > P46, (P47 - t) / (P47 - P46), 0], 0)
O48 = If(t < P48, If[t > P47, (P48 - t) / (P48 - P47), 0], 0)
O49 = If(t < P49, If[t > P48, (P49 - t) / (P49 - P48), 0], 0)
O50 = If(t < P50, If[t > P49, (P50 - t) / (P50 - P49), 0], 0)
O31 = If(t < P31, If[t > P30, (P31 - t) / (P31 - P30), 0], 0)
O2 = If(t < P2, If[t > P1, (P2 - t) / (P2 - P1), 0], 0)
O51 = If(t < P51, If[t > P50, (P51 - t) / (P51 - P50), 0], 0)
O52 = If(t < P52, If[t > P51, (P52 - t) / (P52 - P51), 0], 0)
O53 = If(t < P53, If[t > P52, (P53 - t) / (P53 - P52), 0], 0)
O54 = If(t < P54, If[t > P53, (P54 - t) / (P54 - P53), 0], 0)
O55 = If(t < P55, If[t > P54, (P55 - t) / (P55 - P54), 0], 0)
O56 = If(t < P56, If[t > P55, (P56 - t) / (P56 - P55), 0], 0)
O57 = If(t < P57, If[t > P56, (P57 - t) / (P57 - P56), 0], 0)
O58 = If(t < P58, If[t > P57, (P58 - t) / (P58 - P57), 0], 0)
O59 = If(t < P59, If[t > P58, (P59 - t) / (P59 - P58), 0], 0)
O60 = If(t < P60, If[t > P59, (P60 - t) / (P60 - P59), 0], 0)
O61 = If(t < P61, If[t > P60, (P61 - t) / (P61 - P60), 0], 0)
O62 = If(t < P62, If[t > P61, (P62 - t) / (P62 - P61), 0], 0)
O63 = If(t < P63, If[t > P62, (P63 - t) / (P63 - P62), 0], 0)
O64 = If(t < P64, If[t > P63, (P64 - t) / (P64 - P63), 0], 0)
O65 = If(t < P65, If[t > P64, (P65 - t) / (P65 - P64), 0], 0)
O66 = If(t < P66, If[t > P65, (P66 - t) / (P66 - P65), 0], 0)
O67 = If(t < P67, If[t > P66, (P67 - t) / (P67 - P66), 0], 0)
O68 = If(t < P68, If[t > P67, (P68 - t) / (P68 - P67), 0], 0)
O32 = If(t < P32, If[t > P31, (P32 - t) / (P32 - P31), 0], 0)
phaseL1: (0.3 - u, -3.2)
phaseL2: phaseL1 + (u, 0)
phaseL4 = Midpoint(phaseL1, phaseL2)
phaseL7: phaseL4 + (0, u)
phaseL3: phaseL1 + (0, u)
phaseL8: phaseL3 + (u + 0.2, 0)
phaseL6: phaseL1 + (u + 0.2, 0)
R_S = 0.6 (slider from 0.1 to 1.5, step 0.05)
D_{TS} = 2.9 (slider from 0.5 to 5, step 0.1)
S: D + (0.3, 9.5) - ((x(F) - x(D))^(3), (y(F) - y(D))^(3)) * 0.1
Soleil1 = Circle(S, R_S)
theta = 3.14159 (slider from -6.28 to 6.28, step 0.1)
phi = 1 (slider from 0 to 4.712388, step 0.1)
Dplc: (theta, phi)
C13 = 3.141592653589793 rad
C12 = 1.5707963267948966 rad
C10 = 0 rad
C14 = 4.71238898038469 rad
B4 = 0.017
B8 = 1.017
B9: B8 * (1 - B4)
B13: B9 / (1 + B4 * cos(C13))
D13: B13 * cos(C13)
E13: B13 * sin(C13)
T4: Point3D[Dplc, D_{TS} * D13, D_{TS} * E13, 0] + S
B12: B9 / (1 + B4 * cos(C12))
D12: B12 * cos(C12)
E12: B12 * sin(C12)
T3: Point3D[Dplc, D_{TS} * D12, D_{TS} * E12, 0] + S
B10: B9 / (1 + B4 * cos(C10))
E10: B10 * sin(C10)
D10: B10 * cos(C10)
T1: Point3D[Dplc, D_{TS} * D10, D_{TS} * E10, 0] + S
B14: B9 / (1 + B4 * cos(C14))
D14: B14 * cos(C14)
E14: B14 * sin(C14)
T5: Point3D[Dplc, D_{TS} * D14, D_{TS} * E14, 0] + S
C11 = 0.7853981633974483 rad
B11: B9 / (1 + B4 * cos(C11))
D11: B11 * cos(C11)
E11: B11 * sin(C11)
T2: Point3D[Dplc, D_{TS} * D11, D_{TS} * E11, 0] + S
Trajectoire_T = Conic(T5, T1, T2, T3, T4)
ArcTrajectoire1_1 = Arc(Trajectoire_T, T1, T4)
temps: {967, 10, 55, 14, 9, 3, 2012, "March", "Friday", 6}
inclinaison_L = 5.145 (slider from -5.145 to 5.145, step 0.1)
F11 = 0.7853981633974483 rad
G4 = 0.0549
G8 = 0.00271
G9: G8 * (1 - G4)
G11: G9 / (1 + G4 * cos(C11))
B5 = 365.25636
Ajustement_T = 176.058 (slider from 0 to 360, step 0.0010)
C15: (t * 360 / B5)° + Ajustement_T°
B15: B9 / (1 + B4 * cos(C15))
Date9: 2455927.5 + t
Date2: floor((floor(Date9 + 0.5) - 1867216.25) / 36524.25)
Date1: floor(Date9 + 0.5) + 1 + Date2 - floor(Date2 / 4)
Date3: Date1 + 1524
Date4: floor((Date3 - 122.1) / 365.25)
Date5: floor(365.25 * Date4)
Date7: (Date9 + 0.5 - floor(Date9 + 0.5)) * 24
Date8: floor(Date7)
Date13: floor((Date7 - Date8) * 60)
Date10 = Mod(floor(Date9 - 0.5) + 0.5 + 1.5, 7)
Date6: floor((Date3 - Date5) / 30.6001)
Date11: Date3 - Date5 - floor(30.6001 * Date6) + Mod[Date9 + 0.5, 1]
Date12: Date6 - If[Date6 < 13.5, 1, 13]
Date14: Date4 - If[Date12 > 2.5, 4716, 4715]
B3 = 1.49597870691E8
B7 = 0.983
E15: B15 * sin(C15)
D15: B15 * cos(C15)
i: (1, 0) * cos(theta) + (0, 1) * sin(theta) * cos(phi)
j: -(1, 0) * sin(theta) + (0, 1) * cos(theta) * cos(phi)
k: (0, 1) * sin(phi)
ArcTrajectoire1 = Arcd'éllipse(Trajectoire_T, T4, T1)
ArcTrajectoire2_1 = Arcd'éllipse(Trajectoire_T, T4, T1)
ArcTrajectoire2 = Arc(Trajectoire_T, T1, T4)
T: Point3D[Dplc, D_{TS} * D15, D_{TS} * E15, 0] + S
v = 0.06
Trajectoire_{T1} = Segment(T1, T4)
G7 = 0.00243
lune_P = 27.149999999999977 (slider from 26 to 29, step 0.01)
G5: lune_P
G10: G9 / (1 + G4 * cos(C10))
G12: G9 / (1 + G4 * cos(C12))
G13: G9 / (1 + G4 * cos(C13))
G14: G9 / (1 + G4 * cos(C14))
I10: G10 * sin(C10)
I12: G12 * sin(C12)
I13: G13 * sin(C13)
I14: G14 * sin(C14)
G3: inclinaison_L°
H13: G13 * cos(C13) * cos(G3)
H10_1: G10 * cos(C10) * cos(G3)
H12_1: G12 * cos(C12)
H14: G14 * cos(C14)
J12 = 0
J14 = 0
J10: G10 * cos(C10) * sin(G3)
J13: G13 * cos(C13) * sin(G3)
I11: G11 * sin(F11) * cos(G3)
J11: G11 * cos(F11) * sin(G3)
H11_1: G11 * cos(F11) * cos(G3)
D_{TL} = 500 (slider from 100 to 900, step 10)
L1: Point3D[Dplc, D_{TL} * H10_1, D_{TL} * I10, D_{TL} * J10] + T
L2: Point3D[Dplc, D_{TL} * H11_1, D_{TL} * I11, D_{TL} * J11] + T
L3: Point3D[Dplc, D_{TL} * H12_1, D_{TL} * I12, D_{TL} * J12] + T
L4: Point3D[Dplc, D_{TL} * H13, D_{TL} * I13, D_{TL} * J13] + T
L5: Point3D[Dplc, D_{TL} * H14, D_{TL} * I14, D_{TL} * J14] + T
Trajectoire_{L1} = Conic(L4, L3, L2, L1, L5)
Trajectoire_1 = Segment(L1, L4)
R_T = 0.5 (slider from 0.1 to 1.5, step 0.1)
Terre = Circle(T, R_T)
Ajustement_L: S100
C16: (t * 360 / G5)° + Ajustement_L°
G15: G9 / (1 + G4 * cos(C16))
H15: G15 * cos(C16) * cos(G3)
I15: G15 * sin(C16) * cos(G3)
J15: G15 * cos(C16) * sin(G3)
L: Point3D[Dplc, D_{TL} * H15, D_{TL} * I15, D_{TL} * J15] + T
R_M = 0.11 (slider from 0.1 to 1, step 0.01)
Lune1 = Circle(L, R_M)
b = 81.6527
c = 81.6527
d = 81.6527
e = 81.6527
h = 81.6527
l = 583.994
m = 583.994
n = 583.994
o = 583.994
tt: t + 15
F17 = IsInteger(floor(2 * tt / 365.25) / 2)
P_{NORD}: Point3D[Dplc, -0.97 * R_T * cos(66.5°), 0, 0.97 * R_T * sin(66.5°)] + T
P_{SUD}: Point3D[Dplc, 0.97 * R_T * cos(66.5°), 0, -0.97 * R_T * sin(66.5°)] + T
Equateur1 = Conic(Point3D[Dplc, 0, -R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(23.5°), 0, R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(°) * cos(23.5°), R_T * sin(°) * cos(23.5°), R_T * sin(23.5°)] + T, Point3D[Dplc, 0, R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, -R_T * cos(23.5°), 0, -R_T * sin(23.5°)] + T)
P_{NORD1}: Point3D[Dplc, -1.4 * R_T * cos(66.5°), 0, 1.4 * R_T * sin(66.5°)] + T
axeN1 = Segment(P_{NORD}, P_{NORD1})
P_{SUD1}: Point3D[Dplc, 1.4 * R_T * cos(66.5°), 0, -1.4 * R_T * sin(66.5°)] + T
AxeS1 = Segment(P_{SUD1}, P_{SUD})
Eq = Segment(Point3D[Dplc, R_T * cos(23.5°), 0, R_T * sin(23.5°)] + T, Point3D[Dplc, -R_T * cos(23.5°), 0, -R_T * sin(23.5°)] + T)
Eq1: Point3D[Dplc, -R_T * cos(23.5°), 0, -R_T * sin(23.5°)] + T
Eq2: Point3D[Dplc, R_T * cos(23.5°), 0, R_T * sin(23.5°)] + T
P_{Cancer1}: Point3D[Dplc, R_T * cos(47°), 0, R_T * sin(47°)] + T
P_{Cancer2}: Point3D[Dplc, -R_T, 0, R_T * 0] + T
EqC1 = Midpoint(P_{Cancer1}, P_{Cancer2})
Cancer1 = Conic(Point3D[Dplc, 0, -cos(23.5°) * R_T * cos(23.5°), 0] + EqC1, Point3D[Dplc, cos(23.5°) * R_T * cos(23.5°), 0, cos(23.5°) * R_T * sin(23.5°)] + EqC1, Point3D[Dplc, cos(23.5°) * R_T * cos(°) * cos(23.5°), cos(23.5°) * R_T * sin(°) * cos(23.5°), cos(23.5°) * R_T * sin(23.5°)] + EqC1, Point3D[Dplc, 0, cos(23.5°) * R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + EqC1, Point3D[Dplc, -cos(23.5°) * R_T * cos(23.5°), 0, -cos(23.5°) * R_T * sin(23.5°)] + EqC1)
EqC2 = Mirror(EqC1, T)
Capricorne1 = Conic(Point3D[Dplc, 0, -cos(23.5°) * R_T * cos(23.5°), 0] + EqC2, Point3D[Dplc, cos(23.5°) * R_T * cos(23.5°), 0, cos(23.5°) * R_T * sin(23.5°)] + EqC2, Point3D[Dplc, cos(23.5°) * R_T * cos(°) * cos(23.5°), cos(23.5°) * R_T * sin(°) * cos(23.5°), cos(23.5°) * R_T * sin(23.5°)] + EqC2, Point3D[Dplc, 0, cos(23.5°) * R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + EqC2, Point3D[Dplc, -cos(23.5°) * R_T * cos(23.5°), 0, -cos(23.5°) * R_T * sin(23.5°)] + EqC2)
Equateur1_1 = Conic(Point3D[Dplc, 0, -R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(23.5°), 0, R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(°) * cos(23.5°), R_T * sin(°) * cos(23.5°), R_T * sin(23.5°)] + T, Point3D[Dplc, 0, R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, -R_T * cos(23.5°), 0, -R_T * sin(23.5°)] + T)
axeN2 = Segment(P_{NORD}, P_{NORD1})
Equateur2 = Conic(Point3D[Dplc, 0, -R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(23.5°), 0, R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(°) * cos(23.5°), R_T * sin(°) * cos(23.5°), R_T * sin(23.5°)] + T, Point3D[Dplc, 0, R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, -R_T * cos(23.5°), 0, -R_T * sin(23.5°)] + T)
AxeS2 = Segment(P_{SUD1}, P_{SUD})
Cancer2 = Conic(Point3D[Dplc, 0, -cos(23.5°) * R_T * cos(23.5°), 0] + EqC1, Point3D[Dplc, cos(23.5°) * R_T * cos(23.5°), 0, cos(23.5°) * R_T * sin(23.5°)] + EqC1, Point3D[Dplc, cos(23.5°) * R_T * cos(°) * cos(23.5°), cos(23.5°) * R_T * sin(°) * cos(23.5°), cos(23.5°) * R_T * sin(23.5°)] + EqC1, Point3D[Dplc, 0, cos(23.5°) * R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + EqC1, Point3D[Dplc, -cos(23.5°) * R_T * cos(23.5°), 0, -cos(23.5°) * R_T * sin(23.5°)] + EqC1)
CancerS = Segment(P_{Cancer1}, P_{Cancer2})
Capri1 = Mirror(P_{Cancer1}, Eq)
Capri2 = Mirror(P_{Cancer2}, Eq)
CapricorneS = Segment(Capri1, Capri2)
ArcTrajectoire2_2 = Arc(Trajectoire_T, T4, T1)
ArcTrajectoire1_2 = Arc(Trajectoire_T, T1, T4)
Soleil2 = Circle(S, R_S)
Soleil3 = Circle(S, R_S)
Soleil4 = Circle(S, R_S)
Lune2 = Circle(L, R_M)
F18: floor(C16 / 180) / 2
texte9: "" + (Element[temps, 9]) + " " + (Element[temps, 5]) + " " + (Element[temps, 8]) + " " + (Element[temps, 7]) + ""
Equat = Conic(Point3D[Dplc, 0, -R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(23.5°), 0, R_T * sin(23.5°)] + T, Point3D[Dplc, R_T * cos(°) * cos(23.5°), R_T * sin(°) * cos(23.5°), R_T * sin(23.5°)] + T, Point3D[Dplc, 0, R_T * cos(23.5°), 0 * R_T * sin(23.5°)] + T, Point3D[Dplc, -R_T * cos(23.5°), 0, -R_T * sin(23.5°)] + T)
Lune3 = Circle(L, R_M)
Lune4 = Circle(L, R_M)
F19 = IsInteger(floor(10 * t / 180) / 2)
F20 = IsInteger(floor(10 * (t + 10) / 180) / 2)
F21 = IsInteger(floor(10 * (t + 20) / 180) / 2)
F22 = IsInteger(floor(10 * (t + 30) / 180) / 2)
F23 = IsInteger(floor(10 * (t + 40) / 180) / 2)
F24 = IsInteger(floor(10 * (t + 50) / 180) / 2)
F25 = IsInteger(floor(10 * (t + 60) / 180) / 2)
F26 = IsInteger(floor(10 * (t + 70) / 180) / 2)
F27 = IsInteger(floor(10 * (t + 80) / 180) / 2)
F28 = IsInteger(floor(10 * (t + 90) / 180) / 2)
F29 = IsInteger(floor(10 * (t + 100) / 180) / 2)
F30 = IsInteger(floor(10 * (t + 110) / 180) / 2)
F31 = IsInteger(floor(10 * (t + 120) / 180) / 2)
F32 = IsInteger(floor(10 * (t + 130) / 180) / 2)
F33 = IsInteger(floor(10 * (t + 140) / 180) / 2)
F34 = IsInteger(floor(10 * (t + 150) / 180) / 2)
F35 = IsInteger(floor(10 * (t + 160) / 180) / 2)
F36 = IsInteger(floor(10 * (t + 170) / 180) / 2)
F37 = IsInteger(floor(10 * (t + 180) / 180) / 2)
aa: sgn(sin(phi))
A0: Point3D[Dplc, R_T * cos(23.5°), 0, R_T * sin(23.5°)] + T
A: Point3D[Dplc, -R_T * cos(23.5°), 0, -R_T * sin(23.5°)] + T
f = Arc(Equateur1_1, A0, A)
g = Arc(Equateur1_1, A, A0)
p = Arc(Equateur1_1, A, A0)
q = Arc(Equateur1_1, A0, A)
r = Arc(Capricorne1, Capri1, Capri2)
s = Arc(Capricorne1, Capri1, Capri2)
c_1 = Arc(Capricorne1, Capri2, Capri1)
d_1 = Arc(Capricorne1, Capri2, Capri1)
e_1 = Arc(Cancer1, P_{Cancer1}, P_{Cancer2})
f_1 = Arc(Cancer1, P_{Cancer1}, P_{Cancer2})
g_1 = Arc(Cancer1, P_{Cancer2}, P_{Cancer1})
h_1 = Arc(Cancer1, P_{Cancer2}, P_{Cancer1})
logo = 2 (slider from 0 to 2, step 0.0010)
o_1 = false
Trajectoire_{L2} = Conic(L4, L3, L2, L1, L5)
poly1, a_11, a_12, a_13, a_14 = Polygon((-199.82011506398368, 199.19910680862094), (200.17988493601632, 199.19910680862094), (200.17988493601632, -200.80089319137906), (-199.82011506398368, -200.80089319137906))
texte1 = RotateText("Sun-Earth", 270*°)
texte2 = RotateText("Earth-Moon", 270*°)
texte3 = RotateText("rotation", 270*°)
trajL4 = Arcd'éllipse(Trajectoire_{L2}, L1, L4)
Lune5 = Circle(L, R_M)
Lune6 = Circle(L, R_M)
Soleil5 = Circle(S, R_S)
Soleil6 = Circle(S, R_S)
trajL3 = Arcd'éllipse(Trajectoire_{L2}, L4, L1)
trajL1 = Arcd'éllipse(Trajectoire_{L2}, L1, L4)
trajL2 = Arcd'éllipse(Trajectoire_{L2}, L4, L1)
trajL5 = Arcd'éllipse(Trajectoire_{L2}, L1, L4)
trajL6 = Arcd'éllipse(Trajectoire_{L2}, L4, L1)
vitesselogo = If(logo < 1, 1, 10)
PS_{NORD}: Point3D[Dplc, 0, 0, R_S] + S
PS_{SUD}: Point3D[Dplc, 0, 0, -R_S] + S
PS_{NORD1}: Point3D[Dplc, 0, 0, 0.2 + R_S] + S
axesoleil1 = Segment(PS_{NORD}, PS_{NORD1})
PS_{SUD1}: Point3D[Dplc, 0, 0, -0.2 - R_S] + S
axesoleil2 = Segment(PS_{SUD1}, PS_{SUD})
axesoleil11 = Segment(PS_{NORD}, PS_{NORD1})
axesoleil22 = Segment(PS_{SUD1}, PS_{SUD})
A_2: Point3D[Dplc, aa * R_T * cos((10 * (t + 10))°) * cos(23.5°), aa * R_T * sin((10 * (t + 10))°) * cos(23.5°), aa * R_T * cos((10 * (t + 10))°) * sin(23.5°)] + T
A_1: Point3D[Dplc, aa * R_T * cos((20 * t)°) * cos(23.5°), aa * R_T * sin((20 * t)°) * cos(23.5°), aa * R_T * cos((20 * t)°) * sin(23.5°)] + T
A_3: Point3D[Dplc, aa * R_T * cos((10 * (t + 20))°) * cos(23.5°), aa * R_T * sin((10 * (t + 20))°) * cos(23.5°), aa * R_T * cos((10 * (t + 20))°) * sin(23.5°)] + T
A_4: Point3D[Dplc, aa * R_T * cos((10 * (t + 30))°) * cos(23.5°), aa * R_T * sin((10 * (t + 30))°) * cos(23.5°), aa * R_T * cos((10 * (t + 30))°) * sin(23.5°)] + T
A_5: Point3D[Dplc, aa * R_T * cos((10 * (t + 40))°) * cos(23.5°), aa * R_T * sin((10 * (t + 40))°) * cos(23.5°), aa * R_T * cos((10 * (t + 40))°) * sin(23.5°)] + T
A_6: Point3D[Dplc, aa * R_T * cos((10 * (t + 50))°) * cos(23.5°), aa * R_T * sin((10 * (t + 50))°) * cos(23.5°), aa * R_T * cos((10 * (t + 50))°) * sin(23.5°)] + T
A_7: Point3D[Dplc, aa * R_T * cos((10 * (t + 60))°) * cos(23.5°), aa * R_T * sin((10 * (t + 60))°) * cos(23.5°), aa * R_T * cos((10 * (t + 60))°) * sin(23.5°)] + T
A_8: Point3D[Dplc, aa * R_T * cos((10 * (t + 70))°) * cos(23.5°), aa * R_T * sin((10 * (t + 70))°) * cos(23.5°), aa * R_T * cos((10 * (t + 70))°) * sin(23.5°)] + T
A_9: Point3D[Dplc, aa * R_T * cos((10 * (t + 80))°) * cos(23.5°), aa * R_T * sin((10 * (t + 80))°) * cos(23.5°), aa * R_T * cos((10 * (t + 80))°) * sin(23.5°)] + T
A_{10}: Point3D[Dplc, aa * R_T * cos((10 * (t + 90))°) * cos(23.5°), aa * R_T * sin((10 * (t + 90))°) * cos(23.5°), aa * R_T * cos((10 * (t + 90))°) * sin(23.5°)] + T
A_{11}: Point3D[Dplc, aa * R_T * cos((10 * (t + 100))°) * cos(23.5°), aa * R_T * sin((10 * (t + 100))°) * cos(23.5°), aa * R_T * cos((10 * (t + 100))°) * sin(23.5°)] + T
A_{12}: Point3D[Dplc, aa * R_T * cos((10 * (t + 110))°) * cos(23.5°), aa * R_T * sin((10 * (t + 110))°) * cos(23.5°), aa * R_T * cos((10 * (t + 110))°) * sin(23.5°)] + T
A_{13}: Point3D[Dplc, aa * R_T * cos((10 * (t + 120))°) * cos(23.5°), aa * R_T * sin((10 * (t + 120))°) * cos(23.5°), aa * R_T * cos((10 * (t + 120))°) * sin(23.5°)] + T
A_{14}: Point3D[Dplc, aa * R_T * cos((10 * (t + 130))°) * cos(23.5°), aa * R_T * sin((10 * (t + 130))°) * cos(23.5°), aa * R_T * cos((10 * (t + 130))°) * sin(23.5°)] + T
A_{15}: Point3D[Dplc, aa * R_T * cos((10 * (t + 140))°) * cos(23.5°), aa * R_T * sin((10 * (t + 140))°) * cos(23.5°), aa * R_T * cos((10 * (t + 140))°) * sin(23.5°)] + T
A_{16}: Point3D[Dplc, aa * R_T * cos((10 * (t + 150))°) * cos(23.5°), aa * R_T * sin((10 * (t + 150))°) * cos(23.5°), aa * R_T * cos((10 * (t + 150))°) * sin(23.5°)] + T
A_{17}: Point3D[Dplc, aa * R_T * cos((10 * (t + 160))°) * cos(23.5°), aa * R_T * sin((10 * (t + 160))°) * cos(23.5°), aa * R_T * cos((10 * (t + 160))°) * sin(23.5°)] + T
A_{18}: Point3D[Dplc, aa * R_T * cos((10 * (t + 170))°) * cos(23.5°), aa * R_T * sin((10 * (t + 170))°) * cos(23.5°), aa * R_T * cos((10 * (t + 170))°) * sin(23.5°)] + T
A_{19}: Point3D[Dplc, aa * R_T * cos((10 * (t + 180))°) * cos(23.5°), aa * R_T * sin((10 * (t + 180))°) * cos(23.5°), aa * R_T * cos((10 * (t + 180))°) * sin(23.5°)] + T
SS: Point3D[Dplc, 0, 0, D_{TL} * J15] + S
z = Vector(S, T)
a_1 = Line(SS, z)
Pleine_L, Nouvelle_L = Intersect(Trajectoire_{L1}, a_1)
PO51 = 0
O69 = If(t < P69, If[t > P68, (-(P68 - t)) / (P69 - P68), 0], 0)
O100 = Sum(O2:O69)
mvt: phaseL1 + (u * O100, u / 2)
Date: "" + (If[Date10 ≟ 0, "Sunday", If[Date10 ≟ 1, "Monday", If[Date10 ≟ 2, "Tuesday", If[Date10 ≟ 3, "Wednesday", If[Date10 ≟ 4, "Thursday", If[Date10 ≟ 5, "Friday", "Saturday"]]]]]]) + " " + (floor(Date11)) + "/" + Date12 + "/" + Date14 + " - " + Date8 + " h " + Date13 + " min"
N2 = If(P2 ≟ 0, 1, If[t < P2, 1, 0])
N3 = If(P3 ≟ 0, 1, If[t < P3, 1, 0])
N4 = If(P4 ≟ 0, 1, If[t < P4, 1, 0])
N5 = If(P5 ≟ 0, 1, If[t < P5, 1, 0])
N6 = If(P6 ≟ 0, 1, If[t < P6, 1, 0])
N7 = If(P7 ≟ 0, 1, If[t < P7, 1, 0])
N8 = If(P8 ≟ 0, 1, If[t < P8, 1, 0])
N9 = If(P9 ≟ 0, 1, If[t < P9, 1, 0])
N10 = If(P10 ≟ 0, 1, If[t < P10, 1, 0])
N11 = If(P11 ≟ 0, 1, If[t < P11, 1, 0])
N12 = If(P12 ≟ 0, 1, If[t < P12, 1, 0])
N13 = If(P13 ≟ 0, 1, If[t < P13, 1, 0])
N14 = If(P14 ≟ 0, 1, If[t < P14, 1, 0])
N15 = If(P15 ≟ 0, 1, If[t < P15, 1, 0])
N16 = If(P16 ≟ 0, 1, If[t < P16, 1, 0])
N17 = If(P17 ≟ 0, 1, If[t < P17, 1, 0])
N18 = If(P18 ≟ 0, 1, If[t < P18, 1, 0])
N19 = If(P19 ≟ 0, 1, If[t < P19, 1, 0])
N20 = If(P20 ≟ 0, 1, If[t < P20, 1, 0])
N21 = If(P21 ≟ 0, 1, If[t < P21, 1, 0])
N22 = If(P22 ≟ 0, 1, If[t < P22, 1, 0])
N23 = If(P23 ≟ 0, 1, If[t < P23, 1, 0])
N24 = If(P24 ≟ 0, 1, If[t < P24, 1, 0])
N25 = If(P25 ≟ 0, 1, If[t < P25, 1, 0])
N26 = If(P26 ≟ 0, 1, If[t < P26, 1, 0])
N27 = If(P27 ≟ 0, 1, If[t < P27, 1, 0])
N28 = If(P28 ≟ 0, 1, If[t < P28, 1, 0])
N29 = If(P29 ≟ 0, 1, If[t < P29, 1, 0])
N30 = If(P30 ≟ 0, 1, If[t < P30, 1, 0])
N31 = If(P31 ≟ 0, 1, If[t < P31, 1, 0])
N32 = If(P32 ≟ 0, 1, If[t < P32, 1, 0])
N33 = If(P33 ≟ 0, 1, If[t < P33, 1, 0])
N34 = If(P34 ≟ 0, 1, If[t < P34, 1, 0])
N35 = If(P35 ≟ 0, 1, If[t < P35, 1, 0])
N36 = If(P36 ≟ 0, 1, If[t < P36, 1, 0])
N37 = If(P37 ≟ 0, 1, If[t < P37, 1, 0])
N38 = If(P38 ≟ 0, 1, If[t < P38, 1, 0])
N39 = If(P39 ≟ 0, 1, If[t < P39, 1, 0])
N40 = If(P40 ≟ 0, 1, If[t < P40, 1, 0])
N41 = If(P41 ≟ 0, 1, If[t < P41, 1, 0])
N42 = If(P42 ≟ 0, 1, If[t < P42, 1, 0])
N43 = If(P43 ≟ 0, 1, If[t < P43, 1, 0])
N44 = If(P44 ≟ 0, 1, If[t < P44, 1, 0])
N45 = If(P45 ≟ 0, 1, If[t < P45, 1, 0])
N46 = If(P46 ≟ 0, 1, If[t < P46, 1, 0])
N47 = If(P47 ≟ 0, 1, If[t < P47, 1, 0])
N48 = If(P48 ≟ 0, 1, If[t < P48, 1, 0])
N49 = If(P49 ≟ 0, 1, If[t < P49, 1, 0])
N50 = If(P50 ≟ 0, 1, If[t < P50, 1, 0])
N51 = If(P51 ≟ 0, 1, If[t < P51, 1, 0])
N52 = If(P52 ≟ 0, 1, If[t < P52, 1, 0])
N53 = If(P53 ≟ 0, 1, If[t < P53, 1, 0])
N54 = If(P54 ≟ 0, 1, If[t < P54, 1, 0])
N55 = If(P55 ≟ 0, 1, If[t < P55, 1, 0])
N56 = If(P56 ≟ 0, 1, If[t < P56, 1, 0])
N57 = If(P57 ≟ 0, 1, If[t < P57, 1, 0])
N58 = If(P58 ≟ 0, 1, If[t < P58, 1, 0])
N59 = If(P59 ≟ 0, 1, If[t < P59, 1, 0])
N60 = If(P60 ≟ 0, 1, If[t < P60, 1, 0])
N61 = If(P61 ≟ 0, 1, If[t < P61, 1, 0])
N62 = If(P62 ≟ 0, 1, If[t < P62, 1, 0])
N63 = If(P63 ≟ 0, 1, If[t < P63, 1, 0])
N64 = If(P64 ≟ 0, 1, If[t < P64, 1, 0])
N65 = If(P65 ≟ 0, 1, If[t < P65, 1, 0])
N66 = If(P66 ≟ 0, 1, If[t < P66, 1, 0])
N67 = If(P67 ≟ 0, 1, If[t < P67, 1, 0])
N68 = If(P68 ≟ 0, 1, If[t < P68, 1, 0])
N69 = If(P69 ≟ 0, 1, If[t < P69, 1, 0])
N100 = IsInteger((68 - Sum[N2:N69]) / 2)
Tfixe: E + (14, 7) - ((x(G) - x(E))^(5), (y(G) - y(E))^(5))
Rt2 = 0.65 (slider from 0.4 to 1.5, step 0.05)
Terre_{fixe} = Circle(Tfixe, Rt2)
G17: floor((t + 11) / 365.25)
H20: t + 11 - G17 * 365.25
H21: 0.2569 * H20 - 23.467
H22: -0.258 * H20 + 70.79
H23 = If(H20 ≤ 183, H21, H22)
Sud = Rotate(Tfixe + (0, -Rt2 - 0.15), H23°, Tfixe)
Nord = Mirror(Sud, Tfixe)
l_2 = Segment(Nord, Sud)
Sud_{fixe1}: Tfixe + (0, -Rt2)
Nord_{fixe1} = Mirror(Sud_{fixe1}, Tfixe)
nuit = Semicircle(Nord_{fixe1}, Sud_{fixe1})
jour = Semicircle(Sud_{fixe1}, Nord_{fixe1})
t_2 = OrthogonalLine(Tfixe, l_2)
R = Intersect(Terre_{fixe}, jour)
Eq_{fixe1}, Eq_{fixe2} = Intersect(Terre_{fixe}, t_2)
Equateur_F = Segment(Eq_{fixe1}, Eq_{fixe2})
Tropique_F2 = Segment(Rotate[Eq_{fixe1}, 23.5*°, Tfixe], Mirror[Rotate[Eq_{fixe1}, 23.5*°, Tfixe], l_2])
Tropique_F1 = Mirror(Tropique_F2, Equateur_F)
k_1 = Segment(Nord_{fixe1}, Sud_{fixe1})
Nord' = Mirror(Nord, k_1)
Sud' = Mirror(Sud, k_1)
q_2 = Segment(Sud', Nord')
r_2 = OrthogonalLine(Tfixe, q_2)
Eq_{fixe4}, Eq_{fixe3} = Intersect(Terre_{fixe}, r_2)
s_2 = Segment(Rotate[Eq_{fixe3}, 336.5*°, Tfixe], Mirror[Rotate[Eq_{fixe3}, 336.5*°, Tfixe], q_2])
a_3 = Segment(Eq_{fixe3}, Eq_{fixe4})
i_3 = Mirror(s_2, a_3)
Rayon1: Tfixe + (-Rt2 - 1, 0)
Rayon2: Tfixe + (-Rt2 - 1, Rt2 / 4)
Rayon3: Tfixe + (-Rt2 - 1, Rt2 / 2)
Rayon_4: Tfixe + (-Rt2 - 1, 3 * Rt2 / 4)
Rayon2': Tfixe + (-Rt2 - 1, (-Rt2) / 4)
Rayon3': Tfixe + (-Rt2 - 1, (-Rt2) / 2)
Rayon4': Tfixe + (-Rt2 - 1, -3 * Rt2 / 4)
j_3 = OrthogonalLine(Rayon_4, k_1)
k_3 = Line(Rayon3, j_3)
l_3 = Line(Rayon2, k_3)
m_3 = Line(Rayon1, l_3)
n_3 = Line(Rayon2', m_3)
p_3 = Line(Rayon3', n_3)
q_3 = Line(Rayon4', p_3)
Q = Intersect(jour, j_3)
U = Intersect(jour, k_3)
V = Intersect(jour, l_3)
W = Intersect(jour, m_3)
Z = Intersect(jour, n_3)
B_1 = Intersect(jour, p_3)
C_1 = Intersect(jour, q_3)
r_3 = Segment(Rayon_4, Q)
s_3 = Segment(Rayon3, U)
t_3 = Segment(Rayon2, V)
a_4 = Segment(W, Rayon1)
b_4 = Segment(Rayon2', Z)
c_4 = Segment(B_1, Rayon3')
d_4 = Segment(Rayon4', C_1)
D_1 = Midpoint(r_3)
e_4 = OrthogonalLine(D_1, r_3)
E_1 = Intersect(e_4, s_3)
F_1 = Intersect(e_4, t_3)
G_1 = Intersect(e_4, a_4)
H_1 = Intersect(e_4, b_4)
I_1 = Intersect(e_4, c_4)
J_1 = Intersect(e_4, d_4)
w = Vector(Rayon_4, D_1)
u_1 = Vector(Rayon3, E_1)
v_1 = Vector(Rayon2, F_1)
w_1 = Vector(Rayon1, G_1)
z_1 = Vector(Rayon2', H_1)
u_2 = Vector(Rayon3', I_1)
v_2 = Vector(Rayon4', J_1)
test = If(phi < 2.9, 1, 0)
Mvt0: test * (Eq_{fixe1} + (Eq_{fixe1} - Eq_{fixe2}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 0 * R_T / 6 + floor(t - 0 * R_T / 6)]) + (1 - test) * (Eq_{fixe4} + (Eq_{fixe4} - Eq_{fixe3}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 0 * R_T / 6 + floor(t - 0 * R_T / 6)])
Mvt4: test * (Eq_{fixe1} + (Eq_{fixe1} - Eq_{fixe2}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 4 * Rt2 / 6 + floor(t - 4 * Rt2 / 6)]) + (1 - test) * (Eq_{fixe4} + (Eq_{fixe4} - Eq_{fixe3}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 4 * Rt2 / 6 + floor(t - 4 * Rt2 / 6)])
Mvt5: test * (Eq_{fixe1} + (Eq_{fixe1} - Eq_{fixe2}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 5 * Rt2 / 6 + floor(t - 5 * Rt2 / 6)]) + (1 - test) * (Eq_{fixe4} + (Eq_{fixe4} - Eq_{fixe3}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 5 * Rt2 / 6 + floor(t - 5 * Rt2 / 6)])
Mvt6: test * (Eq_{fixe1} + (Eq_{fixe1} - Eq_{fixe2}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + Rt2 + floor(t - Rt2)]) + (1 - test) * (Eq_{fixe4} + (Eq_{fixe4} - Eq_{fixe3}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + Rt2 + floor(t - Rt2)])
Mvt1: test * (Eq_{fixe1} + (Eq_{fixe1} - Eq_{fixe2}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + Rt2 / 6 + floor(t - Rt2 / 6)]) + (1 - test) * (Eq_{fixe4} + (Eq_{fixe4} - Eq_{fixe3}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + Rt2 / 6 + floor(t - Rt2 / 6)])
Mvt2: test * (Eq_{fixe1} + (Eq_{fixe1} - Eq_{fixe2}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 2 * Rt2 / 6 + floor(t - 2 * Rt2 / 6)]) + (1 - test) * (Eq_{fixe4} + (Eq_{fixe4} - Eq_{fixe3}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 2 * Rt2 / 6 + floor(t - 2 * Rt2 / 6)])
Mvt3: test * (Eq_{fixe1} + (Eq_{fixe1} - Eq_{fixe2}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 3 * Rt2 / 6 + floor(t - 3 * Rt2 / 6)]) + (1 - test) * (Eq_{fixe4} + (Eq_{fixe4} - Eq_{fixe3}) * If[(floor(t) ≤ t) ∧ (t ≤ 1 + floor(t)), -t + 3 * Rt2 / 6 + floor(t - 3 * Rt2 / 6)])
arc_{CacheL} = CircumcircleArc(phaseL7, mvt, phaseL4)
phaseL8' = Mirror(phaseL8, phaseL7)
phaseL6' = Mirror(phaseL6, phaseL4)
lune1 = Point(arc_{CacheL})
lune2 = Point(arc_{CacheL})
lune3 = Point(arc_{CacheL})
lune4 = Point(arc_{CacheL})
lune5 = Point(arc_{CacheL})
lune6 = Point(arc_{CacheL})
lune7 = Point(arc_{CacheL})
lune8 = Point(arc_{CacheL})
lune9 = Point(arc_{CacheL})
lune10 = Point(arc_{CacheL})
lune11 = Point(arc_{CacheL})
lune12 = Point(arc_{CacheL})
lune13 = Point(arc_{CacheL})
lune14 = Point(arc_{CacheL})
lune15 = Point(arc_{CacheL})
lune16 = Point(arc_{CacheL})
lune17 = Point(arc_{CacheL})
lune18 = Point(arc_{CacheL})
lune19 = Point(arc_{CacheL})
lune20 = Point(arc_{CacheL})
poly2, phasel8, phasel7, lune1_1, lune2_1, lune3_1, lune4_1, lune5_1, lune6_1, lune7_1, lune8_1, lune9_1, lune10_1, mvt_1, lune11_1, lune12_1, lune13_1, lune14_1, lune15_1, lune16_1, lune17_1, lune18_1, lune20_1, phasel4, phasel6 = Polygon(phaseL8, phaseL7, lune1, lune2, lune3, lune4, lune5, lune6, lune7, lune8, lune9, lune10, mvt, lune11, lune12, lune13, lune14, lune15, lune16, lune17, lune18, lune20, phaseL4, phaseL6)
poly3, phasel8', phasel7_1, lune1_2, lune2_2, lune3_2, lune4_2, lune5_2, lune6_2, lune7_2, lune8_2, lune9_2, lune10_2, mvt_2, lune11_2, lune12_2, lune13_2, lune14_2, lune15_2, lune16_2, lune17_2, lune18_2, lune20_2, phasel4_1, phasel6' = Polygon(phaseL8', phaseL7, lune1, lune2, lune3, lune4, lune5, lune6, lune7, lune8, lune9, lune10, mvt, lune11, lune12, lune13, lune14, lune15, lune16, lune17, lune18, lune20, phaseL4, phaseL6')
B: Pleine_L
C: Nouvelle_L
a = false
b_1 = Segment(I, J)
Q100 = 800
LogoP2: (-2.4, -1.5) + If[logo ≥ 1, (logo - 1) * ((2, -3.6) - (-2.4, -1.5)), (0, 0)]
LogoP1: (-7.5, -1.5) + If[logo ≥ 1, (logo - 1) * ((1, -3.6) - (-7.5, -1.5)), (0, 0)]

Tool: Point 3D (command Point3D)
Dplc,x,y,z
Input: Dplc, a, b, c
Output: M

phi: y(Dplc)
theta: x(Dplc)
i: (1, 0) * cos(theta) + (0, 1) * sin(theta) * cos(phi)
j: -(1, 0) * sin(theta) + (0, 1) * cos(theta) * cos(phi)
k: (0, 1) * sin(phi)
M: a * i + b * j + c * k

Tool: Arc d'éllipse (command Arcd'éllipse)
Elipse, Point1,Point2
Input: Trajectoire_T, T1, T4
Output: ArcTrajectoire

Périhélie

Aphélie

ArcTrajectoire = Arc(Trajectoire_T, T1, T4)