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