Aspect anneaux saturne
Author: phm

Rotation

Déphasage
rot. Terre

Limites dates

Rotation
planètes

point vernal

Anneau

Saturne

Ligne Terre-Saturne

upx = 0.0854788318610717
upy = 0.462606617188416
upz = 0.882433276252181
Ω = 1.9847011589053518 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
i = 0.04336200713409832 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
si: sin(i)
ci: cos(i)
delta = 0 (slider from 0 to 366, step 0.1)
H = (-0.01793834130508549, 0.5136868850878152)
tb = 12000 (slider from 1 to 12000, step 1)
ta = 0 (slider from 0 to 12000, step 0.1)
ta2 = If(ta < tb, ta, tb)
tb2 = If(tb > ta, tb, ta)
j = 0 (slider from 0 to 1, step 0.0010)
jj: (tb2 - ta2) j + ta2
x_T: x(H) + cos(2.0 * 3.14159 * (jj + delta) / 365.25 + 100.575 * 3.14159 / 180.0)
y_T: y(H) + sin(2.0 * 3.14159 * (jj + delta) / 365.25 + 100.575 * 3.14159 / 180.0)
a: 9.537 cos(2.0 * 3.14159 / 10760.0 * (jj + 1453.0))
b: 9.537 sin(2.0 * 3.14159 / 10760.0 * (jj + 1453.0))
y_s: sin(Ω) a + cos(Ω) b
yst: -1.0 * 1.0 * (y_s - y_T)
x_s: (cos(Ω) a - sin(Ω) b) ci
xst: -1.0 * 1.0 * (x_s - x_T)
z_s: si b
zst: -1.0 * 1.0 z_s
dst: sqrt(xst xst + yst yst + zst zst)
uty: yst / dst
utx: xst / dst
utz: zst / dst
thet: acos(upx utx + upy uty + upz utz) 180.0 / 3.141592653589793 - 90.0
thet2: round(thet 1000.0) / 1000.0
S: (x_s, y_s)
a_e: dst / 5.0
b_e: abs(a_e sin(thet 3.141592653589793 / 180.0))
c_e: sqrt(a_e a_e - b_e b_e)
I = (-10.9575, 13)
K = (10.9575, 13)
p: 0.08548 * (1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 x - x(S)) + 0.46244 * (1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 y - y(S)) - 0.88252 z_s = 0.0
T: (x_T, y_T)
c = Circle(H, 1.0)
f: 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 y = y(H)
f' = Rotate(f, 113.715°, H)
O' = (0, 0)
e = Circle(O', 9.537)
N', N = Intersect(e, f')
test: sin(jj + 1453.0)
SH = Segment(H, S)
deltax: x(S) - x(H)
deltay: y(S) - y(H)
année: 2009.0 + jj / 365.25
texte1: "Date : " + année + "
n jour : " + jj
u = Vector(H, (x(H) + 12.0, y(H)))
vthet: thet / 12.0
einc = Ellipse((x(S) - c_e, y(S)), (x(S) + c_e, y(S)), a_e)
d = true
g = true
h = Line(S, T)
k = false
texte5: "Inclinaison : " + thet2 + "°"