Calcolo di un’orbita astronautica mediante l’integrazione numerica con la calcolatrice grafica
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