Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.mso.anu.edu.au/pfrancis/phys1101/Lectures/L09/iteration_BH_vector.py Дата изменения: Thu Mar 10 02:40:20 2011 Дата индексирования: Tue Oct 2 17:12:56 2012 Кодировка:

from visual import *

scene.background = color.white
scene.range = 200000000

Delta_t = 0.2
Mbh = 1.0E31
G = 6.67E-11

t = 0

spaceship = sphere(pos = vector(100000000,0,0), radius = 5000000, color=color.blue)
spaceship.vel = vector(2000000, 0, 0)

black_hole = sphere(pos = vector(0,0,0), radius = 10000000, color=color.black)

while mag(spaceship.pos) > black_hole.radius:
rate(50)
t = t + Delta_t
spaceship.pos = spaceship.pos + spaceship.vel*Delta_t
Agrav = - G*Mbh*norm(spaceship.pos-black_hole.pos)/(mag(spaceship.pos-
black_hole.pos)**2)
spaceship.vel = spaceship.vel + Agrav*Delta_t