gravitation_turtle.py from math import * from turtle import * # valori di G e M massa Terra G=6.67E-11; M=5.97E24 # valori orbita geostazionaria x=4.2164E7; y=0 vx=0; vy=3066 # integrazione numerica t=0; dt=60 while t<86400: setheading(degrees(atan2(vy,vx))) goto(x/1E6, y/1E6) r=sqrt(x**2+y**2) ax=-G*M*x/r**3; ay=-G*M*y/r**3 vx=vx+ax*dt; vy=vy+ay*dt x=x+vx*dt; y=y+vy*dt t=t+dt
https://my.numworks.com/python/gianfranco-oddenino/gravitation_turtle
gravitation_kandinsky.py from math import * from kandinsky import * # valori di G e M massa Terra G=6.67E-11; M=5.97E24 # valori orbita geostazionaria x=4.2164E7; y=0 vx=0; vy=3066 # integrazione numerica t=0; dt=60 fill_rect(159,110,3,3,'blue') while t<86400: set_pixel(160+round(x/1E6),111-round(y/1E6),'black') r=sqrt(x**2+y**2) ax=-G*M*x/r**3; ay=-G*M*y/r**3 vx=vx+ax*dt; vy=vy+ay*dt x=x+vx*dt; y=y+vy*dt t=t+dt
https://my.numworks.com/python/gianfranco-oddenino/gravitation_kandinsky