4 g constante de gravitation
5 k coefficient de frottement
6 dt intervalle de temps de la simulation
8 x0, y0 position initiale
15 fx(t) composante x de la force de frottement à l'instant t
16 fy(t) composante y de la force de frottement à l'instant t
18 ax(t) composante x de l'accélération à l'instant t
19 ay(t) composante y de l'accélération à l'instant t
21 x(t) composante x de la position à l'instant t
22 y(t) composante y de la position à l'instant t
24 vx(t) composante x de la vitesse à l'instant t
25 vy(t) composante y de la vitesse à l'instant t
27 vrx(t) composante x de la vitesse relative à l'instant t
29 vry(t) composante y de la vitesse relative à l'instant t
31 |vr(t)| vitesse relative à l'instant t
32 = SQRT(vrx(t)^2 + vry(t)^2)
49 |f(t)| = -k |vr(t)| vr(t)
50 fx(t) = -k |vr(t)| vrx(t)
51 = -k |vr(t)| (vx(t) - wnd)
52 = -k SQRT((vx(t) - wnd)^2 + vy(t)^2) (vx(t) - wnd)
53 fy(t) = -k vr(t) vry(y)
55 = -k SQRT((vx(t) - wnd)^2 + vy(t)^2) vy(t)
58 = -k SQRT((vx(t) - wnd)^2 + vy(t)^2) (vx(t) - wnd)
60 = -k SQRT((vx(t) - wnd)^2 + vy(t)^2) vy(t) - g
62 x(t+1) = x(t) + vx(t) dt
63 y(t+1) = y(t) + vy(t) dt
65 vx(t+1) = vx(t) + ax(t) dt
66 vy(t+1) = vy(t) + ay(t) dt
90 kvr <- -k × SQRT(vxr × vxr + vy0 × vy0)